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Section  1 
INTRODUCTION 


1.1  BACKGROUND 

Increasing  threat  of  global  warfare  has  created  a  significant  need  for  large  space 
platforms  both  for  surveillance  and  for  deployment  of  space  based  defensive  weapon 
systems.  Because  these  systems  may  operate  in  optical  or  mm-wave  regimes,  the 
structures  must  maintain  accurate  geometric  shapes.  The  presence  of  active  systems 
onboard  the  spacecraft  also  may  require  the  minimization  of  vibration  from  random 
disturbances.  Self-defense  may  demand  that  the  platform  be  maneuvered  without 
structural  damage  or  loss  of  performance. 

Insuring  that  such  performance  can  be  met  using  passive  techniques  (e.g.  ,  struc¬ 
tural  stiffening,  material  selection)  generally  poses  unacceptable  weight  penalties  or 
performance  uncertainties  even  when  such  techniques  may  be  feasible.  For  high 
performance  maneuvering,  passive  methods  can  totally  fail  to  provide  adequate  per¬ 
formance.  Active  control  techniques,  however,  may  be  effectively  used  to  create 
so-called  control  configured  spacecraft  where  the  structural  dynamic  properties  in¬ 
cluding  structural  damping  are  strongly  influenced  by  active  closed-loop  control 
systems  placed  on  the  basic  structure.  In  this  way,  requirements  on  passive  stiff¬ 
ness  and  damping  may  be  considerably  relaxed  with  potential  savings  in  weight, 
material  costs,  and  payload  subsystem  complexity.  The  technology  of  active  struc¬ 
tural  control  however  is  still  immature  and  poorly  understood. 

The  performance  potential  for  control  configured  spacecraft  is  the  motivation  for 
the  DARPA/ACOSS  (Active  control  of  Space  Structures)  program  .  This  report, 
covering  the  activities  of  the  LMSC  Phase  I  effort ,  summarizes  the  analytical 
investigations,  control  synthesis  methods,  and  experimental  results  (funded  by 
IRAD)  achieved  during  this  1978  program.  Contributions  of  individual  team  mem¬ 
bers,  including  Lockheed,  Palo  Alto  Research  Lab,  the  LMSC  Space  Systems 
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Division,  Systems  Control,  Inc. ,  Synergistic  Technology,  Inc.  ,  and  Prof. 

R.  E.  Skelton  (Purdue) ,  were  critical  to  the  progress  made  during  Phase  I.  To 
keep  this  report  to  a  reasonable  length,  many  techniques  are  only  summarized 
or  briefly  discussed  as  indicated  in  the  sections  below. 

1.2  OUTLINE  OF  REPORT 

Section  2  is  a  brief  exposition  of  the  technical  approach  taken  in  Phase  I  of 
ACOSS. 

Section  3  deals  with  the  important  subject  of  modeling  techniques  for  flexible 
spacecraft  and  the  relative  merits  and  pitfalls  of  continuum  or  discrete  models 
and  the  effects  of  closed-loop  controls. 

Section  4  comprises  the  bulk  of  the  control  synthesis  methodology  in  both  low- 
and  high  authority  control  techniques.  Section  5  provides  summaries  of  the 
analytical  designs  for  stability  augmentation  systems  for  the  DARPA  Strawmen 
configurations  studied  under  Phase  I.  Finally,  experimental  results  developed 
under  IRAD  which  are  relevant  to  Phase  I  objectives  and ,  in  particular ,  to 
follow-on  efforts  are  included  in  Section  6. 


Section  2 
APPROACH 


2.1  INTRODUCTION  AND  METHODOLOGY 

Stability  augmentation  or  the  addition  of  closed-loop  damping  to  large  flexible 
space  structures  has  been  the  subject  of  numerous  recent  papers  (Refs  [1-5]). 
These  publications  have  focused ,  for  the  most  part ,  on  single  methodologies  or 
techniques  for  synthesizing  a  control  law.  This  report  is  a  synopsis  of  research 
work  performed  over  the  past  year  and  endeavors  to  match  various  synthesis 
techniques  and  controller  forms  to  specific  system  performance  requirements.  In 
so  doing,  it  evolves  a  general  approach  to  structural  stability  augmentation.  More¬ 
over,  this  combination  of  techniques  permits  guarantees  of  system  robustness, 
if  not  performance,  which  would  otherwise  be  unattainable.  The  methodology 
advocated  here  requires  several  complex  techniques  indicated  below. 

The  general  requirements  for  stability  augmentation  are  listed  in  Table  1.  It 
should  be  noted  that  shape  control  is  a  different  process  from  removal  of  vibra¬ 
tional  energy  from  the  structure  as  noted  in  the  table.  The  first  two  categories 
require  large  increases  in  structural  damping  at  low  frequency  either  because  the 
attitude  control  system  is  interacting  with  the  structure  or  because  maneuver  in¬ 
duced  structural  ringing  must  be  removed  to  meet  mission  requirements.  Modify¬ 
ing  structural  admittance  functions  at  specific  low  frequencies  generally  re¬ 
quires  substantial  knowledge  of  mode  shapes  and  frequency  quantities  which 
are  known  less  accurately  at  higher  frequency.  The  third  category  arises  from 
the  need  to  control  high  frequency  settling  times  and  to  suppress  vibration  prop¬ 
agation  from  sources  imbedded  in  the  structure  (e.g.  ,  coolers,  CMG's,  etc.). 
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2.2  BASIC  CONTROLLER  FORMS 


Two  basic  classes  of  controls  emerge  from  these  requirements.  Low- Authority 
Structural  Controls  (LASC)  which  provide  limited  (broadband)  modal  damping 
may  be  synthesized  using  perturbation  techniques  as  shown  by  Aubrun  [Ref  2 
and  Section  4.1].  For  systems  restricted  to  collocated  rate  sensing  and  comple¬ 
mentary  actuation  (decentralized,  no  cross- couplin g)  ,  Aubrun  shows  that  stability 
is  guaranteed  over  the  system  bandwidth  and,  with  suitable  precautions  as  dis¬ 
cussed  below,  the  system  is  always  stable.  High- Authority  Structural  Control 
(HASC)  which  provides  arbitrary  damping  ratios  in  low  frequency  modes  by  ex¬ 
ploiting  knowledge  of  system  dynamics  is  generally  synthesized  as  a  multi-input/ 
multi-output  (centralized)  control.  A  direct  application  of  standard  control  design 
procedures  leads  to  serious  stability  problems.  Within  the  HASC  bandwidth,  spill¬ 
over  in  modeled  and  unmodeled  modes  must  be  bounded  by  reduction  of  parameter 
sensitivity  or,  ultimately,  by  system  identification  prior  to  imposing  the  control. 
Beyond  the  HASC  bandwidth  uncontrolled  phase  shifts  introduced  by  lightly  damped 
modes  can  cause  system  instability.  Synthesis  of  the  HASC  system  which  forces 
rapid  control  effort  roll-off  with  increasing  frequency  coupled  with  an  overlapped 
LASC  system  to  lower  modal  q  in  the  crossover  region  significantly  increases  sys¬ 
tem  robustness.  This  relationship  is  illustrated  in  Fig.  1.  The  methods  for  guaran¬ 
teeing  stability  are  summarized  in  Table  2.  Fundamentally,  HASC  stability  due  to 
interaction  with  high  frequency  modes  is  assured  by  limited  high  frequency  HASC 
authority  and  a  decentralized  LASC  system  for  which  critical  modes  at  HASC  cross¬ 
over  are  controllable  (fixed  by  LASC  actuator  placement). 

The  selection  of  a  control  law  using  a  centralized  form  (LASC  or  HASC)  may  proceed 
either  by  using  perturbation  techniques  based  on  a  formula  by  Jacobi  or  by  mini¬ 
mization  of  a  quadratic  cost  functional  with  state  and  control  weighting  which  are 
frequency  dependent.  This  frequency  dependence  forces  a  dynamic  controller 
form  with  an  augmented  state  vector  dictated  by  the  frequency  roll-off  selected 
by  the  designer. 


Fig.  1  Integrated  LASC/HASC  Stability  Augmentation 


2.3  FEATURES  OF  THE  CONTROL  SYNTHESIS  METHODOLOGY 


The  following  features  of  the  control  design  methodology  provide  a  systematic 
approach  to  meet  system  performance  goals. 

(a)  The  basic  method  consists  of  three  steps,  (i)  model  reduction, 

(ii)  filter  /control  law  specification,  and  (iii)  control  law  simplifi¬ 
cation.  This  provides  a  good  model  in  the  control  design  stage. 

The  resulting  controller  may  be  further  simplified  in  the  last  step. 

(b)  Specific  attention  is  given  to  poles ,  zeros  and  residues  in  the  model 
reduction  stage.  This  resulting  reduced  order  model  then  provides 
a  good  approximation  of  system  behavior  in  the  frequency  range  of 
interest . 

(c)  For  collocated  systems  (sensors  and  actuators  at  the  same  location) , 
a  general  control  design  approach  is  proposed  based  both  on 
modern  control  theory  and  root  perturbation  techniques. 

(d)  A  unique  linear  quadratic  -  Gaussian  design  method  is  developed  in 
which  the  performance  index  is  shaped  with  frequency.  This  results 

in  control  laws  which  have  minimum  influence  on  high  frequency  modes. 
This  provides  a  mechanism  to  minimize  the  possibility  of  high  frequency 
modes  going  unstable. 

(e)  The  controller  is  treated  as  a  dynamic  model  with  specified  inputs 
and  outputs.  The  model  reduction  methods  are  then  used  to  simplify 
the  controller. 

(f)  Techniques  for  optimal  actuator /sensor  locations  are  studied. 
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Section  3 

STRUCTURAL  MODELING  TECHNIQUE 


3. 1  FINITE-ELEMENT  MODELS 

A  central  issue  in  the  active  control  of  space  structures  is  the  development  of 
"correct”  mathematical  models  for  the  open  and  closed  loop  dynamical  plants. 

While  a  great  deal  of  polemic  has  been  generated  in  recent  times  against  finite- 
element  structural  programs  and  their  assorted  pitfalls,  programs  such  as 
NASTRAN,  SPAR,  ASTRO,  etc.,  are  never-the-less  the  primary  tools  for 
generating  dynamical  models  of  conceptual  spacecraft  whose  structure  cannot 
be  idealized  by  simple  models  of  beams,  plates,  beams  with  lumped  masses,  and 
so  on. 

Finite-element  structural  programs  were  used  extensively  during  Phase  I  of  the 
ACOSS  Program  in  order  to  model  the  DARPA  Strawmen  dynamical  plants 
(HALO,  MM  WAVE,  ADOPT’s  5/4  and  12,  CSDL  example)  to  which  various  con¬ 
trollers  were  applied.  While  this  modeling  was  not  heavily  emphasized  by  the 
contract  requirements,  it  turned  out  to  be  an  important  activity  in  the  sense 
that  insights  were  gained  into  the  various  aspects  of  the  process  and  the  emer¬ 
gence  of  a  general  methodology  was  identified.  The  essential  observation  was 
made  that  a  single  finite-element  model  of  a  structure,  postulated  in  an  a  priori 
fashion,  was  not  necessarily  the  best  modal  representation  suitable  for  the  closed- 
loop  model.  After  insight  was  gained  into  the  control  policy  to  be  applied,  and, 
in  particular,  into  the  type  of  actuctors  (and  their  locations)  required  to  im¬ 
plement  it,  it  became  very  plausible  that  a  more  refined  open-loop  model  would 
improve  modal  convergence  of  the  closed-loop  model  and  hence  alleviate  the 
difficulties  of  modal  truncation.  Implicit  in  this  approach  is  the  requirement  that 
a  structure  can  be  examined  geometrically  by  computer- graphic  animation  tech¬ 
niques  which  clarify  the  physical  interaction  between  actuators  and  the  local 
structural  deformations  which  they  impose.  In  this  sense,  a  structure-to-be- 
controlled  becomes  much  more  than  a  set  of  constant  coefficient  matrices  given 
once  and  for  all. 


The  structural  design  iteration  procedure  identified  at  the  end  of  Phase  I  of  the 
ACOSS  program,  could  typically  involve  the  following  steps: 

1.  Design  the  global  structure,  and  '’tune”  local  substructures.  In  the 
HALO  design ,  for  example ,  the  tuning  of  the  flexible  mirror  sub¬ 
structures  was  to  a  large  extent  independent  of  the  overall  design 
which  would  treat  these  mirrors,  in  a  first  approximation,  as  rigid 
plates  interconnected  by  elastic  beams. 

2.  Examine  the  elastic  behavior  (generally  via  the  undamped,  free /free 
eigenmodes)  using  computer- graphic  animation  such  as  available  on 
the  ADAGE  computer  in  LMSC's  STARLAB.  What  is  required  here  is 
the  observation  of  the  physical  motions  occurring  locally  in  the  neigh¬ 
borhood  of  nodal  and  antinodal  points  of  the  dominant  structural  modes. 

3.  Infer  from  these  geometrical  motions  (translations,  rotations,  and  com¬ 
binations  of  both)  preliminary  actuator  locations  and  actuator  types  to 
be  used.  These  could  consist  of  forces,  force  couples,  torques,  torque 
couples,  etc.  Depending  upon  the  structural  design  and  the  mission 
requirements,  trade-offs  can  be  made  between  inertial  and  intra- 
structural  actuators.  This  step  is  particularly  important  for  decentral¬ 
ization  colocated  LAC  sensor /actuator  pairs,  and  serves  indirectly  to 
corroborate  the  information  derived  from  Aubrun’s  controller  effective¬ 
ness  coefficients* 


or,  equivalently,  from  the  Skelton-Hughes  controllability  maps  when 
these  can  be  constructed. 

4.  Synthesize  preliminary  controls  and  test  performance  via  closed-loop 
eigenanalysis .  Performance  criteria  may  be  derived  from  the  insensitivity 
of  the  zeroes  of  closed-loop  transfer  functions  to  modal  truncation  (see 
next  section),  from  convergence  of  phase  and/or  gain  measures  vs 
number  of  modes,  or  from  other  measures.  If  performance  is  not  satis¬ 
factory,  proceed  to  next  step. 

5.  Modify  the  boundary  conditons  of  the  structure  by  introducing,  at 
actuator  locations,  "fictitious"  inertias  (translational  or  rotational) 
whose  D'Alembert  reactions  (forces  or  torques)  will  mimic  the  effect  of 
the  corresponding  actuators.  In  structural  dynamics  this  procedure 
is  equivalent  to  the  use  of  certain  classes  of  quasi-static  modes.  For 
example,  a  beam  which  is  torque-controlled  at  one  end  is  better  repre¬ 
sented  by  modes  which  correspond  to  the  presence  of  a  suitable  rotary 
inertia  at  the  controlled  end  rather  than  by  free-free  modes,  even  if 

a  lumped  mass  is  added  at  the  controlled  end  to  model  the  actuator. 


♦See  eg  (8)  in  Section  4.1,  Low  Authority  Control  Synthesis: 


6.  Repeat  eigenanalysis  of  the  modified  structure  to  generate  faster  con¬ 
verging  modes  vis-a-vis  performance  criteria. 

The  above  design  iteration  steps  represent,  in  general,  extensive  work,  but 
experience  gained  in  Phase  I  of  the  ACOSS  program  clearly  indicates  that  the 
integration  between  controls  and  structural  dynamics  expertise  is  far  from  com¬ 
plete.  The  effectiveness  of  procedures  such  as  outlined  above  can  only  be  eval¬ 
uated  in  detailed  applications. 
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3.2  DISTRIBUTED  MODELS 


Such  models  may  be  obtained  for  structures  that  can  be  conveniently  described 
by  partial  differential  equations  (PDEs).  Although  it  seems  unlikely  that  a  com¬ 
plex  space  structure  will  ever  fall  into  that  category,  interest  in  this  type  of 
modeling  is  threefold:  1)  for  laboratory  experiments,  where  the  test  structures 
may  be  simple  enough,  an  analytical  model  would  be  the  most  accurate  and  there¬ 
fore  of  greater  value  for  evaluating  the  experimental  results;  2)  PDEs  may  be 
used  directly  for  control  synthesis;  and  3)  PDEs  offer  a  direct  way  for  evaluat¬ 
ing  the  finite -element  models  and  understand  their  deficiencies. 

Distributed  models  may  be  used  in  two  different  ways,  both  starting  with  partial 
differential  equations  which  are  assumed  to  be  a  complete  description  of  the  dy¬ 
namical  characteristics  of  the  system. 

The  first  method  converts  these  PDEs  and  associated  boundary  conditions  into  an 
infinite  set  of  ordinary  differential  equations.  A  small  subset  of  these  ODEs  is  then 
chosen  as  a  model,  and  this  is  strictly  equivalent  to  the  finite-element  approach, 
except  that  a  finite  set  of  modes  is  obtained  analytically  and  not  via  a  numerical 
eigensystem  decomposition.  This  method  was  used  in  all  the  beams  and 
experiments  described  in  Section  6  and  were  found  to  be  very  accurate  and 
useful. 

The  second  method  starts  with  the  same  PDEs  but  does  not  expend  the  solution 
into  modes.  Better,  it  directly  obtains  an  exact  analytical  transfer  function  which 
can  be  used  in  classical  control  theory,  or  in  the  more  elaborate  -  distributed  op¬ 
timal  control  theory  if  desired.* 

Since  approximate  transfer  functions  can  also  be  obtained  from  model  expansion 
of  the  first  method,  comparisons  can  be  made  between  finite-element  representa¬ 
tions  (identical  to  the  first  method)  and  the  exact  model  (second  method).  The 


*J.  A.  Breakwell,  ’’Automatic  Control  of  Flexible  Spacecraft,"  Ph.D.  Thesis, 
Stanford  University,  to  appear  Dec  1979. 


particular  problem  of  a  free-free  beam  was  treated  using  these  two  approaches 
and  analytical  expressions  obtained  for  the  continuous  and  discrete  (modal) 
transfer  functions  are  shown  in  Fig.  1.  Two  types  of  transfer  function  appear 
in  Fig.  1,  i.e.,  displacement /force  and  rotation /torque  at  one  end.  All  symbols 
are  defined  in  Fig.  1.  The  discrete  transfer  function  T  (/3;n)  is  a  function  of 
the  number  n  of  modes  retained  in  the  series  expansion. 

Figure  2  shows  a  graphical  comparison  of  the  exact  rotation-to-torque  transfer 
function  Tc(/3)  (obtained  from  the  continuous  model)  and  the  truncated  (discrete) 
transfer  functions  T^(/3;n)  when  7  terms  or  20  terms  are  retained,  i.e. ,  the 
beam  is  modeled  by  one  rigid  body  mode  and  either  6  or  19  structural  modes. 

The  graph  clearly  shows  the  error  in  the  zero-crossings ,  while  the  other  parts 
of  the  curves  are  in  excellent  agreement.  Because  the  zeroes  of  the  open-loop 
system  play  an  important  part  in  the  closed-loop  system,  their  insensitivity  to 
modal  truncation  is  of  prime  importance  to  the  control  design.  Zero  crossing 
errors  can  be  reduced  by  increasing  the  number  n  of  modeled  modes,  but  a 
satisfactory  result  will  require  a  very  large  value  for  n.  These  curves  indicate 
that  a  much  better  approximation  could  be  obtained  by  adding  a  constant  bias 
term  (negative  in  this  case)  to  the  discrete  transfer  function. 

A  more  systematic  view  of  this  problem  can  be  seen  in  Fig.  3  which  shows  the 
partitioning  of  the  infinite  series  representing  the  exact  transfer  function  for 
the  case  of  a  general  undamped  structure.  The  0.'s  are  the  mode  shapes  at  the 
location  where* the  displacement-to-force  ratio  is  measured  (or  rotation-to-torque, 
etc.).  The  number  of  modes  to  be  controlled  is  nc,  the  number  of  modes  re¬ 
tained  in  the  model  of  the  structure  is  n  . 

r 

In  the  straight  truncation  case  one  has  np  =  nQ ,  and  a  substantial  error  is 
committed  in  representing  the  zeroes  of  the  system,  leading  to  erroneous  control 
design.  The  expanded  truncation  uses  a  larger  model,  but  still  a  significant 
error  may  remain.  The  "rounded"  transfer  function  approximates  the  remainder 
of  the  infinite  series  by  a  constant  bias  term.  (For  these  terms,  s  is  indeed 
much  less  than  the  w's,  hence  the  approximation.) 
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BEAM  DENSITY  S:  LAPLACE  TRANSFORM  VARIABLE  y:  BEAM  DEFLECTION  AT  ONE  END 


Transfer  Functions  for  Beam  (II) 


T6  e>e  e.erMwiED  Ktor  betiumed 


Fig.  3  Modal  Partitioning  and  Approximate  Transfer  Functions 


Because  all  these  approximations  depend  on  the  convergence  of  the  series,  there 
are  open  questions  regarding  the  use  of  other  kinds  of  mode  shape  functions 
which  may  improve  convergence.  For  instance  mode  shape  derivatives  (used  to 
express  rotation  and  torque  properties)  converge  much  slower  than  the  usual 
displacement  mode  when  free-free  modes  are  used.  However,  using  mode  shapes 
which  reflect  better  the  boundary  conditions  imposed  by  a  torquer  may  drastically 
improve  convergence. 

Figure  4  shows  an  example  of  a  more  general  6th  order  system  transfer  function 
and  three  different  ways  to  obtain  reduced  order  transfer  functions.  Numerical 
values  are  displayed  in  Fig.  5.  The  best  results  are  obtained  with  the  rounded- 
off  T.F.  where  not  only  the  zeroes  are  well  represented  but  also  the  residues  of 
the  retained  poles  (which  is  not  the  case  for  the  zero-retaining  T.F.). 


EXAMPLE 


Fig.  4  Transfer  Function  Simplifications  (Truncation  vs  Round-Off) 


Transfer  Function  Simplifications  (Truncation  vs  Round-Off) 


Section  4 

CONTROL  SYNTHESIS  TECHNIQUES 

4.1  LOW  AUTHORITY  CONTROL  SYNTHESIS 

The  following  controller  design  synthesis  is  based  on  J.  N.  Aubrun's  Low 
Authority  Control  (LAC)  Theory  (Ref  [*]).  To  describe  the  methodology,  only 
an  outline  of  the  required  theory  will  be  given;  for  proofs  and  additional  details, 
see  Ref  [  *] . 

Consider  the  system  described  in  state-space  form  by: 


=>X  =  (F  +  GCH)  X  (1) 

(Closed-Loop  Dynamics) 

For  sufficiently  "small"  controls  C  (see  Ref  [*]),  let 

GCH  2  dF  (2) 

be  considered  as  a  perturbation  of  F ,  so  that  the  closed-loop  dynamics  may  be 
written  as 

X  =  (F  +  dF)  X  (3) 

Let  Xn  denote  the  n-th  (complex)  root  of  F,  and  denote  by  Ln  ,  Rn  respectively 
the  corresponding  left  and  right  eigenvectors  of  F,  i.e. 

[*]  J.  N.  Aubrun,  "Theory  Of  The  Control  Of  Structures  By  Low  Authority 

Controllers",  AIAA  Conference  On  Large  Space  Platforms:  Future  Needs  And 
Capabilities,  Los  Angeles,  Calif/Sept  27-29,  1978,  Page  No,  78-1689;  also, 
AIAA  J. ,  Guidance  and  Control,  revised  version  (in  press),  1979. 


(4) 


F  L  =  \  L 

n  n  n 


F  R  =  X  R 
n  n  n 


T 

normalized  so  that  L  R  =1. 

n  n 

For  "small"  C,  let  d\  denote  the  eigenvalue  shift  corresponding  to  dF. 
Then,  as  shown  in  egs  (19)  and  (27)  of  Ref  [*] ,  Jacobi's  formula  leads  to  the 
fundamental  LAC  root  shift  formula 


d*n  a  dFRn 

*  GCHR„  (5) 

*  K)cK) 

AT  R 

where  <t>  =  L  G  and  <t>  =  HR _  are  generalized  actuator  and  sensor 

n  n  n  n  ■** - 

modes  respectively.  The  fundamental  LAC  formula  expresses  then  the  eigenvalue 
shifts  by  the  bilinear  form 

dA  s  C  <t>A  0R  (6) 

n  ar  an  rn 

a,r 

where  a  =  1,  2,  ....  N  and  r  =.  1,  2,  N  are  indices  describing,  re- 
spectively,  actuator  and  sensor  locations  (nodes). 

Consider  next  the  case  where: 

.  *) 

(i)  sensors  and  actuators  are  physically  colocated 

(ii)  sensors  and  actuators  are  of  "corresponding  types"  (i.e.  translation/ 
force ,  rotation  /torque ) 

(iii)  sensors  measure  rates 


*)  This  does  not  necessarily  require  that  the  sensor- to-actuator  feedback  is  also 
colocated . 
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then  it  can  be  shown  (Ref  [*])  that  H  =  G1 .  If,  furthermore,  there  exists  a 
right  and  left  eigenvector  for  root  n,  it  can  also  be  shown  that 


and  thus 


G  L  =  HR 
n  n 


GTLn  =  HRn  =  **  <»<•„> 


which  is  an  equivalent  re-statement  of  the  classical  reciprocity  relations  between 
actuators  and  sensors  for  colocated  systems  of  the  type  considered  above. 

From  (8),  the  fundamental  LAC  formula  (6)  becomes  then  the  following  quadratic 
form  (for  each  n) : 


d\>  a  £c 


0  0 

ar  an  yrn 


The  above  formula  is  the  basic  LAC  prediction  formula  for  the  root  shifts  (hence 

closed-loop  damping)  produced  by  sufficiently  "small'’  control  gains  C  (see 

Ref  [*]).  To  synthesize  the  gains  C  ,  let  (dx  )p  denote  the  predicted  root 

ar  n  jf 

shifts  given  by  the  above  formula,  and  let  (dxn)j)  (given  numbers)  denote  the 
desired  root  shifts  imposed  by  the  LAC  controller  design.  Then  the  gains  C 

ar 

are  chosen  so  as  to  minimize  the  weighted  quadratic  cost  functional 


J(C> 


(dAn)P 


(d  An)D 


in  which  the  modal  weights  WR  help  specify  pole  locations,  and  the  term^  C&  2 

a,r 

improves  robustness  of  the  controller.  Since  the  cost  functional  J(C)  is  quadratic 
in  C ,  the  gains  can  be  obtained  algebraically  by  solving  the  linear  equations 


it  A  T 

Jr  -  °  f°r  c  • 
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In  order  to  carry  out  this  procedure,  we  need  to  (i)  relabel  the  gain  matrix  C 

-  ar 

(a  =  1,  ....  Na,  r  =  1,  ....  Np)  as  a  gain  vector  Ca  (a  =  11,  12,  ...lNr, 

21,  22, . -N  N  ) ,  (ii)  relabel  the  multiply-indexed  modal  coefficients 

a  1  ( o) 

0an  0rn  as  a  niodal  coefficient  matrix  ^nQ;  ,  using  the  same  correspondence 

between  the  single  index  a  and  the  pair  of  indices  a,  r  as  used  in  (i),  and 

finally  (ii)  define  a  diagonal  weighting  matrix  W  whose  elements  are  W  .  For 

n 

simplicity,  denote  the  desired  root-shifts  (d Xn)D  by  dR  ,  considered  as  com¬ 
ponents  of  a  vector  d. 

With  the  above  relabeling,  the  weighted  quadratic  cost  functional  (10)  can  be 
rewritten  as 


J(C)  = 


2>n£ 


n 


$ 


(2) 


na 


Ca- 


n 


£Ca 


"a 


a 


or  in  matrix  form 


(ID 


J(C)  =  (4>( 2)C  -  d)T  W  (*(2)C  -  d)  +  CT  C 


(12) 


From  (11)  we  can  calculate  the  partial  derivatives  (where  p  is  an  index  having 
the  same  range  as  a ) : 


2  aCy3  —  n  |  n  n 


a 


*n(S  +  S 


(13) 


or  equivalently  in  matrix  form  (I:  identity  matrix) 


1  aj 

2  Tc 


($(2))T  w$(2)  + 1 


C  -  ($(2))T  Wd 


(14) 


aj 


Thus  the  condition  =  0  leads  to  the  LAC  gain  synthesis  formula 
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The  above  LAC  gain  systhesis  procedure  was  carried  out  for  the  CSDL 
tetrahedral  truss  example,  using  rate-feedback  only  and  an  upper  triangular 
form  for  the  gain  matrix  C.  The  CSDL  example  will  be  discussed  in  Section  5.1 
below . 


4.2  HIGH-AUTHORITY  CONTROL  METHODS 

4.2.1  Control  Problem 

Introduction 


The  application  of  advanced  control  design  techniques  requires  the  use  of 
models.  Dynamics  of  continuous  structures  is,  in  general,  described  by  a  set 
of  partial  differential  equations  of  the  wave  type.  These  equations  may  be 
derived  for  simple  structures  like  beams  and  plates  but  are  difficult  for 
practical  structures  used  in  space.  In  addition,  even  though  control  laws 
can  be  designed  for  such  structures,  in  concept,  the  solutions  are  too  diffi¬ 
cult  for  implementation. 

A  more  practical  approach  for  modeling  the  structures  uses  the  modal 
approach.  In  this  approach  the  basic  functions  describing  the  dynamic 
system  behavior  are  written  in  terms  of  spatial  expansions  and  the  dynamics 
of  the  basis  function  is  described  by  a  coupled  ordinary  differential  equation. 
This  differential  equation  forms  the  basis  of  control  design  approaches. 

4. 2. 1.1  Modal  Representation.  Finite  order  modal  representations  are  obtained 
either  from  simplification  of  partial  differential  equations  representing  the 
dynamics  of  the  structures  or  a  direct  derivation  based  on  finite  element  analysis. 
In  the  finite  element  analysis,  the  structure  is  divided  into  a  set  of  lumped 
masses  and  the  differential  equations  are  derived  with  appropriate  boundary 
conditions  at  each  lumped  mass.  The  details  of  the  finite  element  approach  may 
be  found  in  other  references. 

A  finite  order  lumped  mass  model  for  the  structure  is  written  as 

M £  +  Ct  +  Kf  =  Gu  (1) 
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where  f  are  translational  or  rotational  deflections  and  u  is  the  input.  M,  C, 

K  and  G  are  mass,  damping,  spring  and  input  distribution  matrices.  This 
model  is  typically  of  high  order,  depending  upon  the  number  of  elements  used 
in  finite  element  analysis  or  the  grid  size  used  to  convert  the  partial  differential 
equation  into  an  ordinary  differential  equation. 


The  model  may  be  converted  into  a  state  space  formulation  by  defining  the  state 
vector  as  follows: 


1/2 

The  states  may  alternately  be  defined  as  follows,  where  M  is  a  symmetric 
square  root  of  M. 


X  = 


(4) 


The  corresponding  state  equations  are 


With  this  selection  of  the  state,  the  damping  and  the  spring  matrices  are  sym¬ 
metric.  These  equations  are  written  as: 


_d_ 

dt 


(6) 


The  following  transformation  could  convert  A  into  a  diagonal  matrix. 


(7) 
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where  T  is  any  matrix  which  makes  TAT  diagonal.  The  state  equations 
then  become 


_d 

dt 


where  Aq  is  a  diagonal  matrix. 


(8) 


Inertial  devices  measure  absolute  values  of  translational  or  rotational  positions, 
velocities  or  accelerations  at  several  points  on  the  structure.  Optical  devices 
or  other  systems  may  measure  relative  position,  velocity  or  acceleration  at 
two  or  more  points  on  the  structure.  The  measurements  are  written  as 


position : 

y  =  Cji 

II 

1  1 - 1 

O 

i-* 

V 

t, 

r*  -i 

velocity : 

y  -  c2i 

II 

'  o 

o 

to 

L 

' — i — 

vTrr|-|<rrr| 

l  1 

acceleration : 

y  =  c  3t 

=  -C„  A 

3  o 

:  -C 

3  D°]  [  i 

[-]  +C3B„U 

In  general,  therefore,  the  measurements  are  expressed  as: 


y 


(9) 


4.2. 1.2  Control  Problems  and  Requirements.  The  basic  problem  with  many 
structural  dynamics  problems  is  the  low  damping  ratios  of  system  modes.  The 
systems  of  interest  to  us  must  perform  several  tasks  which  require  both  model 
following  and  stability  augmentation. 


The  structure  models  are  often  known  accurately  at  low  frequencies.  In  addition, 
the  low  frequency  modes  have  large  settling  times.  Therefore,  control  activity 
is  more  efficiently  utilized  in  controlling  low  frequency  modes.  A  model  reduction 


approach  is  required  to  reduce  the  partial  differential  equation  model  or  high 
order  lumped  parameter  model  into  a  low  order  model.  These  techniques  must 
not  only  represent  low  frequency  behavior  faithfully  but  must  also  include  the 
impact  of  high  frequency  behavior  in  low  frequency  models. 


The  control  design  problem  is  complicated  because  the  model  being  lumped  mass 
and  of  low  order,  approximates  the  dynamic  behavior  of  the  structure.  The 
control  design  procedure  must  appreciate  this  limitation  of  the  model.  A  direct 
application  of  the  linear-quadratic  Gaussian  methods  produces  the  control  and 
observation  spillover  problem,  where  the  control  law  adversely  affects  one 
or  more  unmodeled  modes.  The  control  structure  must  avoid  or  minimize  the 
spillover  problem. 


The  general  control  law  resulting  from  the  above  developments  is  likely  to  be  of 
high  order.  Techniques  are  developed  to  reduce  controller  order  without 
compromizing  closed  loop  performance,  if  possible. 

Figure  1  shows  the  overall  control  design  procedure.  The  overall  procedure 
consists  of  three  major  steps: 

(1)  Model  Reduction:  Development  of  control  design  models  from 
distributed  parameter  or  finite  element  descriptions. 

(2)  Control  Design:  Design  of  a  controller /filter  or  output  feedback 
logic  to  specify  the  control  structure. 

(3)  Controller  Simplification :  Reduction  of  the  controller  obtained  in 


the  previous  step  to  simplify  implementation  and  for  rebustification. 

The  steps  in  the  overall  procedures  must  be  carefully  integrated  to  support 
assumptions  and  approximations  made  in  each  of  the  steps.  This  integration 
and  the  development  of  specific  design  methodologies  which  support  this 
integration  is  the  basis  of  our  approach  for  design  of  control  system  for  space 
structures. 


r 
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Fig.  1  Control  Synthesis  Approaches 


4.2.2  Model  Reduction 


4. 2. 2.1  Requirements.  Partial  differential  operations  describe  the  dynamics  of 
large  space  structures.  Lumped  mass  approximation  of  this  dynamics  via  finite 
element  analysis  or  direct  discretization  leads  to  high  order  ordinary  differential 
equation  models.  Though  this  model  form  is  suitable  for  the  application  of  many 
well  developed  control  design  methods,  the  high  order  creates  a  problem  both  of 
computation  and  of  implementation.  In  addition,  higher  order  modes  in  continuous 
structures  are  more  sensitive  to  microscopic  changes  in  the  physical  structure 
than  low  order  modes.  The  description  of  high  order  modes,  therefore,  has  large 
errors  resulting  from  the  finitie  element  analysis  and  approximations  in  modeling 
the  physical  structure.  Lumping  of  two  truss  elements  may  have  a  minor  effect 

on  the  first  mode  but  may  completely  alter  the  mode  shape  of  a  high  frequency 
system  mode. 

Model  reduction  is  required  for  robustness  and  to  simplify  control  law  computation 
and  implementation. 

Criteria  and  computation  procedures  for  model  reduction  are  discussed  in  this 
chapter.  The  reduced  model  must  contain  essential  elements  of  the  dynamics 
such  that  the  resulting  control  law  is  suitable  for  the  continuous  structure  des¬ 
cribed  by  partial  differential  equations.  Roles  of  poles,  zeros  and  residues  of 
transfer  functions  in  closed  loop  control  design  are  discussed  first  followed  by 
selection  of  specific  criteria  and  numerical  procedures. 

4. 2. 2. 2  Poles,  Zeros  and  Residues.  In  linear  systems,  the  closed  loop  behavior 
is  dictated  significantly  by  several  transfer  function  properties.  Three  of  the 
most  significant  properties  are  the  poles,  the  zeros,  and  the  residues. 
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If 


The  roles  of  transfer  function  properties  are  best  explained  by  considering  a 
single  input  single  output  (SISO)  system  with  output  y  and  input  u. 

y(s)  =  T(s)  u(s)  (1) 

The  characteristic  values  Aj,  i  =  1,  2  •••  n  of  the  T(s)  denominator  are  the  open 
loop  system  eigenvalues  and  indicate  system  stability  properties.  The  transfer 
function  may  be  written  in  terms  of  zeros  Z.f  i  =  l,  2  •  •  .  n  or  residues  r., 
i  =  1,  2  •  •  •  n  as  follows 


T(s) 


K(s  -  Zp  (s  -  Z2) 
(s  -Xj)  (s  -  X2) 


<s  -  Zm) 


(s-x„) 


:  Zeros  Representation 


r. 

l 


"*i 


:  Residue  Representation 
(2) 


Consider  now  a  feedback  matrix  C(s)  with  gain  a.  The  closed  loop  transfer 
function  is 


Tc(s) 


T(s) 

1  +  «C(s)  T (s) 


(3) 


which  may  be  written  as 


K(s  -  z-)  •  •  •  (s  -  z— ) 

T  (si  =  - ■=-= - = - !_ m 

c1  (s  -  \1>  (s  -  \2)  •••  (s  -  X-) 


(4) 


It  is  easy  to  see  that  for  smaller  ot. 


=  aC(A.)  r. 


(5) 


Thus  r.  dictates  the  behavior  of  the  pole  for  small  gain.  When  a  is  large,  the  finite 
closed  loop  poles  are  the  zeros  of  C(s)  and  T(s). 
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To  summarize,  the  residues  of  the  transfer  function  describe  the  low  gain  proper¬ 
ties  and  the  zeros  the  high  gain  properties.  Both  zeros  and  residues  are  important 
in  closed  loop  control  design. 

4. 2. 2. 3  Numerical  Procedures.  Consider  a  Nth  order  linear  system  with  x, 
input  u  and  output  y 


x  =  Fx  +  Gu  (6) 

y  =  Hx  +  Du  (7) 


The  input-output  transfer  function  is  given  by 

y  =  (H(sl  -  F)_1  G  +  D)u 


(8) 


or 


y(s)  1_ 

u(s)  A(s) 


N 

£  ^(i)s1 

i=0 


(9) 


The  eigenvalues  of  A(s)  define  system  poles,  while  B(i)  are  matrices  of  numerator 
polynomial.  The  zeros  of  y^  and  Uj  transfer  function  is 


yk(s) ■  i  f 

Uj(s)  A(S)  Ajj 


Bkl(i)si 


Kkl  <s 


Zlkl)  (s 


Z2kl>  • 


(s  -  Z 


Nkl) 


-  (s  -  Aj)  (s  -  *.) 


(s  ^N^ 


N 


a  E 

i=l 


rikl 

(s  -  X.) 


(10) 


Different  approaches  are  useful  for  maintaining  zeros  or  residues  in  reduced 
order  models. 
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Maintaining  Residues  and  Poles 

A  modal  decomposition  provides  the  framework  for  reducing  arbitrary  linear  models 
to  design  models  with  same  residues  and  poles  as  the  high  order  models.  The  pro¬ 
cedures  are  well  known  [  2]  and  their  application  to  model  reduction  is  described 
below . 

The  linear  Eqs.  (6)  and  (7)  can  be  transformed  to  block  diagonal 
the  n  x  n  dynamics  matrix,  F,  has  no  repeated  eigenvalues: 

x  =  Tz  S 
z  =  Az  +  u 
y  =  HTz  +  Du 

where  A  is  an  n  x  n  block  diagonal  matrix,  T  is  an  n  x  n  matrix  composed  of 
the  column  eigenvectors  of  F ,  Z  is  an  n  x  1  modal  coordinate  vector ,  and  E  is 
the  n  x  m  modal  control  distribution  matrix.  Also, 

FT  =  TA  (14) 


form  assuming 

(11) 

(12) 

(13) 


The  system  of  Eqs.  (11)  through  (13)  can  be  partitioned  into  a  set  of  q  states 
and  q  eigenvalues  (time  constants)  and  n-q  states  and  eigenvalues  as  follows: 


f  xll  _  [*T11  !  T12lKl 

LkV]  "  1*22!  vlhJ 


(15) 


(16) 


where  x^  and  z^  are  n  x  1  vectors  partitioning  the  states  and  modes  and  Xg  and 
Zg  are  (N-n)  x  1  vectors  partitioning  the  remaining  states  and  modes. 

In  the  representation  of  the  above  equation  it  is  easy  to  verify  that  the  residues 
of  the  poles  are  appropriate  products  of  terms  in  matrices  HT  and  H.  The  residue 
maintaining  model  reduction  procedure  is  simply  derived  as  in  the  following. 
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(17) 


N-n  poles  are  eliminated  by  setting  Z2  =  0.  This  gives 


Z2  ~  ~  A2*  *  u 


The  reduced  order  model  is 


Z1  ~  A1Z1  +  *iu 


(18) 


Y  =  HT 


-A'1  -u 


+  Du 


=  H’Z1  +  D'u 


(19) 


Many  model  reduction  procedures  do  not  modify  the  coefficient  of  u  in  the  measure¬ 
ment  equation.  It  is  clear  from  the  above  equations  and  subsequent  examples  that 
such  a  change  is  essential. 


Maintaining  Zeros  and  Poles 

If  the  transfer  function  between  output  y^  and  input  u^  is  simplified  to  keep  the 
first  w  poles  and  m  zeros ,  the  following  form  is  obtained . 


If  the  first  n  poles  and,  m  zeros  the  simplified  transfer  function  is 


yk(s) 

"Ks) 


N 

n 

i=m+l 


(-Zikl> 


TT-  (~V 

i=n+l  1 


(s  "  Zlkl)  ***  (s  '  Zmkl) 


(s  -  \.) 


(s  -  Xn) 


j 

(20) 


Note  that  the  gain  has  to  be  modified  because  high  frequency  poles  and  zeros 
have  steady  state  effects. 


The  reduced  model  is  obtained  by  performing  a  similar  reduction  on  each  transfer 
function  and  then  detaining  a  state  space  realization  from  the  complete  transfer 
function.  The  process  of  obtaining  the  minimum  state  space  realization  is  not 
simple  (except  for  single  input  multi  output  or  multi  input  and  single  output  sys¬ 
tem).  In  fact  in  multi  output  systems,  this  procedure  may  even  lead  to  an  increase 
in  model  order. 

Model  reduction  procedures  which  maintain  zeros  have  not  been  used  extensively 
except  in  single  input  single  output  systems. 


4. 2. 2.4  Selections  of  Significant  Poles.  The  procedures  of  the  previous  section 
may  be  used  once  poles  to  be  retained  in  the  reduced  order  model  are  specified. 
This  specification  has  usually  been  based  on  two  considerations: 

1 .  Frequency  separation 

2.  Performance  evaluation 

The  frequency  separation  approach  keeps  poles  in  a  prespecified  frequency  range 
depending  on  the  actuator /sensor  bandwidth,  control  requirements  and  modeling 
accuracy.  In  most  large  space  structures  control  design  problems,  all  these  con¬ 
siderations  lead  to  retention  of  low  frequency  poles.  One  exception  is  the  elimina¬ 
tion  of  forced  vibration  at  specified  frequencies  (e.g.,  motors). 

The  performance  evaluation  approach  evaluates  the  effect  of  elimination  of  each 
combinations  of  modes  on  the  performance  of  the  closed  loop  system.  Significant 
work  has  been  done  on  this  problem  by  Skelton  [Sec.  5.1.6] .  The  concept  of  open 
loop  modal  cost  analysis  provides  an  integrated  framework  for  selecting  modes 
which  will  effect  the  feedback  control  law. 
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The  selection  of  the  poles  to  be  eliminated  will  consist  of  two  steps.  In  the  first 
step,  all  poles  and  zeros  not  in  the  frequency  range  of  interest  will  be  eliminated. 
In  the  second  step,  modal  cost  analysis  will  be  used  to  selectively  eliminate  those 
modes  which  have  minimum  effect  on  performance. 

With  the  simplified  model  for  large  space  structures,  we  are  now  ready  for 
control  design. 
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4.2.3  Control  Design 


4. 2. 3.1  Control  Approach.  Design  of  feedback  control  laws  for  large  space 
structures  is  dictated  by  the  following  considerations: 


(i)  Model:  Because  of  modeling  accuracy  as  well  as  computational  consider¬ 
ations,  the  model  is  of  low  order.  The  neglected  modes  have  low  damping 
and  can  be  unstable  if  not  properly  considered  in  control  design. 

(ii)  Sensors  and  Actuators:  Though  structures  have  significant  response 
up  to  infinite  frequency,  physical  actuators /sensors  have  finite  bandwidth. 

In  addition  actuators  and  sensors  have  dead  zeros,  hysterisis  and 
non-linearities . 

(iii)  Model  Uncertainties:  The  behavior  of  structures  is  known  only  approxi¬ 
mately  at  high  frequency.  Therefore,  it  is  difficult  to  assess  the  behavior 

of  high  frequency  poles. 

(iv)  Space  Application:  Because  of  weight  and  volume  considerations  in 
space  applications,  optimization  of  sensors  and  actuators  is  critical  in  control 
of  space  structures. 


The  control  design  approach  developed  addresses  these  issues  in  the  feed¬ 
back  law  specification.  The  integrated  approach  is  based  on  an  extension  of  the 
work  of  Aubrun  [sec  4.1]  on  collocated  sensors /actuators  and  on  frequency 
shaping  filter /controller  design  methods. 


The  collocated  sensor /actuator  based  control  design  approach  has  been  called  the 
low-authority  structural  control  (LASC)  because  of  the  limited  amount  of  structural 
damping  which  it  can  produce.  The  LASC  theory  has  been  generalized  to  include 
many  sensor  /actuator  pairs  with  possible  cross  feedback  among  pairs. 

The  general  feedback  controller /filter  design  problem  considers  general  sensor/ 
actuator  location  with  general  signal  flow  structure.  The  validity  of  the  model 
over  a  small  frequency  range  is  used  directly  in  the  control  design  approach. 


The  resulting  controller  provides  maximum  robustness  against  truncated  modes.  ] 

High  damping  ratios  can  be  obtained  in  selected  modes.  The  stability  of  unmodeled  i 

i 

modes  cannot  be  guaranteed  even  with  perfect  actuators  and  sensors. 

t 

4. 2. 3. 2  Generalization  of  Colocated  Control.  Collocated  control  is  generalized 
for  multi-input  multi-output  systems.  It  is  shown  that  the  solution  is  closely  \ 

related  to  an  inverse  optimal  control  solution.  A  design  procedure  is  discussed 
followed  by  controller  robustness  properties. 

Inverse  Optimal  Control  Problem 

This  section  considers  the  control  problem  with  collocated  actuators  and  sensors. 

In  addition,  all  sensor /actuator  pairs  are  complimentary  to  each  other  (e.g. , 
velocity  sensor /force  actuator,  rate  sensor /moment  actuator).  The  output  vector, 
y,  is  obtained  then  through 

y (t)  =  GT  x(t) ,  yeRm  (1) 

Neglecting  the  low  natural  damping  of  the  flexible  system,  the  finite  order*  linear 
oscillatory  system  is  described  through 

x(t)  =  F  x(t)  +  G  u(t) ,  xeR2n  ,  ueR™  (2) 

where 


and  A  is  an  n  x  n  symmetric  positive  definite  (p .  d . )  matrix ,  and  Bq  is  an  n  x  m 
matrix  of  full  rank  m  (m<n). 


•Refer  to  chapter  3  for  further  details  on  model  derivation. 


The  state  vector,  x,  of  the  system  described  in  (2)  is  defined  by 


x  4  [q,  q]  (4) 

where  q  and  d[  are  generalized  displacements  and  velocities  respectively.  Using 
the  output  of  ( 1)  in  implementing  an  output  feedback  control  law  will  result  in 
a  velocity  (rate)  feedback  only  due  to  the  structure  of  G  in  (3).  This  is  an  im¬ 
portant  property  of  collocated  controls  for  systems  described  as  in  ( 2) ;  it  will 
be  exploited  in  deriving  the  analytical  foundations  for  the  design  procedure. 

Consider  now  the  case  where  the  system  in  (1)  and  (2)  is  a  single-input, 
single-output  (SISO)  system.  Then,  it  can  be  shown  that  the  use  of  a  collocated 
reate  feedback  control  will  result  in  an  eigenvalue  shift  to  the  left  in  the  complex 
plane  (increased  damping) .  The  proof  of  this  property  is  based  on  Jacobi's  for- 

y  * 

mula  for  eigenvalues  perturbations  and  the  modal  structure  of  the  system  in  (2) 
and  (3)  (see  Sec.  4.1). 

Can  collocated  rate  feedback  control  provide  increased  damping  in  multi-input, 
multi-output  (MIMO)  systems  too? 

The  answer  is  provided  through  the  application  of  inverse  optimal  control  theory. 

We  consider  an  objective  function  quadratic  in  state  and  control  over  an  infinite 
time  period ,  i .  e . , 

J(x,  u)  =  f°°  [xT(t)Qx(t)  +  uT(t)Ru(t)]dt  (5) 


and 

Q>0,  R > 0  (6) 

The  problem  of  inverse  optimal  control  to  be  considered  is  the  following.  Given 
a  stabilizing  output  feedback  control  law  described  by 
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(7) 


u(t)  =  Ky  (t)  =  KGTx(t) 

for  what  weighting  matrices  Q  and  R  in  (5)  is  this  control  also  optimal? 

Appendix  A  derives  the  solution  to  this  problem;  a  summary  of  results  is  provided 
in  Table  1. 

The  solution  to  the  inverse  optimal  control  problem  as  formulated  above  will  allow 
us  to  generalize  the  stabilization  properties  of  a  collocated  rate  feedback  control 
to  MIMO  systems.  This  is  discussed  next. 

Properties  of  Solution 


The  stabilization  of  a  linear  system  by  a  state  feedback  control  law  requires  the 
system  in  (2)  to  be  controllable*  [5] .  This  is  an  important  structural  property 
that  will  be  assumed  for  the  system  under  consideration. 

T 

It  is  well  known  [  5]  that  if  Q  =  L  L,  where  (F,  L)  is  an  observable  pair,  the 
optimal  control  obtained  through  the  solution  of  the  ARE  yields  a  stable  closed- 
loop  system.  In  the  case  discussed  in  Appendix  A  we  have 

Q  =  LTL  (8) 

where  the  m  x  2n  matrix  L  is  given  by 

L  =  [o!aB"1/2  b£]  (9) 

It  is  simple  to  show  that  if  (F,  G)  is  a  controllable  pair,  R  and  Aq  are  symmetric 
p.d  matrices  anda#0,  then  (F,  L)  is  an  observable  pair  too.  This  shows  that 
selecting  any  symmetric  negative  definite  matrix  for  the  output  feedback  gain 
matrix  will  result  in  a  stabilizing  collocated  control.  In  addition,  this  control  is 
also  optimal  for  a  specific  optimization  problem  as  discussed  in  Appendix  A.  This 
result  provides  the  required  generalization  to  MIMO  systems. 

*This  requirement  can  be  relaxed  to  that  of  stabilizability  [  5] . 

+Algebraic  Ricatti  Equation. 
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Table  1 
SUMMARY 


Given 


x(t)  =  Fx(t)  +  Gu(t)  ,  xeR2n  ,  ueRm 

u(t)  =  KG^  x  (t)  (rate  feedback) 

T 

K  =  K  <  0  (a  given  gain  matrix) 


Problem : 


For  what  Q  and  R  matrices  is  u  the  optimal  control  for 
linear  quadratic  objective  function  given  by 


J*  =  min  /“  (xT(t)Qx(t)  +  uT(t)Ru(t))  dt 

U€  ° 

Solution  of  ARE:^ 


P  = 


Input  Penalty  Matrix 
R 

State  Penalty  Matrix: 


A  l  0 

_  o_,  _  _ 

0  .  I 

I  n 


=  -aK  x  ,  a>  0 


Q  = 


L°i  _  o  _  J 

f°  'Qf2B  R_1BT  | 
L  o  o  J 


Optimal  Cost: 

=  *oPxo  =  “  (\Vo  +  I  2) 

Algebraic  Riccati  Equation 


The  results  presented  in  Table  1  can  be  used  to  yield  a  physical  interpretation 
of  the  solution  to  the  inverse  optimal  control  problem. 


The  optimal  performance  is  written  as 

J*  =  min[J1  +  JJ  (.10) 

u 


where 


J.  =  a2  /■"  qTB  R-1BT  q  dt  (11) 

1  Jo  o  o  M 

J  2  =  [°° uTRu  dt  (12) 

■'o 


The  interpretation  suggested  above  is  that  the  cost  functional  is  the  total  energy 
of  the  system.  The  term  contributed  by  the  state  is  the  kinetic  energy  where 

E  =  qTMq  (13) 


and  the  "mass,"  M,  is  given  by 

M  =  a2B  R_1BT  (14) 

o  o 

Also,  the  use  of  velocity  feedback  where  the  actuators  and  sensors  are  collocated 
results  in  an  L<5  problem  where  only  kinetic  energy  is  weighted. 

The  solution  of  the  inverse  optimal  control  problem  as  presented  above  has  a 
theoretical  and  practical  significance.  On  the  theoretical  side  we  have  been  able 
to  show  that  this  solution  results  in  a  stable  closed-loop  system.  This  solution  can 
also  be  used  to  answer  the  question  of  control  "spillover".  The  practical  aspect 
of  the  inverse  optimal  control  problem  lends  itself  to  an  implementable  design 
procedure  discussed  in  the  following  section. 


Design  Procedure 


It  was  shown  above  that  the  solution  of  the  inverse  optimal  control  problem  results 
in  a  state  weighting  matrix  given  by 

«  “  p?f— <15) 

0  a  a  R  c 1 
L  1  o  oj 

In  practice,  however,  one  starts  the  design  procedure  by  specifying  the  weighting 
matrices  rather  than  obtaining  them  as  an  end  result.  This  specification  results 
from  the  performance  levels  the  system  is  required  to  achieve. 

The  use  of  low-authority  controllers  (collocated  sensors  and  actuators)  to-date 
was  concerned  mainly  with  the  question  of  stability.  The  results  obtained  above 
allow  us  to  extend  the  scope  of  these  controllers  beyond  that  of  stability. 

Let  the  state  weighting  matrix  Q ,  be  specified  as 


'o  ;  o 


where  Qq  is  any  n  x  n  symmetric  p.s.d  matrix. 


Given  Q  as  in  (16),  the  design  is  concerned  with  solving 


-aB  kb;  =  <3 
o  o  o 


for  the  symmetric  gain  matrix  K. 


If  the  n  x  m  matrix  Bq  is  square  and  nonsingular,  i.e. ,  n  =  m,  then  the  required 


gain  is  obtained  from  ( 17)  as 


K  -  -&)  n;1  % 


1 


In  practice,  however,  one  finds  that  men  in  which  case  (18)  will  not  be 
available.  A  valid  assumption  that  can  still  be  made  is  that  Bq  is  of  full  rank 
equal  to  m. 

In  this  case  we  can  solve  ( 17) ,  so  that  the  solution  is  approximate  in  the  least 
squares  sense. 

Let 

1  =  <r[K  +«B0KB^)  (q0  ^B0KbJ)T]  (19) 

then,  setting 

a  j 
a  k 

we  find 

K>  =  -(s-)(Bo  Bo)_1  BXBo(Bo  Bo)_1  <2» 


=0  (20) 
K=K* 


Remarks : 


1.  Since  Q0  is  p.d,  and  a  >70,  the  gain  matrix  obtained  in  (21)  is  symmetric, 
negative  definite  as  required. 

2.  Since  the  objective  function  of  (19)  is  convex,  the  gain  matrix  of  (21) 
is  globally  optimal. 


Using  the  gain  as  given  in  (21)  for  the  control  of  the  system  in  (2)  through 
the  application  of  (7)  will  result  in  a  solution  of  an  inverse  optimal  control  prob¬ 
lem  where  the  actual  state  weighting  matrix  is  given  by 


A 

Q 


(22) 


and 


Q*  =  Bo(Bo  Bc)‘lBXBo(BJ Bo)'1bJ  (23) 

The  actual  weighting  matrix  9*  of  (23)  can  be  rewritten  as 

Q*  =  CQqC  (24) 

where 

C  =  Bo(BoBo)‘1eo  -  cT  (25) 

The  design  procedure  is  concerned  with  specifying  Qq  so  that  the  resulting  Q* 

as  given  in  (24)  will  approach  a  desired  weighting  structure  i.e. , 

Q*  *  Qd  (26) 

The  above  discussion  can  be  summarized  into  a  numerical  procedure. 

The  numerical  procedure  is  motivated  by  equation  (26).  The  matrix  Qq  is  the 
design  matrix  to  be  chosen  in  an  iterative  manner  until  system's  closed-loop 
performance  is  acceptable. 
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Step  4:  If  eigenvalue  distribution  is  acceptable  (damping  ratios,  frequencies, 
etc . )  go  to  5  ,  otherwise  find  the  diagonal  difference  of  Q  *  and  Qk 

y.  =  Q*(i,i)  -Qk  (i,i>,  l<i<n 

then 

Qk+i  (M)  =  Qk  (i.i)  -£.sgn(y.),  f.>0 

k  — ►  k  +  1 
go  to  3  . 

Step  5:  STOP. 

Remark :  The  choice  of  the  parameter  £.  can  be  varying  depending  present  eigen¬ 
value  distribution. 

This  numerical  procedure  will  be  used  in  section  5  to  design  a  controller  for  a 
flexible  structure  (a  pyramid). 

Robustness  Properties  of  Generalized  Collocated  Control 

A  problem  associated  with  control  design  for  large  space  structures  is  that  of 
truncated  modes.  The  dimension  of  the  model  described  in  Eq.  (2)  depends  on 
the  fineness  of  grid  used  in  the  application  of  the  finite  element  method.  Usually, 
however,  the  dimension  is  quite  large  and  some  reduction  in  order  will  result 
in  modes  that  are  truncated  from  the  final  model  used  for  design.  The  effect  of 
the  control  derived  from  the  reduced-order  model  on  the  truncated  modes  ("control 
spillover,"  Ref.  4  in  sec.  2)  is  of  extreme  importance  in  the  overall  evaluation  of 
the  control  system. 

Another  problem  to  be  considered  is  the  effect  of  parameter  variations  (or  uncer- 
tanity)  on  performance  quality.  The  solution  to  the  inverse  optimal  control  problem 
will  be  used  to  demonstrate  design  robustness  in  regard  to  the  following  two 
specific  problems. 


Problem  1: 


To  what  extent,  if  any,  is  the  solution  to  the  inverse  optimal  control  problem 
sensitive  to  truncated  modes? 


Problem  2: 

How  are  closed-loop  stability  properties  affected  by  parameter  variations 
(or  uncertainty)  in  system  matrices? 


The  problem  of  truncated  modes  will  be  considered  first.  The  matrix  describing 
the  dynamics  of  the  flexible  system,  given  by  Eq.  (2)  has  eigenvalues  distrib¬ 
uted  on  the  imaginary  axis.  The  process  of  mode  truncation  is  carried  out  by 
ignoring  eigenvalues  that  are  far  from  the  origin  on  the  imaginary  axis  (high- 
frequencies) .  Analytically,  the  problem  of  mode  truncation  is  as  follows.  Let 
the  full-order  finite  linear  system*  be  described  by 


z(t) 
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where  the  subscript  T  stands  for  the  part  of  the  system  to  be  truncated.  The 
output  used  for  control  is  given  by 

y (t)  =  B^x(t)  +  B^xT(t) ,  ye  Rm  (28) 


The  matrix  AT  contributes  to  the  full  order  system  frequencies  higher  than  those 
contributed  by  Aq  and ,  therefore ,  can  be  truncated  to  yield  the  reduced  order 
model  described  by 


x(t) 


(29) 


*The  real  flexible  system  is  described  through  a  second  order  PDE  (wave  equation). 
In  linear  form,  the  flexible  system  is  of  infinite  dimension.  Equation  (1),  there¬ 
fore,  considers  only  a  finite  number  of  these  modes. 


Robustness  Property  1: 

The  full  order  system  described  in  Eq.  (27)  under  the  control 

u(t)  =  Ky (t) ,  -K  =  -KT  >  0  (30) 


and  y  given  by  Eq.  (28)  will  remain  stable  for  every  AT  provided  that 
(Ao,  Bq),  (At>  Bt)  are  controllable  pairs. 

The  proof  of  this  property  is  by  using  the  control  as  derived  from  the  reduced 
order  system  and  noting  that  the  full  order  system  assumes  the  form  given  by 
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Appendix  B  indicates  that  this  system  will  be  asymptotically  stable  provided 
(Aq,  Bq)  and  (AT>  B^,)  form  controllable  pairs. 

This  result  shows  how  the  inverse  optimal  control  solution  can  be  put  to  use  to 
answer  the  question  of  control  spillover  (Ref.  4  in  sec.  2). 

Observing  the  proof  of  Appendix  B  we  can  arrive  at  the  second  robustness  result. 
This  result  sljows  how  stability  properties  of  the  closed-loop  system  are  affected 
by  parameter  variations  in  the  system  matrix  Aq. 

Robustness  Property  2: 

Consider  the  dynamic  system  described  in  Eq.  (28)  and  where  the  n  x  n  matrix 
Aq  can  have  uncertain  parameters  that  belong  to  certain  set.  Applying  the  output 
feedback  control  of  Eq.  (4)  will  result  in  a  stable  closed-loop  system  for  every 
Aq  in  the  give  range  of  possible  parameter  variation  provided  AQ  is  always 
symmetric  and  p.d.  over  its  entire  range. 
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Remark :  In  the  case  when  Aq  is  a  diagonal  n  x  n  matrix,  the  above  result 
guarantees  stability  of  the  closed-loop  system  provided  the  diagonal  of  AQ 
remains  positive  over  all  possible  variations. 

4. 2. 3. 3  General  Controller /Filter  Design.  With  collocated  sensors  and  actuators 
output  feedback  may  be  used  to  increase  the  damping  ratio  of  all  modes  up  to  a 
certain  level.  In  many  systems  it  is  not  possible  to  place  actuators  and  sensors 
at  the  same  place.  Often  also  higher  damping  ratio  is  desired  in  certain 
selected  modes.  General  controller /filter  design  theory  is  the  desired. 

Design  Procedure 

Since  the  models  are  of  high  order  with  multi-inputs  and  -outputs ,  I, inear-Quadratic- 
Gaussian  LQG  techniques  are  useful.  LQG  methods  must  be  modified  such  that 
they  account  for  the  following  factors. 

(i)  Model  unvalidity  at  high  frequencies 

(ii)  Bandlimited  actuators /sensors ,  and 

(iii)  Control  requirements  at  low  frequency 

This  section  describes  a  completely  new  theoretical  development  to  extend  LQG 
methods  for  large  space  structure  applications. 

LQG  techniques  are  based  on  state  description  models  of 

x  =  Fx  +  Gu  (31) 

y  =  Hx 

The  function  to  be  optimized  is  quadratic  in  states  x  and  input  u.  A  direct  optimi¬ 
zation  of  the  performance  index  gives  the  desired  control  law.  These  methods 
require  many  assumptions. 


1.  The  model  is  valid  for  all  values  of  inputs  and  states.  In  addition,  the 
dynamics  is  well  described  at  all  input  and  state  frequencies. 

2.  The  filter  uses  values  of  F,  G  and  H  matrices  explicitly  in  the  implementa¬ 
tion  stage  (in  addition  to  the  gain  K  which  also  depends  on  the  state 
definition  matrices).  The  filter  design  also  assumes  that  the  dynamics  is 
known  equally  accurately  at  all  frequencies.  This  may  make  the  combina¬ 
tion  of  the  filter  and  the  control  law  extremely  sensitive  to  errors,  see 
also  Doyle  [2]. 

3.  The  optimality  of  the  filter  is  strongly  dependent  on  the  accuracy  of 
noise  statistics. 

Classical  control  design  methods  are  based  on  frequency  domain  descriptions  of 
systems  and  account  for  model  uncertainty  at  high  frequencies  by  the  use  of  gain 
and  phase  margin  comcepts.  Gain  and  phase  margins  provide  a  measure  of  uncer¬ 
tainty  in  the  model  at  high  frequency  for  which  stability  can  be  guaranteed. 

It  has  been  shown  by  Anderson  and  Moore  [5]  and  Athans  and  Safonov  [6]  that 
the  LQG  controller  has  a  60°  phase  margin  and  50%  to  infinite  gain  margin.  The 
phase  and  gain  margin  properties  provide  for  a  constant  phase  error  in  all  chan¬ 
nels  and  for  individual  variations  in  gain.  The  system  may  be  extremely  sensitive 
if  two  gains  change  in  opposite  directions.  In  general,  the  phase  and  gain  margin 
property  does  not  relate  directly  to  parameter  sensitivity.  The  performance  of 
the  LQG  controllers  degrades  further  when  a  filter  is  used  for  state  estimation. 

The  gain  and  phase  margin  properties  are  no  longer  valid.  In  addition,  the  filter 
dynamics,  dictated  by  specified  noise  characteristics,  may  be  too  fast,  leading  to 
interaction  with  unmodeled  terms. 

To  understand  the  concept  of  frequency  shaping,  it  is  necessary  to  write  the 
standard  LQG  cost  functional  in  the  frequency  domain.  With  inifinte  time  horizon 
and  no  weighting  on  the  final  state,  the  cost  functional  may  be  written  as: 

J  =  1/2/  [x*(icJ>Ax(jw)  +  u*(ju>)Bu(ja>)]du>  (32) 

00 

where  *  implies  complex  conjugate.  Clearly,  in  this  formulation,  the  weighting 
matrices  are  not  functions  of  frequency,  i.e. ,  the  state  and  control  excursions 
at  all  frequencies  are  equally  unacceptable.  In  many  systems,  on  the  other  hand, 
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inputs  in  the  neighborhood  of  a  particular  frequency  are  not  desirable  because 
of  poor  sensor  or  actuator  characteristics  at  that  frequency.  Historically,  this 
constraint  of  constant  weighting  at  all  frequencies  has  resulted  because  of  the 
difficulty  of  shaping  the  weighting  functional  with  frequency  in  the  conventional 
time  domain  LQG  formulation.  Representation  of  the  cost  functional,  Eq.  (32), 
provides  a  clue  to  the  frequency  domain  approach  for  using  frequency  shaping 
ideas  in  modern  control  theory  techniques.  Matrices  A  and  B  in  Eq.  (32)  may  be 
made  functions  of  frequency  to  give  a  generalized  cost  function  of  the  form 


=  1/2  /  [x*(jcti)A  (ja;)x(ju>>  +  u*(ju))B(ja>>  u(ju>)]dw  (32) 

OO 


A(joj)  and  B (jco)  are  Hermitian  matrices  at  all  frequencies. 


A  design  procedure  based  on  this  weighting  is  given  in  Appendix  C . 

The  frequency  shaped  weighting  matrices  are  used  in  large  space  structure  control 
as  follows.  Maximum  control  activity  is  desired  at  low  frequency  with  minimum  on  no 
control  activity  at  high  frequency.  This  may  be  achieved,  for  example,  by  select¬ 
ing  B (jco)  small  at  low  frequency  and  large  at  high  frequencies.  If  the  model  is  not 
valid  at  high  frequency,  we  may  also  want  to  minimize  state  excursions  in  that  region. 
This  may  require  large  state  weighting  at  high  frequency.  Table  2  gives  a  set 
of  weighting  functions  which  may  be  useful  in  space  structure  control.  As  shown 
in  Appendix  C  and  with  examples  in  Section  5,  such  weighting  functions  lead  to 
general  compensation  structures  in  the  feedback  loop. 

Filter  Design 

States  required  for  feedback  control  are  often  not  measured  directly.  Reliable 
estimators  should  then  be  designed  to  extract  the  state  time  history  from  the  meas¬ 
ured  data.  In  general  states  corresponding  to  the  modes  retained  in  the  model  must 
be  estimated  with  minimium  contamination  from  high  frequency  model  states. 
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A  USEFUL  SET  OF  FREQUENCY  SHAPED  WEIGHTING  FUNCTIONS 


INTEGRAL  CONTROL 


Skelton  has  developed  an  orthogonal  filter  approach  for  state  estimation  in  systems 
with  unknown  non- white  noise.  Consider  the  dyanmic  system 


x  =  Fx  +  Gu  +  w  (34) 

where  w  is  a  vector  of  nonwhite  and  nonstationary  noise.  The  concept  involves 
approximating  w  by  a  vector  of  onthogonal  functions ,  E 

w  =  n/e  (35) 

Chekychev  polynomials  could  be  selected  as  the  set  of  approximating  orthogonal 
function. 

In  Kalman  filter  or  observer  formulation ,  it  is  advantageous  to  select  approximat¬ 
ing  functions  which  can  be  generated  by  a  set  of  differential  equations.  Fourier 
series  and  orthogonal  polynomials  both  satisfy  this  requirement ,  E  then  follows 
a  differential  equation 

E  =  F'E  +  W'  (36) 

where  W  is  white  noise  source  introduced  to  account  for  changes  in  the  error  terms. 

A  nonwhite  error  term  may  also  be  introduced  in  the  measurement  equation 

y  =  Hx  +  v  (37) 

and  v  may  be  approximated  by  the  same  set  of  orthogonal  functions  and  a  pure 
white  noise  term  ’v* 

v  =  DE  +  v'  (38) 

Equations  (34)  through  (38)  may  be  formulated  into  a  Kalman  filter  for  estimat¬ 
ing  the  states.  As  a  by  product,  the  orthogonal  error  knows  E  will  also  be  estimated. 
If  necessary  r  and  D  may  also  be  estimated  in  realtime.  This  method  is  described 
later  in  greater  detail. 


Another  approach  for  state  estimation  is  based  on  frequency  shaping  of  state  and 
measurement  noise.  In  Eq.  (34)  ,  for  example,  noise  W  is  primary  high  frequency 
corresponding  to  truncated  modes.  The  noise  spectrum  is  then  selected  to  be  a 
function  which  is  small  in  the  model  frequency  range  and  high  in  regions  of  trun¬ 
cated  modes.  The  measurement  noise  is  treated  in  a  similar  manner. 

The  negative  log  likelihood  function  for  state  estimation  becomes 

OO  v 

NLLF  =  J  (w*  Q~1(ju>>  W  +  v*  R~1(]w)v)  da) 

OO 

Much  work  has  been  done  on  guidelines  for  selecting  the  weighting  functions. 

Once  appropriate  shaping  functions  have  been  selected ,  the  filter  design  procedure 
is  similar  to  the  control  design  procedure  outlined  in  Appendix  C. 

Relationship  Between  Orthogonal  Filter  and  Frequency  Shaping  Methods 

Orthogonal  filters  and  frequency  shaping  state  estimation  methods  try  to  correct 
for  inadequacy  of  the  model  in  certain  frequency  ranges.  The  orthogonal  filter 
approach  models  the  error  by  a  series  of  orthogonal  functions  in  time  domain 
while  the  frequency  shaping  approach  approximates  the  error  by  a  series  of 
functions  in  the  frequency  domain.  Both  methods  require  additional  states  to 
model  the  nonwhite  error.  If  the  set  of  orthogonal  functions  and  the  frequency 
shaped  weights  are  appropriately  selected ,  the  two  methods  should  give  identical 
results.  Figure  2  compares  the  frequency  domain  and  time  domain  algorithms. 

4. 2. 3. 4  Robustness  Properties.  In  analyzing  the  effect  of  control  on  unmodeled 
modes  we  need  to  look  at  their  motion  in  the  complex  plane.  Let  us  introduce  a 
scaling  parameter  Y  in  the  feedback  loop.  For  small  Y,  a  pole  at  ck  moves  from 
open  loop  value  at  the  rate 


Frequency  Weighted  Cost  Functional  and  Orthogonal  Filter  Approaches  to  Estimator  Design 
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where  R.  is  the  residue  of  the  pole  in  the  open  loop  system.  If  the  reduced  model 
maintains  residues  of  the  modeled  modes ,  the  small  gain  behavior  of  retained  modes 
is  not  affected  by  model  reduction.  The  behavior  of  the  unmodeled  modes  is  dictated 
by  the  residue  at  the  pole  (a  property  of  the  system  and  actuator  /sensor  locations) 
and  the  controller  transfer  function  at  the  pole  (a  property  of  the  controller).  The 
sign  of  the  residue  is  often  not  known  precisely  but  its  magnitude  can  be  bounded. 
This  provides  a  bound  on  the  rate  at  which  each  eigenvalue  moves  for  small  gains. 

The  high  gain  properties  of  the  closed  loop  system  are  difficult  to  analyze.  All 
poles  will  ultimately  move  into  a  zero  or  go  to  infinity.  The  closed-loop  behavior 
is  extremely  sensitive  to  zero  locations  whose  computed  position  depends ,  in  a 
sensitive  way,  on  the  modal  functions  selected  for  the  truncated  model. 

Frequency  shaping  approach  ensures  robustness  of  unmodeled  modes  by  designing 
control  laws  such  that  the  rate  of  movement  of  unmodeled  modes  is  much  smaller 
compared  to  the  rate  at  which  modeled  poles  move.  This  is  achieved  by  designing 
controllers  such  that  Tc(s)  is  high  at  low  frequency  and  small  in  the  high  fre¬ 
quency  region.  This  is  realized  by  frequency-shaped  weighting  functions  in  which 
the  weights  increase  with  frequency.  If,  in  addition,  the  low  frequency  poles  do 
not  move  in  the  high  frequency  region,  the  residues  Rj  on  unmodeled  modes  will 
not  change  significantly.  If  the  residue  of  the  ith  mode  is  of  magnitude  smaller 
than  | Rj  |  ,  the  closed-loop  pole  a.  will  be  contained  within  a  circle  with  center  at 
a.  and  radius  |  R- 1  abs  ^  X  — or . )  ^  .  To  guarantee  stability,  this  circle  must  lie 
in  the  left  half  plane. 
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4.2.4  Orthogonal  Filters 


Introduction 


The  orthogonal  filter  seeks  to  compensate  for  errors  in  the  mathematical  model 
of  the  structure  in  the  state  estimation  problem.  The  potential  of  the  method 
lies  in  the  fact  that  the  total  model  error  is  estimated  rather  than  just  the  para¬ 
meter  errors  available  in  present  parameter  adaptive  and  identification  methods. 
First  we  show  that  modeling  errors  can  always  be  decomposed  into  (1)  parameter 
errors,  (2)  truncated  modes,  (3)  disturbances,  and  (4)  nonlinearities.  Then 
the  idea  of  the  orthogonal  filter  is  presented. 

Model  Errors  for  Matrix- Second-Order  Systems 


Let  the  nonlinear  evaluation  model  of  the  structure  be  described  by 


+  (D1  +  Gl)qi  +  Klql  =  B1U  +  Qy  u>  t)  +  lw1 

N1  m 

Yj  =  Pi<Ii  +  Riqi  +  Aiqi  Qj  £  R  ,  ueR  (1) 

zi  =  piqi  +  Riqi  +  Aiqi  +  g(qr  qr  °  +  vi’  zi  eR* 


where  q^  has  dimension  N^,  and  g(q^,  q^,  t)  represents  any  nonlinearities 

in  the  measurenlent  equations.  Parameters  of  the  matrices  P^,  R and  A ^ 

are  fixed  by  the  location  of  displacement  sensors,  rate  sensors,  and  accelerometers, 

k 

respectively,  y^  e  R  represents  the  k  variables  we  wish  to  control.  For 

example,  if  we  wish  to  control  the  entire  state  vector  then  y  becomes 

T*  T  T  ^ 

yi  =  (ql’  q^)  by  aPProPriate  choice  of  P^,  Rj,  Aj  . 


P 


1 


(2) 
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The  performance  objectives  might  be  stated  in  terms  of  y^(t)  and  u(t),  as 

in  Ref  2,  or  they  might  be  stated  in  terms  of  closed  loop  eigenvalue  locations, 

as  in  Ref  3.  This  report  is  not  specifically  concerned  with  the  control  policy. 

Instead,  it  is  concerned  with  the  question  of  whether  the  model  error  vectors 

e  (t),  e  (t),  or  e  (t)  associated  with  any  design  model  S9 
q  y  z  i 

n2 

^2^2  +  ^2  +  ^2^2  +  ^2^2  =  ®2U  +  eq’  ^  £  R 

y  1  =  ^2^2  +  ^2^2  +  ^2^2  +  6y  ’  U  e  ^  (3) 

k  t 

Z1  =  P2^2  +  ^2q2  +  ^2q2  +  ez  ’  yi  £  R  >  Z1  e  R 

T  T  T 

can  be  reconstructed  on-line  at  the  same  time  the  state  vector  x2  =  (q2,  q2> 
is  reconstructed  from  the  measurements  z^(t)  .  (By  definition  of  the  evaluation 
model  (1)  it  is  assumed  throughout  that  z^(t)  represents  the  actual  measure¬ 
ments  sufficiently  accurately  for  controller  evaluation  purposes . ) 

No  k  . 

Model  error  vectors  e  (t)  t  R  ,  e  ( t )  e  R  ,  e  (t)  e  R1  exist  such  that  the 

q  y  z 

output  vectors  y^t),  and  z^(t)  of  the  reduced  order  model  (3)  are  identically 
equal  to  the  corresponding  output  vectors  y^(t)  and  z^(t)  of  the  higher  order 
model  (1).  The  model  error  vectors  are 


where  e  is  that  contribution  to  the  total  model  errors  due  to  parameter  errors 

in  (M2>  D2,  G2,  K2,  P2>  R2,  A2,  P2,  r2,  >12),  et  is  that  contribution  to  the 

total  model  errors  due  to  errors  in  model  order  (truncated  modes) ,  e  is  that 

w 

contribution  due  to  erroneously  modeled  disturbances,  and  e^L  is  that  con¬ 
tribution  to  the  total  model  errors  due  to  neglected  nonlinearities.  Furthermore, 
eA  ’  et’  V  an<i  eNL  ^ave  exPlicit  expressions 
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w 


mm 


The  transformation  matrix  Q-  [  n  ,  fi2]  may  be  any  nonsingular  matrix,  and  the 
parameter  errors  (AIVI^  AD2>  AG^  AK2,  AB2>  AP2>  AR^  AA2,  AP2 ,  A/?2 ,  AA2> 
associated  with  the  reduced  model  (3)  will  be  defined  by 


AP2 

— 

Plfll  "  P2 

ap2 

= 

P1  ni  "  P2 

ar2 

= 

R1  ~  R2 

A  «2 

= 

Vr  p2 

AA2 

= 

A1  Q1  "  A2 

AA2 

= 

AlV*2 

am2 

- 

nIM,V 

AD2 

= 

AG2 

= 

nTGlBl- 

AK2 

= 

°TKlV 

ae2 

= 

bi  -  b 

(7) 


Selection  of  Orthogonal  Functions  for  Model  Error  Approximation 


Under  these  conditions: 

,(k+l)T 


1. 


/ 


e  (t)  e(t)  dt  <  =o  ,  t  >  o,  k  =  0,  1,  2, 


(8) 


for  any  set  of  functions  satisfying 


2. 


/ 


(k+l)r 


-y  (t)  V  (t)  g(t)  dt  =  a  ,  y  e  rc 


(9) 


where  A  is  nonsingular  and  g(t)  is  a  positive  scalar  weighting. 
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there  exists,  during  any  interval  r  >0,  a  set  of  coefficients  such  that 
convergence  in  the  mean  squared  sense  is  guaranteed. 


lim  l|  e.(t)  -  S  F..  Y.(t) 

d-J1  j=l  13  3 


=  0  ,  te  [kr  ,  (k  +  1)t] 


(10) 


Furthermore,  the  functions  Y^.(t),  j  =  1,  2,  ...  d  can  be  generated  by 

y  =  D 

for  some  choice  of  D  .  Thus,  by  augmenting  the  model 

x  =  Ax  +  Bu  +e  e  =  FY 

x  x 

z  =  Mx  +  v 

with  the  model  error  system 


(11) 


Y  =  DY  ,  e  =  FV 
x 


the  estimator  can  be  constructed 


'a 

F 

(!)  . 

'B' 

.0 

D  - 

.  0  . 

u  +  G*(z  -  Mx) 


(12) 


which  produces  estimates  of  the  state  x  and  the  model  error  vector  e  ,  where 
a  a  x 

e  =  FY. 
x 


The  difference  between  the  construction  of  the  orthogonal  filter  and  the  construction 
of  the  estimator  using  frequency  weighted  costs  is  that  the  former  requires  selection 
of  a  set  of  basis  functions  in  the  time  domain,  whereas  the  latter  requires  selection 
of  a  set  of  basis  functions  in  the  frequency  domain.  Further  research  into  the 
equivalence  of  the  two  points  of  view  is  required. 
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4.2.5  Controller  Simplification 


4. 2. 5.1  Problem  Motivation.  The  modeling  process  for  a  large  flexible  structure 
can  result  in  a  dynamic  system  of  prohibitive  dimensions.  Model  order  reduction 
is  performed  in  two  stages.  The  first  stage,  discussing  earlier,  consists  of  model 
reduction  prior  to  control  design.  Further  reduction  in  controller  complexity  is 
obtained  by  a  direct  simplification  of  the  dynamic  controller  itself.  This  section 
discusses  methods  for  control  law  simplification. 

There  are  two  aspects  of  controller  reduction,  (a)  numerical  procedure,  and 
(b)  criteria  function.  It  is  shown  below  that  control  laws  resulting  from  modern 
control  theory  can  be  perceived  as  dynamic  systems.  Numerical  procedures  for 
reduction  are,  therefore,  similar  to  those  for  model  reduction.  Criteria  functions 
which  may  be  used  for  controller  simplification  are  given  subsequently . 

4. 2. 5. 2  Controller  as  a  Dynamic  System.  The  controller  can  be  viewed  as  a 
dynamic  system  to  which  model  reduction  techniques  can  be  applied.  When  a 
state  estimator  is  used  in  implementing  the  control  law  the  controller  is  de¬ 
scribed  by  the  following  dynamic  equations: 


x  =  Fx  +  Gu  +  K  (y  -  Hx) 

A 

u  =  Cx 


(1) 

(2) 


which  can  be  rewritten  as 


This  shows  that  the  controller  can  be  viewed  as  a  dynamic  system  whose  inputs 
are  the  (output)  measurements  and  whose  outputs  are  the  control  to  be  applied 
to  the  system. 

This  dynamic  system  is  of  the  same  order  as  the  original  system  used  in  control 
design.  One  of  the  two  methods  may  be  used  for  reducing  control  law  complexity 
(a)  reduced  order  observers,  or  (b)  dynamic  system  reduction. 

Reduced  Order  Observers 


The  system  described  in  (1)  is  a  full-order  Kalman  filter  on  Luenberger 
observer.  Therefore,  one  simple  way  to  alleviate  the  possible  dimensionality 
problem  mentioned  above  is  to  use  a  reduced-order  observer.  This  is  described 
by: 


where 


z  = 


TA  -  F  T 
o 


+  K  y  +  G  u 

(5) 

o  J  o 

=  Tx 

(6) 

=  K  H 

(7) 

o 

=  TB 

(8) 

and  z  is  an  (n-m)  dimensional  vector  (y  is  m- dimen sional) , 
Defining 
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'  h  T  [  ; 

;  T j  =  Lp ;  QJ 


the  state  estimate  is  given  by 


(9) 


x  =  Py  +  Qz 


(10) 
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Therefore,  the  reduced-order  controller  is  given  by 


(FQ  +  GqCQ)z  +  (Kq  +  CP)y 

(ID 

CPy  +  CQz 

(12) 

Fq  may  be  chosen  arbitrarily 

The  reduction  of  complexity  was  possible  since  the  information  available  through 
the  output,  y  ,  is  not  reconstructed  in  the  reduced-order  observer.  The  order  j 

of  the  system  cannot  be  reduced  by  more  than  the  number  of  outputs. 

Controller  Simplification 

Let  us  consider  the  controller  through  the  following  general  dynamic  equation 

x  =  F*x  +  Ky  (13) 

u  =  Cx  (14) 

This  dynamic  system  may  have  real  on  complex  poles. 

A  simple  structure  for  a  reduced  controller  is  offered  by 

u  =  Fxu  +  Gjy  (15) 

‘ 

This,  however,  requires  that  ' 

CF*  =  FjC  (16)  * 

G :  =  CB  (17) 

Equation  (17)  presents  no  special  problems  but  (16)  is  not  always  solvable.  It 
can  be  shown  that  a  necessary  and  sufficient  conditon  for  (16)  to  be  solvable 
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is  that  the  gain  matrix  C  be  the  left  eigenvectors  of  the  full  order  observer 
matrix  F*  in  (13).  This  condition  could  perhaps  be  met  approximately  in  many 
systems.  The  benefits  in  doing  so  are  clear  but  how  to  actually  carry  it  out  is 
still  an  open  question. 


4. 2. 5. 3  Criteria  for  Model  Reduction.  In  general,  controller  simplification 
should  be  such  that  the  resulting  control  will  yield  satisfactory  performance 
with  not  much  degradation  compared  to  the  following  criteria  are  available  to 
minimize  loss  in  performance  because  of  controller  simplifications. 

(a)  Residues :  Elimate  controller  poles  with  small  residues 

(b)  Pole- Zero  Cancellation:  Eliminate  poles  which  approximately  cancel 
controller  zeros 

(c)  Frequency  Range:  Approximate  controller  poles  and  zeros  outside  the 
frequency  range  of  interest 

(d)  Modal  Cost  Analysis:  Discard  controller  poles  which  have  minimal  effects 
on  system  performance.  This  methodology  has  been  firmalized  by  modal 
cost  analysis ,  developed  by  Skelton  [  ] .  This  method  considers  a  cost 
functional  given  by 


v  =  "<3  +  11  “  "k  )  =Lvi 

i=l 

the  simplification  is  based  on  eliminating  modes  whose  contribution  to 
the  over  all  cost  (V.)  is  small  compared  to  others.  This  is  the  most 
sophisticated  controller  simplification  method. 


Summary 

This  section  shows  methods  for  controller  simplification.  Since  the  controller  can 
be  treated  as  a  dynamic  system,  controller  simplification  is  very  similar  to  model 
reduction . 
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4.3  SELECTION  OF  ACTUATOR  LOCATIONS 

4.3.1  Introduction 

The  problem  of  actuator /sensor  locations  in  dynamic  systems  has  received  little 
attention,  despite  its  importance.  Often,  actuator  locations  are  selected  before 
the  control  design  or,  are  specified  based  on  other  considerations.  In  these  eases 
the  designer  is  faced  with  the  problem  of  constructing  an  acceptable  control  law 
to  best  satisfy  a  predetermined  objective,  with  the  constraints  of  actuator /sensor 
locations. 

As  increasing  demands  are  placed  on  feedback  control,  it  is  necessary  to  manipulate 
actuator  locations  to  meet  design  requirements. 

The  large  space  structure  has  a  low  natural  damping  and  behaves  almost  like  a 
pure  oscillatory  system.  The  actuators  on  this  system  are  required  to  increase  the 
damping  of  the  naturally  flexible  structure.  It  is  up  to  the  designer  to  determine 
actuator /sensor  locations  that  best  meet  the  design  objective. 

This  section  considers  the  problem  of  optimal  actuator  locations  for  an  oscillatory 
dynamic  system.  It  will  be  shown  that  closed-form  expressions  concerning  the  modal 
structure  of  the  dynamic  system  can  be  obtained.  This  modal  information  can  then 
be  used  in  developing  the  design  procedure  of  this  report. 

4.3.2  Problem  Statement 

The  oscillatory  dynamic  system  to  be  considered  in  this  section  is  described  by 
x(t)  =  Fm  x(t)  +  G  u(t),  xeR2n  ,  ueRm  (1) 
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where 


F  = 


(2) 


and  A  is  a  symmetric,  positive-definite  (p.d)  n  x  n  matrix  I  is  the  n  x  n  identity 
matrix,  and  B  is  an  n  x  m  matrix. 


The  system  matirx,  F,  is  determined  by  the  physics  of  the  problem,  but  the  input 
matrix ,  G ,  can  be  selected  subject  to  some  criteria  and  the  structure  given  by 
Eq.  (2).  The  next  section  will  describe  a  particular  criterion  to  be  employed  as  a 
measure  of  optimality  in  actuator  location,  i.e.  ,  the  specification  of  G.  The  pro¬ 
blem  of  sensor  location,  even  though  not  treated  here,  can  be  similarly  defined 
and  solved  by  . duality. 


4.3.3  Optimality  Criteria 


The  state  vector  of  the  oscillatory  system  given  in  Eq.  (1)  contains  the  displace¬ 
ments  and  velocities  at  various  nodes  of  the  system.  In  placing  the  actuators,  it 
is  desirable  that  control  energy  will  be  such  that  after  some  finite  time  both  dis¬ 
placements  and  velocities  will  be  zero,  i.e.  , 

x(T)  =  0  for  some  finite  T  >  0  (3) 


An  additional  requirement  is  that  the  control  energy  exerted  by  the  actuators  will 
be  minimal. 

Kalman,  Ho,  and  Narendra  [1]  have  shown  that  the  minimal  energy  control^  is 
given  by 

u*(t)  =  GT(t)  p  (to,t)W_1(to,T)  [0(tQ,T)x(T)-xo]  (4) 


where  0(.)  is  the  transition  matrix  and  W(.)  is  the  controllability  matrix  defined 

by 

^For  details,  see  Appendix  D 


(5) 


4.3.4  General  Form  For  The  Controllability  Matrix 

The  n  x  n  matrix  A  in  Eq.  (2)  is  a  symmetric  p.d.  matrix.  As  will  be  demonstrated 
shortly ,  the  transition  matrix  as  well  as  the  controllability  matrix  can  be  written 
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rather  easily  in  closed  form  if  A  is  a  diagonal  matrix. 

Therefore,  the  first  step 

in  the  analysis  will  be  to  define  a  transformation  that  will  result  in  an  equivalent 

system  with  this  property.  Let 

(11) 

where 

’  ■  PrivV] 

(12) 

and  U  is  the  n  x  n  matrix  of  eigenvectors  of  A,  i.e. 

AU  =  U  . 

A 

(13) 

Aa  =  diag  (a1,  ...  ,  an) 

(14) 

Then  the  equivalent  system  is  given  by 

z  =  F  z  +  G  u 
o  o 

(15) 

where 

Fo  -  TAoT‘1'  Go  -  TB„ 

=  [|]  (16) 

and 

G  =  U_1B 

(17) 

Furthermore,  if 

sftnA  F«“V 
0X(t’V  “  e 

(18) 

,t  ,  ,  4 

0z(t’V  “  e 

(19) 
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then 


If 


and 


-1 


»x(t’V  =  T  9z(ut0n 


Wx(0,T)  =  /  0x(O,t)G  GT  0T(O,r)dr 


Wz(0’T)  =  /  0z(O)T)Go  Go  0T(°»T>dT 


(20) 


(21) 


(22) 


then,  using  Eqs.  (15)  and  (20) 

W  (0.T)  =  T-1  W  (0,T)(T-1)T 
A  z 


The  measure  of  optimality  will  now  be  given  by 


min 


tr  MW 


1/k 


where 


M  =  Tk  Mo(TT)k 


(23) 


(24) 


(25) 


The  modal  structure  for  the  system  described  by  (15)  is  derived  in  Appendix  A. 
This  structure  results  in  a  closed  form  expresssion  for  the  transition  matrix  given 

by 


4>z(0,t) 


[“1 

s 

0  1 

s  ! 

Bu 

s 

0 

0 

s*ar  | 

0 

'•B 

---1 

n 

“u 

V 

0 

X  1 

■v 

0 

'a  | 

0 

'a 

_ 

n>. 

n  _ 

(26) 
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where 


a.  =  cos  '/a.  t 

h  =  sin 4 

o*j  =  s/itj  sins/a.  t 


i  =  1 ,  2,  . . . ,  n 


(27) 


Using  the  definition  for  the  controllability  matrix  as  given  in  Eq.  (22),  the 
expression  for  the  transition  matrix  as  given  in  Eqs,  (26)  and  (27)  and,  in 
addition ,  the  following  definition 


T 

G 

o  o 


G11  G12  Gln 


G 


nl 


G 


nn 


(28) 


we  obtain  the  following  expression  for  W  (0,T), 

z 


Wz(0,T)  = 


W 


11 


L  W12 


W 


W 


12 


22  J 


(29) 


The  four  n  x  n  symmetric  blocks  are  described  as  follows: 


Wn  =  1 

O 


Gll/31  G12^1^2~ 
G22/32 


Gln¥n 


G  P 
nn  n  j 


dt 


W12  { 
o 


GllPfl  G 1 2^la2 - Gln¥n 


G22/32Qf2 


G  /3q' 

nn  n  n 


dt  =  W 


21 


Using  the  expression  for  the  controllability  matrix  we  obtain  the  main  result. 


Theorem  1:  For  terminal  time,  T,  large  enough  (  10  sec)  the  controllability  matrix 
given  in  Eq.  (29)  becomes  diagonally  dominant  and  can  be  approximated  by 


W..(0,T) 

li 


G..T 


11 

lliin 

2  . 
ai 

G..T 

11 

2 

n<i£2n 

(30) 


Proof : 


The  diagonal  elements  of  W_(o,T)  are  given  by 

z 


W^O.T)  = 


G 

a 


sm 


2T  \T i.\ 


Gii 


.T  + 


2  A 

sin2T  /a.) 
2  , 


l<i<n 

n<i<2n 


for  T  large  enough  the  second  term  becomes  negligible  and  hence  we  get  Eq.  (30) 
The  off-diagonal  elements  are  a  product  of  sine  and  cosine  functions  only  without 
the  additive  term  of  the  terminal  time  as  found  for  the  diagonal  terms.  Therefore, 
for  T  large  enough  the  matrix  W  becomes  diagonally  dominant  with  terms  given  by 
Eq.  (30). 

With  this  expression  for  the  controllability  matrix  we  can  turn  now  toward  the 
numerical  problem  of  solving  for  the  actuator  locations  Eq.  (7)  or  Eq.  (24). 


4.3.5  Numerical  Procedure 


The  optimal  location  of  actuators  will  be  established  through  the  solution  of 


The  steps  of  the  algorithm  are  as  follows: (also  see  enclosed  flow  chart.  Fig.  11: 

Step  1 :  Read  A  ,  Mq  ,  n ,  m  ,  k  ,  e ,  T 
Step  2:  Find  U,A^  such  that 

AU  =  U  A 

Step  3:  Construct  T  as  in  Eq.  (12)  and  find  M  as  in  Eq.  (25) 

Step  4:  Make  an  initial  guess  for  the  parameters  vector  9. 

Step  5:  Find  the  matrix  B  from  known  analytical  relationship  with  0, 
tables,  etc. 

Step  6:  Find  the  diagonal  terms  of 

-IT  -11 
G2  =  U  BB1  (U  V 

Step  7:  Construct  the  diagonal  of  the  controllability  matrix  W^,,  via  Eq.  (30) 
-k  1/k 

Step  8:  min(tr  MWz  )  to  find  an  update,  to  the  parameters  vectors  9. 

Step  9:  If  the  update  Ad  is  small,  i.e.  ,  ||A0||<e  STOP;  otherwise,  GO  TO 
Step  5. 

Remark:  In  step  3  only  the  diagonal  terms  of  M  have  to  be  found. 


4.3.6  Conclusion 


The  problem  of  actuator  location  for  oscillatory  system  was  considered.  It  was 
shown  that  through  the  controllability  matrix  one  can  arrive  at  an  actuator  location 
that  is  derived  via  the  solution  of  the  minimum  energy  control. 

A  similar  approach  can  be  used  for  the  solution  of  the  sensor  location  problem. 


Section  5 

APPLICATIONS  TO  DARPA  STRAWMEN 


5.1  THE  CHARLES  STARK  DRAPER  LABORATORY  (CSDL)  EXAMPLE 

5.1.1  The  CSDL  Example  Structural  Model 

The  model  illustrated  in  Fig.  1  is  a  conceptual  structure  recently  designed  by 
Dr.  Keto  Soosaar  and  R.  Strunce  of  the  Charles  Stark  Draper  Laboratory  (CSDL). 
The  CSDL  example  consists  of  a  (regular)  tetrahedral  truss  supported  on  a 
ground  by  three  right-angled  bipods.  The  bipod  legs  are  pinned  to  ground,  and 
all  other  nodes  are  clamped.  The  individual  truss  members,  including  the  bipods, 
have  elastic  flexibility  in  their  axial  directions  only,  i.e.  ,  they  can  undergo  com¬ 
pression  and  elongation  but  no  lateral  bending.  The  four  vertices  of  the  tetra¬ 
hedron  each  have  three  degrees  of  freedom,  so  that  the  entire  model  can  be 
completely  represented  by  12  structural  modes. 

The  structural  model  is  described  in  consistent,  but  unspecified  units.  The 
edges  of  the  (regular)  tetrahedron  are  each  10  units  long,  and  the  bipod  legs 
are  each  2-72  units  long.  (The  "horizontal"  base  of  the  tetrahedron  is  therefore 
situated  2  units  above  ground,  i.e.  ,  above  the  x,y  plane.)  The  top  vertex  is 
defined  as  LOS  (line-of- sight) ,  and  the  six  bipod  legs  are  defined  to  act  as 
"member  dampers,"  i.e.  ,  active  axial  spring/dashpots,  each  equipped  with  a 
rate  and  position  sensor.  The  three  bipods  represent  therefore  6  colocated 
actuator /sensor  pairs,  each  pair  being  identified  with  a  bipod  leg.  Sensor  meas¬ 
urements  and  actuator  forces  are  thus  confined  to  the  axial  elongation /compression 
motions  and  rates  of  the  bipod  legs. 

The  stiffness  properties  of  the  structure  are  non-uniform,  and  are  given  in  terms 
of  the  truss  members'  cross-sectional  areas.  Two  sets  of  values  are  given,  one 
for  a  "nominal"  model,  and  the  other  for  a  "perturbed"  model,  as  follows:  (see 
Fig.  1  for  truss  element  numbers) 
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Example  Structural  Model 


Truss  Element  No.  Area  (Nominal  Model)  Area  (Perturbed  Model) 


1 

100 

150 

2 

100 

150 

3 

100 

150 

4 

100 

150 

5 

100 

150 

6 

100 

150 

7 

1000 

1200 

8 

1000 

1200 

9 

1000 

1200 

10 

1000 

1200 

11 

100 

150 

12 

100 

150 

where,  for  both  models,  Young's  modulus  is  E  =  1.0.  Finally,  there  are 
lumped  masses  of  2  units  at  each  vertex  of  the  tetrahedron  for  the  nominal  model. 
For  the  perturbed  model,  the  lumped  mass  at  the  LOS  vertex  (see  Fig.  l )  is 
changed  to  4  units,  and  the  others  remain  the  same. 


The  structural  frequencies  (open-loop  eigenvalues)  of  the  two  12-mode  models 
can  be  obtained  from  a  finite-element  structural  program  (e.g. ,  NASTRAN ,  SPAR, 
ASTRO,  etc.)  and  are  as  follows: 


Frequency  (Hz) 


Model  No. 

Nominal  Model 

Perturbed  Model 

1 

0.213588 

0.186316 

2 

0.264949 

0.233445 

3 

0.460071 

0.471833 

4 

0. 470687 

0.566233 

5 

0.540840 

0.612490 

6 

0.669164 

0.819555 

7 

0.741991 

0.903348 

8 

0.756823 

0.908903 

9 

1.35909 

1.42279 

10 

1.47227 

1.63987 

11 

1.63687 

1.73851 

12 

2.05391 

2.22286 
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5.1.2  The  CSDL  Example  Control  Problem 


As  described  in  the  previous  subsection ,  there  are  two  structural  models :  one 
is  the  nominal  model,  and  the  other  is  the  perturbed  model  (nominal  model  with 
parameter  variations).  It  is  assumed  that  all  12  modes  have  0.5  percent  inherent 
structural  damping.  By  definition,  modes  9,  10,  11,  and  12  are  to  be  treated 
as  unmodelled  modes,  so  that  they  cannot  be  used  in  any  controller  design 
process.  However,  they  are  to  be  used  for  the  full-order  evaluation  model  for 
all  designs . 

The  physical  model  to  be  used  in  the  controller  design  process  consists  of  the 
first  eight  modes  of  the  nominal  model.  The  evaluation  of  the  various  controller 
designs  is  then  based  on  the  full  12-mode  models  for  both  the  nominal  and  per¬ 
turbed  structures.  By  "construction",  good  performance  is  obtained  when 
modes  1,  2,  4,  and  5  are  actively  controlled  so  that  they  acquire  at  least  10  per¬ 
cent  damping.  When  that  is  achieved,  the  closed-loop  response  of  the  structure 
to  certain  initial  conditions  is  such  that  the  x-LOS  error  (LOSX)  and  y-LOS 
error  (LOSY)  (i.e. ,  the  x  and  y  components  of  the  LOS  vertex)  have  magnitudes 
less  than  0.0004  and  0.00025  units  of  length,  respectively,  after  20  seconds. 

The  rationale  for  the  statements  above  is  based  on  the  following  simulation  ex¬ 
periment:  10  percent  modal  damping  is  assigned  to  modes  1,  2,  4,  and  5,  and 
0.5  percent  to  the  remaining  modes  of  the  12-mode  nominal  model.  The  structure 
is  then  subjected  to  an  initial  condition  at  t  =  -1  sec  by  giving  a  unit  displace¬ 
ment  along  the  x-axis  to  the  front  base  vertex  of  the  tetrahedron  (see  Fig.  1). 
The  transient  response  is  then  simulated  to  t  =  20  sec,  and  it  is  observed  that 
the  LOS  requirements  are  met. 

To  make  the  problem  more  challenging,  a  "snapshot"  of  the  transient  response  is 
taken  at  t  =  0  sec,  and  values  of  the  modal  coordinates  (displacement  and  velocity) 
are  recorded.  These  values  are  then  modified  by  setting  to  zero  the  modal  coor¬ 
dinates  corresponding  to  the  unmodelled  modes  9  through  12.  The  resulting 
initial  conditions  to  be  used  at  t  =  0  sec  are  then  the  following: 


Mode  No. 


Displacement  (Q) 


Velocity  (Q) 


1 

-0.0001 

-0.003 

2 

0.0006 

0.01 

3 

0.0001 

0.03 

4 

-0.0009 

-0.02 

5 

0.0008 

0.02 

6 

-0.0001 

-0.02 

7 

-0.0002 

0.003 

8 

0.0002 

0.004 

9 

0.0 

0.0 

10 

0.0 

0.0 

12 

0.0 

0.0 

The  CSDL  example  control  problem  consists  then  of  two  parts:  (1)  design  an 
"acceptable"  controller  using  all  6  colocated  sensors  /actuators,  and  (2)  design 
an  "acceptable"  controller  that  works  without  using  colocated  sensors  and  ac¬ 
tuators.  (This  is  to  be  achieved  by  eliminating  either  a  sensor  or  an  actuator  - 
or  both  -  from  each  bipod  leg  so  that  no  colocated  sensor /actuator  pairs  remain 
in  the  system.  Many  combinations  are  obviously  possible.)  In  both  cases,  the 
controller  design  is  "acceptable"  if  the  LOS  requirements  are  essentially  met  for 
both  the  nominal  and  perturbed  12-mode  models,  with  the  controller  design  based 
only  on  the  first  8  (or  fewer)  modes.  The  primary  objective  of  the  CSDL  example 
is  to  provide  a  well-defined,  non-trivial  dynamical  plant  which  serves  as  a  common 
"test-bed"  to  illustrate,  compare,  and  evaluate  various  controller  designs.  In 
that  context,  the  exact  verification  of  the  LOS  requirements. is  only  secondary, 
and  the  given  LOS  specifications  serve  as  approximate  evaluation  guidelines. 

5.1.3  The  CSDL  Example:  Centralized  Low  Authority  Controller  Design 

This  centralized  controller  design  is  based  on  the  Low  Authority  Control  Gain 
Synthesis  procedure  previously  described  in  Section  4.1.  While  sensors  and 
actuators  are  physically  colocated,  the  feedback  in  this  design  is  not  colocated, 
and  hence  (by  definition)  is  called  centralized.  (Note:  A  decentralized  low 
authority  controller,  i.e. ,  using  colocated  feedback,  may  of  course  also  be  con¬ 
sidered.  If,  however,  sensors  are  further  assumed  to  measure  rates  only* 


♦Low  Authority  Control  Theory  is  valid  for  full  state- feedback.  Restriction  of  this 
design  to  rate- feedback  only  was  assumed  here  in  order  to  achieve  the  simplest 
controller  design  leading  to  acceptable  performance  visa  vis  the  LOS  requirements. 


*mc  h*trix  is  rx  «  non-zero  elements  follow  (gain  matrix  c  is  upper  triangular) 


Fig.  2  Centralized  Low  Authority  Control  (LAC)  Design  for  the  CSDL  Example.  Control  Gains 
and  Closed-Loop  Eigenvalues  for  Nominal  and  Perturbed  Models 


the  inherent  nature  of  the  CSDL  example  structure,  together  with  a  a  priori  fixed 
selection  of  sensors  and  actuators,  are  such  that  the  LOS  performance  require¬ 
ments  cannot  be  met.  Indeed,  synthesis  of  such  a  controller  has  shown  that  it  is 
limited  to  approximately  4  percent  (maximum)  damping  in  modes  1,  2,  4,  and  5, 
and  hence  not  acceptable.  This  result  is  physically  not  surprising:  indeed,  in¬ 
finite  gains  for  such  a  controller  would  completely  rigidify  the  bibods  and  the 
base  of  the  tetrahedron  while  leaving  the  elastic  motions  of  the  top  (LOS)  vertex 
essentially  unaffected.  As  will  be  shown  in  Section  5.1.4,  a  modification  of  the 
CSDL  control  problem,  whereby  different  sensor /actuator  types  and  locations 
can  be  used,  leads  to  an  acceptable  decentralized,  rate-feedback  low  authority 
controller  design.) 

For  the  centralized  LAC  rate-feedback  controller  design,  an  upper  triangular 
form  for  the  gain  matrix  C  was  used.  Figure  2  gives  the  values  of  the  gains 
and  shows  the  resulting  damping  ratios,  for  both  nominal  and  perturbed  models, 
obtained  by  a  full  system  eigenanalysis.  The  comparison  between  the  LAC-theory 
predicted  and  the  actual  closed-loop  damping  (for  the  nominal  model)  is  shown  in 
Fig.  3,  where  the  model  weights  are  also  displayed.  (Note:  As  in  most 
gain  synthesis  procedures,  some  trial-and-error  iterations  are  required  to  deter¬ 
mine  acceptable  gains.  In  this  case,  iterations  were  made  on  the  values  of  the 
and  of  the  desired  damping  ratios  (2  Sn  Wn  )^  =  (dA^  =  dn  aPPearing 
in  Eq.  (15)  of  Section  4.1.  A  key  feature  of  the  LAC  theory  is  that  its  closed- 
loop  performance  prediction  formula,  Eq.  (9)  of  Section  4.1,  provides  a  simple 
and  numerically  fast  method  to  evaluate  the  effect  of  these  iterations  prior  to 
final  verification  by  full  system  eigenanalysis.) 

Figures  4,  5,  and  6  show  the  results  of  LOS  transient  response  to  the  CSDL 
initial  conditions,  including  time-histories  of  LOSX,  LOSY,  physical  motion  of 
LOS,  and  the  actuator  force  levels  developed  in  the  bipod  actuators.  These 
levels  are  only  about  1-1/2  times  higher  than  those  obtained  with  the  optimal 
controller  designs  discussed  in  Section  5.1.5. 
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Centralized  Low  Authority  Control  (LAC)  Design  for  the  CSDL  Example.  Comparison 
Between  Predicted  and  Actual  Damping  Ratios  Obtained  for  the  Nominal  Model 
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Fig.  4  Centralized  LAC  Design  Evaluation  (CSDL  Example)  -  Closed-Form  Solution  of  LOS 
Transient  Responses  to  the  CSDL  Initial  Conditions.  (Dashed  Lines  Represent  LOS 
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PERTURBED  MODEL  CENTRALIZED  LAC 


Fig.  5  Centralized  LAC  Design  (CSDL  Example)  -  Projected  LOS  Motion  in  x ,  y  Plane 

Dashed  Lines  Represent  LOS  Requirements  After  20  sec.  LOS  Trajectory  Enters 
Target  Area  at  About  18  sec 


5.1.4  The  CSDL  Example:  Decentralized  Low  Authority  Controller  Design 

This  controller  design  is  not  restricted  to  the  sensor /actuator  pairs  postulated  in 
the  CSDL  example  control  problem  (Section  5.1.2).  Instead,  it  is  based  on  the 
use  of  5  sensor /actuator  pairs,  three  of  which  are  one  "leg"  of  each  bipod,  shown 
in  Fig.  7  as  ,  Fj  ,  F^  ,  and  the  remaining  two  are  inertial  proof-mass  dampers, 
designated  as  F^  ,  Fg  ,  acting  translationally  along  the  x  and  y  directions  , 
respectively. 

Figure  7  shows  the  nominal  model  open-loop  frequencies,  the  synthesized  diagonal 
gain  matrix  C,  the  LAC-predicted  and  the  actual  closed-loop  damping  ratios  ob¬ 
tained  from  full  system  eigenanalysis. 

Figure  8  shows  the  actual  damping  ratios  for  both  the  nominal  and  perturbed 
models,  and  the  LOS  transient  response  evaluations  to  the  CSDL  initial  conditions, 
also  for  both  nominal  and  perturbed  models.  Since  the  latter  were  obtained  by 
"vector  exponential"  closed-form  solution  methods,  LOSX,  and  LOSY  are  tabulated 
directly  for  the  time  interval  19  sec  <  t  <  25  sec.  As  can  be  seen,  the  LOS  require¬ 
ments  are  easily  met  at  t  =  20  sec. 
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Fig.  8  Decentralized  Controllers  -  LOS  Transient  Response  Evaluation 


5.1.5  Generalized  Collocated  Control  Design 


The  CSDL  structure  is  described  analytically  by 


where 

Aq  =  diag  (1.8,  2.7,  8.36,  8.75,  11.5,  17.7,  2.17,  22.6,  72.9,  85.6,  106,  167) 
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0.069 

-0.077 

0.271 

0.046 
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The  matrix  Qq  to  be  used  in  the  design  procedure  outlined  above  will  be  taken 
as  a  diagonal  12  x  12  matrix.  j  j 

In  the  following  iterations  the  diagonal  of  the  matrix  Q  will  be  selected.  Using  \ 

U  k 

the  design  procedure  outlined  in  4.2,  Q  and  K  will  be  found  and,  also,  the  f 

resulting  closed-loop  eigenvalue  system.  > 

High  Authority  Control  Design 


In  the  high  authority  controller  selection  the  following  weighting  matrices  are 
selected 

A(jW)  =  A 

B(jtf)  =  Diag  [b^W2  +  W2),  b2(W2  +  W2).  ... 

bgOV2  +  W2)] 

A(jW)  does  not  depend  on  the  frequency.  B(jW)  can  be  factered  as 
B (jW)  =  Diag  (-jW  +  Wj,  -jW  +  Wg  ....  -jW  +  Wg)  x 

Diag  (bj,  b2  ...,  bg)  x  Diag  (jW  +  Wj,  jW  +  W2  ...,  jW  +  Wg) 

If  we  define 


then 


Diag  (b^  b2 


...  bg)Z 


U*  B  (jW)u 


The  control  design  problem  is  then  described  by  the  following  set  of  equations. 


x  =  Fx  +  Gy  .  ■  16  Equations  (8  modes) 

u.  =  Wj  U.  +  Z.  —  -  16  Equations 

Cost  Functional 


T  T 

x  Ax  +  Z  RZ 

R  Diag  Cb|i  b2  •  •  «bg) 


Control  Law 


Z  =  CjX  +  C2u 


Filter 

X  =  Fx  +  Gu  +  K  (y  -  Hx) 


The  following  weighting  functions  were  used  with  six  actuators  and  six  sensors 
(collocated  system). 


State 

Weighting 


Diag  (0.5,  0.5,  0.5,  0.5,  0.1,  0.1,  0.5,  0.5,  1.0,  1.0,  6*  0.001) 

Modes  in  which  10%  damping 
Ratio  is  desired 


Control 

Weighting 


Process 

Noise 


Diag  (0.01,  0.01,  0.01,  0.01,  0.01,  0.01) 

All  actuators  are 
equally  important 

Diag  (10,  10,  4,  4,  4,  0.4,  0.4,  0.0) 
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Meas 

Noise 


=  Unity  Matrix 


Shaping 

Frequencies 


(3.1,  3.2,  3.3,  3.4,  3.5,  3.6)  Rad/Sec 


Closed  loop  damping  ratio  for  the  nominal  control  design  and  evaluation  modes  are 
shown  in  Table  1  .  Note  the  excellent  behavior  of  high  frequency  modes  not 
considered  in  control  design  model.  Table  2  shows  results  on  the  perturbed 
model . 


The  control  design  procedure  was  repeated  with  the  first  three  actuators  and  the 
last  three  sensors  (noncolocated  system) .  The  following  parameters  were  used  in 
the  design  stage. 


State 

Weighting 


=  (4*1.0,  2*0.1,  4*1.0,  6*0.001) 


Modes  in  which  10%  damping 
Ratio  is  desired 


Control 

Weighting 


(3*0.01)  First  three  actuators  are  used  for  control 


Process 

Noise 


=  (30,  30,  8,  8,  8,  0.8,  0.8,  0.0) 


Meas 

Noise 


Unity  Matrix 


Shaping 

Frequencies 


(3.1,  3.2,  3.3)  Rad/Sec 


Table  3  compares  the  control  design  and  evaluation  models.  Note  again  that  the 
modes  not  considered  in  control  design  are  not  significantly  effected.  When  the 
controller  is  used  on  the  perturbed  modes ,  satisfactory  results  are  predicted 
(see  Table  4  ). 


1 
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Table  1 

NOMINAL  DAMPING  RATIOS 


.00491 


COMPARISON  OF  NOMINAL  AND  PERTURBED  CLOSED-LOOP  SYSTEMS  -  COLOCATED  CONTROL 


NOMINAL  DAMPING  RATIOS  NONCOLOCATED  DESIGN 


COMPARISON  OF  NOMINAL  AND  PERTURBED  CLOSED-LOOP  SYSTEMS  -  NONCOLOCATED  CONTROL 


.00495  '  I.IZ  .00399 


5.1.6  Modal  Cost  Analysis  (MCA)  Results  for  the  CSDL  Example 


Here  we  ask  for  additional  insight  into  the  MCA  (see  Ref  1,  Sec  4.2.4)  methods 
by  solving  a  standard  optimal  control  problem  with  two  different  reduced  models: 
Model  Sg  is  the  16th  order  model  obtained  by  keeping  the  8  lowest  frequency 
modes;  Model  is  the  16th  order  model  obtained  by  keeping  the  8  modes  with 
the  largest  open  loop  modal  cost.  The  nominal  evaluation  model  is  of  order  24 
and  has  the  form 

x  =  Ax  +  Cu  +  w,  E|ww^|  =  [x(0)xT  (0)] 
y  =  Px  line  of  sight  errors  at  tip  (1) 

z  =  Mx  displacement  and  rate  measurements  in  legs 


The  quadratic  cost  for  this  problem  is  defined  by 


yTQy  =  V*  +  njj  +  100  +r?2) 


LOS  Error  Adds  Damping 
to  LOS 
Regulation 


The  open  loop  modal  costs  V  are  defined  by 

ci 

'  T  H 

V  =  lim  y AQy  =  £  V 

t  — 80  i=l  l 


°i  ’  «a4  ’ 


where  6a  are  modal  amplitudes  and  u>a  are  modal  frequencies. 


Figure  1  shows  the  resulting  modal  costs  for  each  mode.  On  the  vertical  scale 
the  frequencies  of  the  modes  are  shown.  Thus,  the  eight  modes  contained  in  the 
design  model  Sg  lie  below  the  dashed  horizontal  line.  The  eight  modes  contained 
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in  the  design  model  lie  to  the  right  of  the  vertical  dashed  line.  The  modes  in 
common  to  both  models  are  1,  2,  4,  5,  7.  8.  Modal  cost  analysis  indicates  modes 
3  and  6  are  least  significant  among  all  modes,  while  modes  3  and  6  rank  among 
the  8  most  significant  in  the  "dominant  mode"  model  Sg.  Table  1  shows  that  the 
damping  added  to  each  optimal  design  has  approximately  the  same  magnitude  but 
is  distributed  differently  among  the  modes.  The  MCA  design  selected  modes  9 
and  11  for  additional  damping,  while  the  dominant  mode  design  (with  model  S^) 
attempted  to  add  damping  to  modes  3  and  6.  The  control  effort  was  approxi¬ 
mately  the  same  for  both  designs  but  the  line  of  sight  regulation  was  written 
better  by  37  percent  with  the  MCA  design. 

The  results  using  MCA  should  be  even  better  than  the  above  results  by  using 
closed-loop  modal  cost  analysis.  Without  such  additional  results,  however,  we 
can  at  least  say  that  open-loop  MCA  has  potential  in  the  reduction  of  large  models 
(obtained  from  finite  element  programs)  to  lower  order  models,  before  any 
closed-loop  analysis  is  done. 

Table  1 

CSDL  EXAMPLE 


Dominant  Mode  Model  _ MCA 


Use  Optimal  Control  for  Sq 

Use  Optimal  Control  for  S. 

to  Drive  S^: 

to  Drive  S^: 

Mode 

Damping 

Mode 

Damping 

1 

0.29 

1 

0.29 

2 

0.34 

2 

0.35 

3 

0.03 

3 

0.005 

4 

0.08. 

4 

0.08 

5 

0. 11 

5 

0.11 

6 

0.005 

6 

0.005 

7 

0.08 

7 

0.08 

8 

0.05 

8 

0.05 

9 

0.005 

9 

0.05 

10 

0.005 

10 

0.005 

11 

0.005 

11 

0.04 

12 

0.005 

12 

0.005 

V  =  £  [yTQy 

T 

+  u  u 

LOS  Regulation 

:  8.84 

6.46 

Control  Effort: 

6.17 

6.34 
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5.2  LBET  SYSTEM 


5.2.1  System  Description 

The  LBET  strawman  system  provides  an  ideal  test  platform  for  active  vibration 
suppression.  Although  on-board  disturbances  (state  noise)  do  not  excessively 
perturb  the  attitude  control  system  nor  do  structural  modes  fall  within  the  atti¬ 
tude  control  bandwidth,  suppression  of  vibration  due  to  on-board  state  noise 
sources  is  essential  to  meet  mission  requirements.  The  system  is  diagramed  in 
Fig.  1  along  with  candidate  locations  for  proof-mass  actuator  systems  used  to 
increase  closed-loop  damping  of  structural  modes.  The  selection  of  structural 
materials  and  design  rationale  are  fully  discussed  in  the  LMSC  ADOPT  TASK  E 
final  report  and  will  not  be  addressed  here.  Some  typical  low-frequency  modes 
from  the  finite-element  model  are  illustrated  in  Fig.  2. 

For  worst  case  values ,  system  pointing  and  wave  front  errors  for  the  unaug¬ 
mented  system  are  approximately  two  orders-of-magnitude  greater  than  specifi¬ 
cation  requirements.  Furthermore,  these  values  assume  1  percent  structural 
damping  in  the  first  20  modes.  This  assumption  is  generally  untenable  for  Gr-Ep 
structural  designs  where  experimental  evidence  indicates  that  0.3  percent  is 
more  realistic.* 

5.2.2  SAS  Performance 

The  principal  modal  interaction  with  pointing  performance  is  due  to  the  tracker 
attachment  to  the  spacecraft.  Several  locations  for  proof-mass  damping  actuators 
were  tried,  as  shown  in  the  system  sketch,  with  the  most  successful  being  loca¬ 
tion  1  at  the  tracker.  Participation  of  other  modes  in  performance  was  down  by 
at  least  a  factor  of  10  and  increased  damping  in  these  modes  did  not  significantly 


♦Measurement  of  modal  structural  damping  is  a  controversial  issue  which  will  be 
addressed  in  the  LMSC  ACOSS  phase  1A  study. 
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Proof-Mass  Damper  Locations  and  Control  Axes 


SPIDER  &  SECONDARY  MIRROR 


LBET  Baseline  Dynamic  Model  and  Dominant  Modes 


enhance  performance.  The  open-  and  closed-loop  damping  performance  obtained 
with  a  single  collocated  (decentralized)  system  is  shown  by  the  computer  gen¬ 
erated  design  in  Table  1.  Open-loop  damping  was  assumed  to  be  0.5  percent  and 
the  significant  modes  are  enclosed  by  the  boxes.  The  synthesis  procedure  was 
asked  to  produce  10  percent  in  the  two  critical  modes  and  9  percent  and  7.7  per¬ 
cent  were  produced  as  shown. 

The  principal  advantages  which  accrue  from  this  mechanization  are  relaxed  sys¬ 
tem  requirements  and  decreased  sensitivity  to  knowledge  of  plant  dynamics  and 
disturbance  sources. 

Specifically,  requirements  to  disable  on-board  noise  sources,  and  to  maintain  high 
bandwidth  figure  control  and  steering  mirrors  may  be  eased  considerably. 
Assumptions  about  passive  structural  damping  no  longer  become  so  significant, 
tolerance  to  non-isolatable  disturbances  is  high,  and  transient  settling  due  to 
large,  on-board  disturbances  becomes  acceptable. 
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5.3  MM-WAVE SYSTEM 


5.3.1  System  Description 

The  MM-wave  strawman  system  is  a  30  m  by  60  m  curved  platform  originally  pro¬ 
posed  by  CSDL.  The  configuration  and  moment  of  inertia  properties  are  shown 
in  Fig.  1.  For  convenience  in  modeling,  the  reflector  back  structure  was  modi¬ 
fied  to  include  the  NASA-LRC  octetruss  configuration  using  Gr-Ep  columns  man¬ 
ufactured  at  LMSC.  The  new  properties  on  the  right  of  the  figure  represent  the 
octetruss  values.  As  shown  in  the  accompanying  Table  1,  the  octetruss  reduces 
both  back  structure  and  feed  support  weight. 

5.3.2  SAS  Performance 

To  show  SAS  feasibility,  the  LASC  synthesis  procedure  (see  Section  4)  was  used 
to  select  feedback  gains  for  actuators  located  by  visually  observing  modal  de¬ 
flections  to  assess  controllability.  A  two-axis  gyro  damper  and  a  single  proof- 
mass  damper  were  installed  on  the  reflector  and  a  three-axis  gyro  damper  was 
used  to  control  the  feed  tower.  These  actuator  locations  are  shown  in  Fig.  2 
along  with  frequency  and  damping  characteristics  of  the  closed-loop  system. 

The  controls  are  all  of  the  colocated  decentralized  low-authority  type  without 
consideration  for  bandwidth  constraints.  Moderate  amounts  of  damping  are  in¬ 
troduced  up  to  20  Hz  with  the  principal  effects  occurring  in  the  low  frequency 
modes  below  1.6  Hz.  In  Figs.  3a  and  3b  the  feed  decenter  transfer  functions  are 
compared  for  open-  and  closed-loop  and  it  may  be  observed  that  peak  responses 
are  reduced  by  a  factor  of  approximately  6  to  1. 


i 


to 

to 

§z 

aiuiC 

ZOq 


to 

to 


-2s 

CM 


Z 

o 


CN I  CM  CM 

E  E  E 


,f  E 

op 

_  CM  K 
O  o  CO 
•  CO  • 
O  I  Tt 


O 

85 

t— 

z 


xxx 

•—  CO  CM 
CM  00  O' 
CN  U-)  O' 

•  •  > 


tZ 

'  to 
:  lu 
>  Q 


2 

o 


u 


_s> 


w 

3 


e  e 

•  •  0 

O  O  Ma¬ 
lt  K  ff 


3 

o 

ca 

< 

1/5 

< 

h“ 

Of 


CM  CM  CM 

SEE 

-2>-S5^> 

XXX 
K  o  00 
CO  lo  CM 
—  M->  K 

9  0  0 

CM  IT)  O 

nun 


X>N 

X>N 


ujQ 

Q  MJ 

y  <x 

03 

t250 

U=  iZ 
ixn—  -r 

u!<  § 

Sxu 

<{£>- 
£  s 

<u.  t 
QOz 

kj3 

£3<J 

^  t/j 

<  oc  < 


•  • 
UJ 

h- 

O 

Z 


bio 

E 


&- 


mm-Wave  Antenna  Characteristics 


30  x  60  mm  mm-WAVE  ANTENNA  MASS  ALLOCATIONS 
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Fig.  2  mm-Wave  Antenna  (12  Modes)  -  LASC  Actuators  and  Performance 
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5.4  HALO /WALRUS  SYSTEM 


5.4.1  System  Description 

The  HALO  (High  Altitude  Large  Optics)  System  examined  for  stability  augmenta¬ 
tion  was  the  HUGHES  WALRUS  (Wide-Angle  Large  Reflecting  Unobscured  System). 
The  optical  characteristics  and  required  element  tolerances  are  described  in  vari¬ 
ous  HALO  documents  contained  in  the  Riverside  Research  Library  and  in  the 
LMSC  proposal  for  An  Actively  Controlled  Structure  Program.  Vol  1.  17  April  78. 
LMSC-L034034(S ) .  The  discussions  here  will  be  restricted  to  the  feasibility  and 
potential  performance  of  active  stability  augmentation  using  only  low-authority 
control.  The  basic  configuration,  shown  in  Fig.  1,  uses  an  LMSC  designed 
structure  in  lieu  of  the  cable  support  system  used  on  the  original  Hughes  Config¬ 
uration.  The  attitude  control  system,  eryo-coolers ,  etc.,  are  housed  in  the 
so-called  "dirty  box"  which  is  not  actively  isolated  but  rather  attached  directly 
to  the  optical  structure  as  shown.  The  mirrors  are  supported  by  pyramid 
structures  which  were  weight  optimized  as  shown  in  Fig.  2.  The  overall  mass 
properties  are  reproduced  in  Table  1. 

5.4.2  SAS  Performance 

A  coarse  finite-element  model,  illustrated  in  Fig.  3,  was  used  to  examine  basic 
feasibility.  Decentralized  low-authority  control  (with  no  bandwidth  constraints) 
of  the  fundamental  low-frequency  structural  modes  was  attempted  while  higher 
frequency  "mirror  modes"  were  left  uncontrolled.  For  a  practical  design,  a  more 
refined  model  is  necessary  and  actuator /sensor  bandwidth  limitations  would  have 
to  be  included  to  insure  that  feedback  compensation  and  passive  low-authority 
actuator  damping  were  appropriately  selected.  The  actuator  locations,  charac¬ 
teristics  and  closed -loop  damping  performance  are  given  in  Fig.  4.  The  modes 
between  1  Hz  and  16  Hz  all  exhibit  damping  greater  than  10%  except  for  a  few  un¬ 
controlled  modes;  this  is  accomplished  by  8  proof-mass  systems  and  1  (rotational) 
gyro  damper.  Individual  mirror  modes,  all  above  24  Hz.  exhibited  less  than  1% 
damping,  as  expected,  since  these  modes  were  uncontrollable  by  the  actuators 


Fig.  2  Halo/Walrus  Structural  Model  (General  View) 
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selected.  Reversal  of  high  q  structural  resonances  in  this  case  eases  significantly 
any  requirements  to  isolate  vibration  sources  on-board  the  spacecraft  and  permits 
greater  freedom  in  synthesis  of  the  attitude  control  system. 
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Fig.  4  Halo/Walrus  (30  Modes)  -  LASC  Actuators  and  Performance 


5.5  TILTING  TELESCOPE  SYSTEM  (ADOPT-12) 

The  structural  description  of  the  system  considered  in  this  section  is  given  in 
Fig.  1.  The  modes  of  this  system  are  shown  in  Fig.  2.  The  model  used  in  the 
control  design  has  26  states,  13  controls,  and  26  measurements;  these  are  de¬ 
tailed  as  follows: 

26  State  Model 

•  2  rigid  body  modes 

•  11  structural  modes 
13  Controls 

•  12  surface  controllers 

•  3  back  body  torquer 
26  Measurements 

•  12  surface  position  96  independent) 

•  12  surface  velocity  (6  independent) 

•  1  relative  angle 

•  1  average  orientation  angle 

1  Disturbance  Source 

•  Lateral  force  on  back  body 

The  control  design  for  this  system  was  concerned  with  the  following  problems: 

1.  Reduced-order  model  derivation 

2.  Comparison  of  LQG  and  frequency  shaped  cost  functional  design  method 

3.  Filter  design 

In  the  derivation  of  the  reduced-order  model  desired  states  are  retained.  In 
addtion  the  reduced  order  model  has  the  same  poles  as  the  high  order  model.  Due 
to  symmetry,  duplicate  sensors /actuators  can  be  eliminated  in  the  control  design 
based  on  the  reduced  model.  Details  of  the  reduction  process  are  shown  in  Fig.  3. 
The  effect  of  the  states  retained  in  the  reduced  model  on  the  zero  location  is  shown 
in  Fig.  4. 
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Fig.  1  Tilting  Telescope  Final  Structural  Design 
Primary  Mirror  Support  Structure 


N-SOOY  NONLINEAR  DYNAMIC  SIMULATION 


Selected  System  Modes 


GLOBAL  MODEL  REDUCED  MODEL 


DUPLICATE  SENSORS/ACTUATORS  ARE  ELIMINATED  IN  CONTROL  DESIGN 

Fig.  3  Tilting  Telescope  System  Definition 


Selection  on  Reduced  System  Zero  Location 


RELATIVE  CONTROL  WEIGHTING 

Fig.  5  Damping  Ratio  of  Controlled  Modes  (LQG  Design) 


The  next  stage  in  the  design  was  the  comparison  between  frequency  shaped  cost 
functionals  and  the  standard  LQG  design  method.  The  basis  of  comparison  was 
taken  as  the  damping  ratios  of  various  modes.  The  results  are  given  in  Figs.  5 
through  8. 

An  important  issue  in  the  application  of  frequency  shaped  cost  functionals  is  the 
choice  of  these  functionals.  The  choice  must  be  adapted  to  the  particular  situation 
at  hand. 

As  mentioned  before,  a  problem  associated  with  large  space  structures  is  that  of 
truncated  modes.  An  approach  toward  dealing  with  this  problem  is  offered  through 
the  application  of  frequency  shaped  cost  functionals.  The  weighting  on  the  control 
is  increased  at  high  frequency  to  account  for  the  model  invalidity  in  that  regim, 
see  Fig.  9,  leaving  the  state  weighting  constant  the  effect  of  control  weighting  can 
be  observed.  In  Fig.  10  the  change  in  damping  raio  as  a  function  of  control  weight¬ 
ing  is  depicted. 

The  last  topic  considered  for  this  structure  was  that  of  filter  design.  Filter  must 
be  designed  for  good  transient  as  well  as  rms  response.  If  only  rms  response  is 
considered  poor  measurement  responses  may  result  (fast  poles). 

Additional  requirements  are: 

•  Very  little  change  in  the  filter  closed-loop  poles  based  on  original  power 
spectral  density  of  the  noise 

•  Noise  matrix  scaled  to  improve  system  behavior 

•  Requires  reduction  of  compensator  for  better  design 

Figures  11  and  12  illustrate  these  points. 

Figure  13  shows  the  rms  response  in  steady  state.  As  the  control  weighting  de¬ 
creases  (a)  the  damping  ratio  increases,  (b)  the  rms  furface  deflection  improves 
(i.e.,  becomes  smaller),  (c)  the  difference  between  low  order  and  high  order 
model  increases,  and  (d)  more  actuator  power  is  required.  Note  that  this  figure 
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does  not  include  effects  of  actuator  sensor  noise.  Those  two  effects  could  sig¬ 
nificantly  degrade  performance.  Control  force  requirements  are  shown  in  Figs. 
14  and  15.  It  can  be  observed  that  with  decreased  control  weighting  associated 
with  increased  damping  ratio,  the  actuator  requirements  are  more  severe. 


I 


RELATIVE  CONTROL  WEIGHTING 

Fig.  6  Damping  Ratio  of  Controlled  Modes  (Frequency  Shaped  Cost  Functional) 


Fig.  7  Damping  Ratio  of  Uncontrolled  Modes  (LQG  Design) 


Fig.  8  Damping  Ra 


Fig.  9  Control  Design  Weighting  Functions 


Fig.  10  Damping  Ratio  as  a  Function  of  Relative  Control  Weighting 


Fig.  11  Filter  Closed  Loop  Poles  with  Original  Noise  PSD 


Fig.  12  Filte 


Section  6 

EXPERIMENTAL  RESULTS 

The  experimental  results  presented  in  this  section  were  created  under  LMSC 
IRAD  funding  in  general  support  of  this  area  of  research  and  are  included  to 
clarify  the  applicability  of  the  theoretical  work  described  in  earlier  sections. 

6.1  PROOF-MASS  DAMPED  BEAM  (MINIBEAM)* 

This  experimental  brass-board  was  developed  under  the  Lockheed  Independent 
Research  Program  to  illustrate  the  principle  of  an  inertial  translational  proof- 
mass  damper  for  structural  vibrations.  While  the  intended  application  is  generally 
for  low  authority  control  (i.e. ,  10  percent  to  20  percent  damping) ,  the  scaling  of 
this  particular  experiment  (which  uses  off-the-shelf  hardware  components  de¬ 
signed  for  other  purposes)  was  such  that  over  50  percent  damping  was  obtained. 

6.1.1  Experimental  Set-Up 

A  40-in.  long  magnesium  bar  is  clamped  on  a  laboratory  table  (Fig.  1)  allowing 
essentially  only  horizontal  cantilevered  structural  oscillations  with  the  first  mode 
around  4.5  Hz.  A  2-lb  Ling  shaker  unit  is  mounted  at  the  free  end  of  the  beam 
in  a  pivoted  plexiglass  cradle  supporting  the  shaker  caging,  while  the  shaker 
stem  is  solidly  attached  to  the  beam.  The  cradle  thus  allows  free,  small  horizon¬ 
tal  motions  of  the  shaker,  but  constrains  its  vertical  motions  which  would  other- 
wise  produce  excessive  bending  moments  on  its  stem.  As  designed,  this  system 
is  intended  to  represent  an  idealized  one-dimensional  proof-mass  damper  in  the 
horizontal  plane. 

To  implement  a  rate  sensing  device,  an  optical  method  is  used.  A  laser  beam 
(emitted  from  a  helium-neon  laser  source  at  the  root  of  the  magnesium  bar)  is 


*This  experiment  is  shown  and  described  on  the  DARPA/ACOSS  videotape.  Since 
its  purpose  is  only  to  demonstrate  a  physical  principle,  no  quantitative  results 
will  be  given. 


6-1 


Fig.  1  Proof-Mass  Damped  Beam  (Minibeam) 


1 


folded  (reflected)  several  times  by  a  system  of  mirrors,  as  shown  in  Fig.  2,  in 
order  to  produce  an  oscillating  spot  corresponding  to  the  beam's  tip  deflections. 
The  spot  is  picked  up  by  a  linear  detector,  and  via  suitable  electronic  compensa¬ 
tion,  produces  a  tip-deflection  rate  signal  used  as  feedback  to  drive  the  shaker 
with  proper  phase  and  amplitude. 

6.1.2  Analytical  Model  and  Actuator  Dynamics 

Let  m  be  the  mass  of  the  shaker  unit  ("proof-mass  damper")  and  let  x,  X 
respectively  denote  the  inertial  displacements  of  the  shaker  and  of  the  beam  tip, 
as  shown  in  Fig.  2. 


Fig.  2  Notation  for  Proof-Mass  Dynamics 

Let  f  denote  the  electromagnetic  force  applied  to  the  shaker,  fg  the  actual 
force  applied  to  the  beam,  and  k  the  (weak)  spring  constant  inherent  to  the 

caging  of  the  shaker  stem.  Since  f  is  proportional  to  the  beam-tip  rate,  the 

equations  of  motion  for  the  system  are 

mx  =  f  =  k  (X  -  x) 

e 

fe  =  ~fe  +  k  (x  -  X) 

f  =  DX 

e 
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In  Laplace  transformed  variables  we  have  then 


(ms2 


+  k)  x 


(Ds  +  k)  X 


and  hence 


x  -  X  =  sX 


D  -  ms 
ms2  +  k 


(2) 


Substituting  (2)  into  the  middle  equation  of  (1)  yields,  after  some  simple  algebra, 

two  equivalent  forms  for  the  force  f  applied  to  the  beam  tip: 

3 


f 


a 


-DsX 


k  (1  -  ms/D) 
k  +  ms2 


(3) 


f 


a 


-Ds2X 


ms/k  -  m/D 
1  +  ms2/k 


(4) 


From  these  two  expressions,  it  can  be  seen  that: 

1)  for  S  »  (k/m)1/2,  f  -  -DsX 

3 

2)  for  S  «  (k/m)1/2  f  ~ms2X 

’  a 

The  above  shows  the  characteristic  property  of  proof-mass  dampers  (such  as  the 
Ling  shaker)  having  a  weak  spring  caging  restraint.  For  sufficiently  large  fre¬ 
quencies,  their  behavior  approximates  that  of  an  idealized  inertial  damper, 
whereas  for  sufficiently  low  frequencies,  their  effect  is  essentially  to  mimic  the 
presence  of  an  added  mass  m  attached  to  the  structure  (beam). 


Finally,  the  colocated  loop-closure  is  produced  by  using  the  optical  position  sen¬ 
sor  with  a  classical  rate  circuit  compensation  analogous  to  the  one  used  for  the 
maxibeam  experiment  described  in  the  next  section. 
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6.2  GYRODAMPED  BEAM  (MAXIBEAM) * 


This  experimental  brass-board  was  developed  under  the  Lockheed  Independent 
Program  to  illustrate  the  principle  of  an  inertial  rotational  gyro-damper  for 
structural  vibrations.  The  major  hardware  component  for  this  experiment  is  a 
pair  of  identical  single-gimbal  control  moment  gyros  (CMGs).  CMGs  are  not  off- 
the-shelf  pieces  of  hardware  and  are  generally  unique,  "few-of-a-kind"  develop¬ 
mental  items,  as  shown  in  Table  1,  which  lists  CMGs  built  by  the  Bendix 
Corporation.  The  two  CMGs  used  in  this  experimental  brass-board  are  model 
no.  MA-5- 100-1,  and  are  part  of  a  set  of  5  identical  units  built  for  LMSC  in  1973. 
(The  shared  developmental  cost  at  the  time  was  $60,000  per  unit.) 

While  gyrodampers  would  normally  be  sized  and  fabricated  for  specific  spacecraft 
structural  stabilization  systems,  this  experimental  brass-board  had  to  be  de¬ 
signed  "in  reverse".  That  is,  given  the  two  available  Bendix  CMGs,  weighing 
approximately  40  lb  each,  a  massive  flexible  support  structure  had  to  be  de¬ 
signed  to  match  the  performance  characteristics  of  these  CMGs.  Approximately 
10  to  20  percent  damping  was  obtained  with  this  experiment. 

6.2.1  Experimental  Set-Up 

A  massive  aluminum  I-beam,  weighing  almost  400  lb,  with  16  in.  web  and  6  in. 
flanges,  25  feet  long,  is  cantilevered  on  a  special  1000  lb  support  structure  so 
that  its  "soft-§xis"  deflections  are  essentially  constrained  to  occur  in  the  hori¬ 
zontal  plane,  with  the  first  mode  around  1  Hz.  At  the  free  end  of  the  beam,  the 
two  MA-5-100-1  Bendix  CMGs  are  symmetrically  mounted  (Fig.  3)  in  a  scissoring 
V-CMG  pair,  with  gimbal  axes  parallel  to  the  long  axis  of  the  beam.  In  the  caged 
position  (i.e.  ,  zero  gimbal  angles),  the  CMG  angular  momenta  are  equal  and  op¬ 
posite  in  the  horizontal  plane,  and  hence  cancel  (zero  stored  momentum).  When 
the  beam  oscillates  horizontally,  the  V-CMG  pair  scissors  so  that  the  resultant 


♦This  experiment  is  shown  and  described  on  the  DARPA/ACOSS  videotape.  Since 
its  purpose  is  only  to  demonstrate  a  physical  principle,  no  quantitative  results 
will  be  given. 


Table  1 

BENDIX  CONTROL  MOMENT  GYROS 


PHOTO  DIODE 
LINEAR  DETECTOR 


I 


momentum  (and  hence  the  output  torque)  always  lies  along  the  vertical  axis  nor¬ 
mal  to  the  horizontal  plane  of  beam  vibrations.  A  summary  of  the  maxibeam  ex¬ 
perimental  data  is  given  in  Table  2. 


To  implement  a  rate-sensing  device,  an  optical  system  is  used.  A  laser  beam, 
emitted  at  the  root  of  the  beam,  is  folded  (reflected)  once  by  a  mirror  colocated 
with  the  CMGs,  and  produces  an  oscillating  spot  picked  up  by  a  linear  detector 
colocated  with  the  laser  source.  The  signal  is  then  processed  electronically  to 
produce  beam-tip  rotational  rates,  and  with  suitable  compensation,  is  used  as 
(rotational)  rate-feedback  to  drive  the  CMG  gimbal  torquer  motors  with  proper 
phase  and  amplitude.  (The  CMG  rotors  have  a  constant  8000  rpm  spin,  and  at 
20  rad /sec  gimbal  rate,  the  output  torque  of  each  CMG  is  about  100  ft -lb;  see 
Table  1). 


Table  2 

MAXIBEAM  EXPERIMENT  DATA 


Beam 

M  =  154.6  kg  (effective  cantilevered  mass) 
L  =  6.96  m 

I  =  4.375  10“6m4 

E  =  7.2  1010  Pa 


Actuators 
Total  mass 

Max  gimbal  sec 
Max  output  torque 


34.  5  kg 
20  rad /sec 
1650  rad/sec^ 
120  nm 


System  Modes  (cantilever  beam  with  tip  mass) 


Mode  No. 


1 

2 

3 

4 


Frequency 


1.000 

7.064 

20.86 

42.22 


Mode  Shape  at  Tip 
(rotation)  rad  kg~l/2 

0.0244 
-0.0884 
0. 1425 
-0.1966 


Control  gain  for  10%  damping  in  first  mode:  Da  -  2105  nm/rad-s 
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6.2.2  Analytical  Model  and  Actuator  Dynamics 

Consider  a  small  rigid  portion  B  of  the  oscillating  maxibeam,  with  an  orthonormal 
frame  je^,  e2>  e^  j  attached  to  it  as  shown  in  Fig.  4.  For  simplicity  of  illus¬ 
trating  the  basic  principles  involved,  we  assume  that  B  can  only  undergo  planar 
oscillatory  rotations  about  e^,  and  that  "e^  maintains  an  inertially  fixed  orienta¬ 
tion.  Then  the  inertial  angular  velocity  fig  of  B  (i.e. ,  of  the  frame  j  e^,  e2>  e^Jj 
is 

=  0^  (1) 

where  0  is  the  angle  of  rotation  of  B  about  e^.  Consider  first  gyro  no.  1,  with 
constant  speed  rotor  having  angular  momentum  h^,  gimbaled  on  e^,  the  gimbal 
angle  o  being  zero  when  h^  is  parallel  to  -e2.  The  total  angular  momentum  H 
the  gyro,  expressed  in  the  moving  frame  ei>  e2*  e3  *s  tJien  gfiven  by 

H  =  h  (sin  cr  e  j  -  cos  a  e,,)  +  J  b  e^  (2) 

where  h  =  Ih^l  =  const,  and  J  is  the  gyro's  inertia  about  the  gimbal  axis 
(direction)  e^. 

Let  T  be  the  torque  applied  externally  to  the  gyro  (e.g.,  by  the  gimbal 
torquer  motor).  The  equation  of  (gimbal)  motion  of  the  gyro  is  then 

5T  =  Tg,  (3) 


where  d(  )/dt  denotes  time  differentiation  w.r.t  inertial  space.  If  (’)  denotes 
time  differentiation  in  the  moving  frame  Je^,  e2>  e^  J  ,  and  if  we  assume  small 
motion  linearized  dynamics,  then  from  (1),  (2),  and  (3)  we  obtain: 


-xr  h  (cre.  -  e9)  +  Jae, 


=  h  (oej  +  x  e1  -!2g  x  e2) 

+  Jcre^  +  Jo-fig  x  e^ 

h  (oej  -  0e3)  +  Jcfeg  =  Tg 
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Fig.  4  Notation  for  CMG  Dynamics 

Assuming  now  that  Tg  =  Tge3  is  an  applied  gimbal  torque,  projection  of  the 
last  vector  equation  (4)  along  results  in 

Jcr  =  h0  +  Tg  (5) 

which  represents  the  actuator  (linearized)  dynamics  for  each  gyro.  If,  further¬ 
more,  T  =  -  k  cr  -  k  <j,  where  k  >  0,  k  >  0  are  spring  and  dashpot  con- 
&  c  g  c  g  _ 

stants,  Eq.  (5)  represents  a  passive  dissipative  system  driven  by  the  input  fl(t). 
This  passive  gyro  system  will  therefore  absorb  beam-rotational  energy  by  oppos¬ 
ing  the  0  motion.  However,  the  transmission  of  this  energy  is  determined  by 
the  magnitude  of  h/J,  which,  in  most  practical  cases  of  oscillating  beams  or 
structures,  is  too  small  for  effective  damping. 

When  two  identical  gyros  are  used  such  that  zero  gimbal  angles  correspond  to 
gyro  momenta  alignments  of  h^  parallel  to  -e9,  Jig  parallel  to  &2<  then  the 
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counter-rotation  of  the  rotors  cancels  the  internally  stored  momentum  for  cr  =  0. 
In  addition ,  because  of  symmetry ,  the  -induced  gyro  gimbal  motions  result  in 
the  classical  "V-scissor"  pairing  of  the  momenta  h^  and  h2,  i.e.,  the  resultant 
total  gyro  momentum  (w.r.t  the  moving  frame  [  e^  e2>  e3j)  always  remains 
aligned  with  the  ±e1  beam  axis.  This  allows  utilization  of  increased  gyro  gimbal 
travel.  This  travel  is  limited  in  practice  due  to  the  nonlinear  relation  between 
gyro  momentum  and  gimbal  angle,  i.e.  , 

|h1  +  h2  *  2h  sincr  (6) 

where  h  =  |  h  1 1  =  |  h2  I  •  Finally ,  the  output  torque  Tg  applied  to  the  beam 
by  the  gyros  is : 

s  -2hoe1 

6.2.3  Control  Equations  for  Active  Augmentation  of  Passive  Gyrodamper 

The  inertial  absorption  of  the  9 -rotational  energy  by  the  passive  gyrodamper  is 
limited  in  practice  by  the  magnitude  of  h/J.  To  increase  this  energy  dissipa¬ 
tion  ,  an  angular  rate  sensor  is  added  to  the  system  to  estimate  the  velocity 
which  is  fed  back  to  the  gyro's  gimbal  motor  with  the  proper  gain  so  that  the 
colocated  sensor /actuator  pair  will  mimic  a  purely  passive  device.  As  is  well 
known,  this  passive  nature  alone  guarantees  structural  stability,  i.e.,  this 
(idealized)  device  cannot  excite  the  structure,  but  the  structure  can  excite  it. 
Such  an  actively  augmented  passive  damper  can  be  caricatured  as  a  passive  de¬ 
vice  with  a  "hearing  aid"  to  magnify  the  0  signal,  and  in  this  sense,  is  quite 
distinct  from  a  purely  active  device  which  does  not  necessarily  require  either  a 
colocated  sensor  or  a  rate  feedback. 

The  maxibeam  brass-board  is  an  implementation  of  this  concept  for  a  system  of 
two  identical,  coaxially  gimbaled  gyros  (so-called  V- gyros)  as  shown  in  Fig.  1. 


Equation  (5)  remains  valid  for  each  gyro,  but  the  gimbal  torque  is  no 
longer  realized  by  passive  springs  and  dashpots,  but  generated  electronically 
by  gimbal  motors  and  will  include  a  feedback  team  proportional  to  the  sensor 
output  0.  The  control  equations  of  each  gyrodamper  are  now  described  by: 


Control  Equation: 


•  -  _D  a 

°c  "  2h  9 


(8) 


CMG  Rate-Servo 
and  Caging  Loops: 


(9) 


rate-loop  caffin£> 
loop 


where  D>0,  k  >0,  k  >0  are  constant  gains,  and  ,  are,  respectively, 
g  c  c 

the  commanded  and  actual  gimbal  rates.  Substituting  Eqs.  (8)  and  (9)  into  (5), 
the  gyrodynamics  of  each  CMG  are  now  given  by: 


Jcf  =  h0  +  (kgD/2h)#  -  k^a  -  kc<r 


(10) 


The  transfer  function  corresponding  to  Eq.  (10)  is  given  by  (l  =  -1): 


iK 


cr/e  = 


io. 


U  +  1W  jW 


where 


K  =  (h  +  (k  D/2h)j/J 

5 

W1  H  kg/J 


(11) 


W, 


=  VJ 


and  the  gyrodamper  bandwidth  is  the  interval  ,  where  “o  *  "o'  "l 
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cr/e 


Fig.  5  Gyrodamper  Transfer  Function 


The  colocated  loop-closure  is  finally  produced  by  using  the  optical  position 
sensor  with  a  classical  rate  circuit  compensation  shown  in  Fig.  6.  The  compen¬ 
sator  generates  a  rate  signal  up  to  7  Hz  to  10  Hz  and  is  then  rolled  off.  Curves 
a  and  b  correspond  to  shifts  in  the  second  pole  so  that  the  relative  damping  of 
the  1  Hz  and  7  Hz  bending  modes  may  be  adjusted.  The  21  Hz  mode  is  not 
controlled . 
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Fig. 


6.3  TOYSAT  EXPERIMENTS 


6.3.1  Introduction 

The  next  three  sections  describe  experiments  performed  in  the  Space  Systems 
Division's  Modal  Lab.  The  experiments  were  designed  to  test  three  of  LMSC's 
basic  control  theories;  optimal  slewing,  Low  Authority  Control,  and  High  Authority 
or  Modern  Modal  Control. 

A  photograph  of  the  test  specimen  is  shown  in  Fig.  1.  The  specimen  consists 
of  a  rigid  block  of  aluminum  in  the  middle,  with  a  long  flexible  bar  of  aluminum 
attached  to  each  side.  The  central  block  is  suspended  from  the  ceiling  by 
thin  wires.  The  control  is  provided  by  two  Electroseis  linear  actuators  which  act 
on  the  rigid  central  body.  For  these  experiments,  the  linear  actuators  are  com¬ 
manded  to  give  equal  but  opposite  forces ,  resulting  in  a  pure  torque  on  the  central 
body.  There  are  two  sets  of  sensors  installed.  On  both  ends  of  the  flexible  bar 
accelerometers  provide  information  on  the  tip  acceleration.  Connected  in  parallel 
with  the  Electroseis  linear  actuators  are  linear  potentiometers  which  give  the 
amount  of  rotation  of  the  central  body.  Small  Ling  shakers  are  also  attached  to 
each  end  of  the  flexible  bar,  however  they  are  not  used  for  control  purposes; 
their  only  function  is  to  provide  a  weight  at  the  end  of  the  flexible  bar  so  as  to 
lower  the  frequencies  of  the  natural  modes.  For  each  experiment  the  control 
commands  are  generated  differently.  They  will  be  described  separately  in  each 
section . 

An  analytical  model  was  developed  in  order  to  implement  the  control  schemes. 

The  model  was  used  to  generate  a  modal  model;  a  truncated  version  is  shown  in 
Table  1.  The  values  predicted  by  the  analytical  model  were  verified  by 
performing  a  sine  sweep  of  the  test  speciment.  Agreement  between  analytically 
predicted  and  experimentally  obtained  frequencies  was  quite  good.  The  biggest 
difference  in  frequency  was  around  5  percent.  The  analytical  model  was  then 
adjusted  to  agree  exactly  with  the  experimental  results. 
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TOYSAT  Test  Vehicle 


Table  1 


SYSTEM  MODAL  EQUATIONS 


0  is  rigid  body  angle 

q.  are  modal  amplitudes 

w.  are  modal  frequencies 

J  is  the  total  moment  of  inertia 

Tj.  are  the  control  influence  coefficients 

T  is  the  control  torque 

Values: 

vil  =  2.62  Hz 

w  =  9.4  Hz 

“  2 

J  =  135.05  lb-in-sec 

7?1(0)  =  1.53 
^(0)  =  0.78 

6.3.2  Optimal  Slew  Maneuver 

The  Optimal  Slew  as  presented  here  is  an  open  loop  torque  command  to  cause  a 
rigid  body  angular  displacement  while  leaving  the  first  bending  mode  quiescent 
at  the  end  of  the  maneuver.  In  general,  more  than  one  mode  can  be  controlled 
at  the  end  of  the  maneuver,  but  in  this  case  only  the  first  mode  is  considered. 
"Optimal"  in  the  title  is  justified  because  the  particular  maneuvers  shown  here 
minimize  the  control  energy  used  for  the  maneuver  as  well  as  accomplishing  the 
desired  objectives.  Briefly,  the  control  profile  is  a  solution  of  the  calculus  of 
variations  problem  shown  in  equation  (1). 
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min 

u 


(1) 


dt 


subject  to  x(0)  =  Xq  ,  x(t^  )  =  0  . 

Here  u  is  the  control  torque,  and  x  is  a  state  vector  composed  of  the  rigid 
body  angle,  the  first  modal  amplitude,  and  both  of  their  rates  of  change.  For 
a  more  detailed  description  see  "Optimal  Feedback  Maneuvering  of  Flexible  Space¬ 
craft"  by  John  A.  Breakwell,  presented  at  the  AIAA/AAS  Astrodynamics  Specialist 
Conference,  Provincetown ,  Massachusetts,  June  1979. 

For  these  experiments,  the  commands  to  the  linear  actuators  were  provided  by  a 
PDP1145  digital  computer.  The  desired  profiles  were  approximated  by  small  seg¬ 
ments  of  constant  output.  Each  command  was  divided  into  one  hundred  such 
segments.  The  output  from  the  computer  was  fed  through  a  digital  to  analog 
converter,  and  then  smoothed  by  an  eight  pole  butterworth  filter  set  at  100 rad/sec 
before  being  input  to  the  actuators. 

The  effectiveness  of  this  slewing  technique  is  demonstrated  in  Figs.  2  and  3.  All 
slews  are  through  10  degrees.  For  comparison  purposes,  the  responses  to  square 
wave  torque  commands  are  also  shown.  The  tip  acceleration  output  is  filtered  to 
show  frequency  content  between  0.1  and  5  Hz.  Both  optimal  slews  show  almost 
no  residuals  of  the  tip  accelerations  at  the  end  of  the  maneuvers.  What  residual 
there  is  can  be  attributed  to  the  1  Hz  symmetric  mode,  which  was  not  controlled, 
and  only  insignificantly  excited. 

6.3.3  Low  Authority  Control 

This  version  of  Low  Authority  Control  takes  the  form  of  a  collocated  rate  feedback. 
The  rate  signal  is  obtained  by  differentiating  the  difference  of  the  linear  poten¬ 
tiometers  attached  to  the  central  body.  The  signal  is  then  rolled  off  twice  for 
smoothing.  Initially  the  roll  off  is  set  at  10.5  Hz,  but  a  case  is  also  described 
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Fig.  2  Comparison  of  Rectangular  and  Optimal  Slew  Command  Results 
Command  Period  =  1.0  Second 


Fig.  4  Low  Authority  Control  Experiment  Block  Diagram 
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Low  Authority  Control  TOYSAT  Angular  Deflection  (Degrees) 
Response  to  Impact  Excitation 


GAIK] 


with  the  double  roll-off  at  18.9  Hz.  The  signal  is  then  multiplied  by  some  gain, 
and  input  to  the  linear  actuators.  A  block  diagram  of  the  system  is  shown  in 
Fig.  4.  Rather  than  picking  the  gains  through  the  use  of  an  analytical  model, 
the  testing  capabilities  of  the  Modal  Lab  were  used.  A  sine  sweep  was  performed 
to  obtain  a  Bode  plot  of  the  transfer  function  between  the  output  (the  smoothed 
rate  signal)  and  the  input  (the  signal  to  the  linear  actuators).  The  Bode  plot 
is  shown  in  Fig.  5.  Classical  techniques  were  then  used  to  predict  the  highest 
gain  which  would  yield  stability.  The  prediction  was  a  gain  of  10.  Figures 
6  and  7  show  the  result  of  the  closed  loop  system  acting  on  an  excited  TOYSAT. 
The  gains  were  varied  from  0  to  20.  Because  of  the  rate-feedback  nature  of  the 
feedback  loop,  the  angular  displacement  history  exhibits  some  drift,  however 
the  tip  acceleration  histories  show  significant  damping.  The  resulting  modal 
dampings  are  shown  in  Table  2. 

6.3.4  High  Authority  Control 

The  final  TOYSAT  experiment  is  a  test  of  both  Modern  Modal  Control  and  digital 
control.  Measurements  of  the  tip  accelerations  and  the  central  angular  deflection 
are  processed  by  the  PDP  1145  digital  computer  to  provide  estimates  of  the  instan¬ 
taneous  modal  amplitudes.  These  estimates  are  then  multiplied  by  gains  and  fed 
back  to  the  linear  actuators.  In  addition  to  the  PDP  1145,  analog  to  digital  and 
digital  to  analog  converters  are  used  to  convert  the  measurement  signals  and  the 
control  commands  respectively.  Before  the  measurements  are  converted  to  digital 
signals,  they  are  passed  through  a  100  Hz  eight  pole  Butterworth  filter  which 
acts  as  an  anti-aliasing  filter.  Also,  before  the  commands  are  fed  to  the  linear 
actuators,  they  are  smoothed  by  a  similar  filter. 

The  formulation  of  the  estimator  and  controller  is  shown  in  Figs.  8  and  9.  The 
method  used  for  dealing  with  the  acceleration  measurements  makes  the  filter  sub- 
optimal.  The  optimal,  or  Kalman  formulation,  is  applicable  only  when  measure¬ 
ment  and  state  disturbances  are  uncorrelated,  which  is  not  the  case  here.  The 
delay  caused  by  the  two  Butterworth  filters  was  quite  significant,  amounting  to 
four  time  steps.  The  usual  way  of  dealing  with  this  sort  of  delay 
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Table  2 


TOYSAT  LOW  AUTHORITY  CONTROL  EXPERIMENTAL  RESULTS 


Mode 

Resonant 

Viscous 

Feedback 

Description 

Frequency  (Hz) 

Damping  (%) 

Gain,  K 

f 

=  10.5  Hz 

c 

First 

2.600 

0.3 

0 

Asymmetric 

2.707 

8.3 

10 

2.714 

18.7 

20 

Second 

8.97 

2.3 

0 

Symmetric 

9.39 

1.4 

10 

9.00 

1.2 

20 

Second 

9.37 

1.6 

0 

Asymmetric 

9.56 

2.5 

10 

9.60 

3.4 

20 

f 

c 

=  18.9  Hz 

First 

2.60 

0.3 

0 

Symmetric 

2.63 

4.4 

10 

Second 

8.97 

2.3 

0 

Symmetric 

8.99 

2.8 

10 

Second 

9.37 

1.6 

0 

Asymmetric 

9.52 

2.7 

10 

requires  augmenting  the  state  vector.  Speed  limitations  on  the  PDP  prevent 
this  alternative.'  A  simpler  approach  is  used;  it's  described  by  Fig.  10. 

The  performance  of  this  digital  regulator  is  shown  in  Figs.  11,  12,  and  13. 

The  TOYSAT  was  manually  excited,  then  the  controller  was  turned  on.  Most  of 
the  tip  acceleration  dies  out  within  one  second  -  much  faster  than  with  the  Low 
Authority  controller.  There  is  some  drift  apparent  in  the  angular  displacement. 
This  is  probably  caused  by  the  poor  DC  response  of  the  linear  actuators. 
Figures  14,  15,  and  16  show  the  output  of  the  estimator  which  was  used  to 
generate  the  actuator  commands. 
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High  Authority  Control  Formulation 
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Fig.  9  High  Authority  Control  Formulation 


Fig.  10  Group  Delay 


Fig.  11  Digital  High  Authority  Control  of  Rigid  Body  and  First 
Two  Asymmetric  Modes  of  TOYSAT  Specimen 
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Fig.  12  Experimental  Angular  Deflection 
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Fig.  13  Torque  Command 


The  preceding  results  were  generated  rather  quickly.  It  is  anticipated  that 
careful  tuning  of  the  gains  could  result  in  even  better  damping  of  the  modal 
disturbances. 


\ 
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6.4  MICROVIBRATION  SENSING 

6.4.1  Introduction 

A  He-Ne  laser  Vibration  Sensor  has  been  developed  which  features  a  digital  out¬ 
put  for  computational  convenience,  and  which  complements  conventional  vibration 
sensors  (i.e.  accelerometers)  by  sensing  vibratory  events  at  low  frequencies  from 
DC  to  beyond  50  Hz.  Vibration  amplitude  resolution  of  the  sensor  is  0.08  n  m; 
maximum  amplitude  and  frequency  product  is  presently  limited  to  0.05  m'Hz  for 
2  MHz  electronic  bandwidth.  For  example,  the  maximum  measurable  vibration  am¬ 
plitude  for  a  25  Hz  vibration  is  2  mm.  The  time  delay  of  the  sensor  output  to  the 
actual  vibration  is  less  than  1  n  sec  which  is  nearly  perfect  for  measuring  the  dy¬ 
namics  of  structures  and  vibration  sensing  for  the  dynamic  damping  of  structures 
(active  control  of  structures).  By  electronically  splitting  the  laser  beam  using  a 
Bragg  cell  it  is  possible  to  simultaneously  sample  and  hence  monitor  a  large  num¬ 
ber  of  points  to  which  retroreflectors  have  been  affixed.  Although  the  laboratory 
Vibration  Sensor  employed  but  two  channels,  it  exhibited  the  basis  for  continu¬ 
ously  sensing  more  than  50  independent  vibrating  targets. 

The  Vibration  Sensor  employs  a  low  power  He-Ne  laser  and  two  Bragg  cells,  one 
to  provide  a  heterodyne  offset  frequency  for  use  as  a  local  oscillator,  the  other 
for  generating  multiple  beams  for  various  targets.  Optical  path  variations  due  to 
the  vibratory  motion  of  each  target  is  measured  by  comparing  zero-crossings  of 
the  local  oscillator  signal  with  those  of  each  target,  identifiable  by  a  specific  Bragg 
frequency.  The  various  differential  zero-crossings  provide  the  basis  for  the  digi¬ 
tal  output  of  the  Vibration  Sensor.  Experimental  results  compare  well  with  theo¬ 
retical  predictions. 
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6.4.2  Theory 


The  Vibration  Sensor  conceptual  layout  diagram  is  illustrated  in  Fig.  1.  The 

first  Bragg  cell  translates  the  laser  output  frequency  from  v  to  v  +  v  .  The 

o  OS 

frequency  translated  beam  is  directed  by  the  second  Bragg  cell  to  two  target 
retroreflectors  which  correspond  to  two  driving  frequencies  u.  and  v0  , 
respectively.  The  laser  beam  incident  upon  the  channel  1  target  at  range  yj 
has  a  frequency  uQ+  which  is  different  from  frequency  vq  +  v&  +  v 

for  the  laser  beam  on  channel  2  target  at  range  y2(t).  The  returning  signals 
vn  +  va  +  v,  +  2 y’./A  and  u  +  vn  +  +  2y /A  (where  X  is  the  laser  wave- 

length)  are  being  shifted  in  frequency,  when  passed  through  the  second  Bragg 
cell  again ,  by  v  ^  and  ^  ,  respectively .  A  photodiode  detector  is  used  to  sense 
the  heterodyne  beat  between  the  optical  local  oscillator  (LO)  beam,  which  is  de¬ 
rived  from  the  laser  output  directly,  and  the  vibration  sensing  beams. 


The  electric  field  of  the  laser  beam  on  the  photo- detector  from  channel  1  and 
channel  2  targets  can  be  expressed  as 


E  =  E  .  +  E  „ 
s  si  s2 

=  exp  ji 
E2  exp  ji 


2jt(  vq  +  va  +  2  rx)  t  +  47ryx  (t)/A 
2H  v  +va  +  2  r  )  t  +  4  ;ry  (t) /A 


The  electric  field  of  the  optical  LO  can  be  expressed  as 

EL  =  Ep  exp  ji2  7T[,otj 


(1) 


(2) 


The  resultant  electric  field  on  the  photo  detector  is  the  sum  of  equations 
( 1)  and  (2)  or 


E  = 


E  +  Et 
s  L 


=  E-  exp  /i  [2  ?r(  v  +  va  +  2  iO  t  +  4?ry  (t)/xl  l 

f  r  J(n 

+  E2  exp  |i  |^2  TT  (  Vq  +  I'a  +  2  1^)  t  +  47ry2  ( t )  /  X  J  ^ 
+  E^  exp  ji  2rr^ot| 


Optical  power  received  by  the  photodetector  is  proportional  to  the  product  of 
equation  (3)  and  its  complex  conjugate  or 


PaEE*  =  Ej2  +  E22  +  E£  2  = 

2E1  E£  cos  ^  2 7T  (  va  +  2^1)t  +  47ry1  (t)/\J 

+2E2  E^cos  2  tt(  u&  +  2^2)t  +  4^y2  (t)Vxj  (4) 

+2Ej  E2  cos  £4  tt(  -  v^)  t  +  4?ry1  (t)/A  -  4?ry2  (t)/xj 


Since  »E^,  E^,  equation  (4)  can  be  simplified  to 


P 


P1  +  P2  +  Pf 


+  PSI  +  P 


s2 


where 


DC  component  of  signal  power 
Optical  LO  power  »Pj,  P2 

cos  ^2 7T(  v&  +  2v^)X  +  4Try1(t)/\j 
2/P2  Pp  cos  |^2ir (  i>a  +  2  v2)  t  +  47ry2(t)/xJ 


(5) 

(6) 


The  AC  component  of  photo  current  is 

Zs  ~  ‘si  +  *s2 

=  Tu7  <psl  +  Ps2) 


(7) 


=  P{  cos  2?r(^a  +  2  ^1  T  +  4 Try x  (t) /xj 

+  /p2  P^  cos  ^2?r(  Uq  +  2  v2)  t  +  4?ry2  (t)/xj 
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Here  we  have  used  equation  (6)  and  (7), 
where 

r)  =  quantum  efficiency  of  the  photo  diode  detector 
q  =  charge  of  electron 
hy  =  laser  photon  energy 


The  mean  square  noise  current  can  be  approximated  as 


(8) 


where 

B  =  electronic  bandwidth 

k  =  Boltzman's  constant 

T  =  equivalent  temperature  of  preamplifier 
d 

=  resistance  of  detector  load  resistor 
The  power  signal-to-noise  ratio  of  the  sensor  is 


(9) 


for  channel  1  target  and 


S_ 

N 


P2P£ 


4kT 


(10) 


for  channel  2  target. 
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Assume  a  0.  4  mW  single  frequency  He-Ne  laser  is  used  to  measure  vibrations  of 
50  targets,  the  signal-to-noise  ratio  of  each  channel  can  be  calculated  according 
to  the  typical  values 


n  =  o.8 

-10 

q  =  1. 6  x  10  3  coul 

hv  =  3. 13  x  10  19  joul 

P 'jj  =  10  4  watt 

Pj  =  5  x  10  9  watt 

B  =  2  x  106  Hz 

-23 

k  =  1.38  x  10  0  j/°k 

Ta  =  5  96°k  (for  an  amplifier  with  3  dB  noise  figure) 

=  50  Ohms 


Yields  the  result  according  to  Equation  (9) 


|0.8x  1.6  x  10~19\ 
\  3. 13  x  10~19  / 


x  5  x  10  9  x  10  4 


=  125 


2  x  10l 


( 


2  x  0.8  x  2.56  x  10  38  x  10~4 


4  x  1.38  x  10  23  x  5  96 


3.13  x  10 


-9 


50 


) 


This  signal-to-noise  ratio  is  more  than  what  is  required  by  the  digital  signal  pro¬ 
cessor  to  process  the  signal.  For  a  constant  signal-to-noise  ratio,  the  maximum 
number  of  channels  can  be  measured  by  the  sensor  is  proportional  to  the  laser 
output  power. 

The  phase  modulation  expression  in  Equation  (1)  for  channel  1  signal 

exp  |i  [ 2n(vQ  +  v&  +  2 u^)  t  +  4-rry1  (t)/A]j  can  also  be  expressed  in  a  frequency 

modulation  form  such  as  exp  ji  [2^(1^  +  v&  +  2v ^  +  2  y.  (t)/\)  t  +  6^]j  . 

where  0^  is  constant.  The  Doppler  frequency,  =  2y^  (t ) / A  ,  for  a  sinusoidal 
vibration  can  be  written  as  v.,  =  4ttv  .y. .  sin  (2-au  .t)/\,  where  v  .  is  the 
modulation  frequency  and  y^  is  the  amplitude  of  vibration.  For  a  2  MHz  electronics 
bandwidth,  the  peak  Doppler  frequency  is  limited  to  1  MHz  or  4  w v  ^y ^ j /A  =  108 
Hz  or  =  0.05  m-Hz.  Therefore,  the  maximum  amplitude  and  frequency  pro¬ 

duct  of  the  vibratory  targets  is  limited  to  0.05  m-Hz  for  2  MHz  electronic  bandwidth. 


In  order  to  separate  two  target  signals  from  the  photodetector  output  we  can  first 

mix  with  cos  [2ir(v  -  v  )  t]  to  obtain  cos  [2^(2^,+  v)t  +  47ry./X] 

a  r  i  r  i 

and  cos  [ 2 7r ( 2  +  v  )  t  +  477y2/A],  and  then  mix  these  with  cos  [277(21^)  t] 

to  obtain  cos  ^2  77  ^  t  +  4  7T  y  ^  /  X  J  or  with  cos  [277  (2  V  t]  to  obtain 

cos  ptfj,  t  +  ixy^lX  j. 

Here  v  is  the  frequency  of  the  reference  signal  which  is  4  MHz  in  this  case.  The 
signal  cos  [277  (  y  )  t  ]  is  a  result  of  mixing  the  reference  signal  with  the 

'“irst  Bragg  cell  driving  signal  and  bandpassing  at  v a  -  v^.  Signals  cos  [ 2 7r( 2  v t] 
and  cos  [2tt  (2  t]  are  derived  from  frequency  doublers  and 
in  turn  from  oscillators  at  frequencies  v ^  and  v 2,  respectively.  Signals  of  both 
channels  are  centered  at  frequency  vv  but  separately  passed  by  each  bandpass 
filter.  The  digital  signal  processor  takes  the  reference  signal  and  the  output  of 
two  channels  from  the  bandpass  filters  to  process  the  vibration  information  of 
targets.  By  comparing  the  zero-crossings  of  cos  (2  77  i>  t)  with  cos  (2  + 

4  77  Yj/A)  and  cos  (2vvv\  +  477y2/A)  the  displacement  Ayx  and  Ay2  can  easily 
be  determined.  For  every  zero-crossing  count  difference  there  is  a  phase  change 
of  2  77  radians  or  a  displacement  of  A/2  =  0.32  (im.  If  four  reference  signals 

COS  (2  77  V  rt)  ,  COS  (2  77  V  _X  +  77/2)  ,  COS  (2  77  V  ft  +  7T  ),  and  COS  (2  77  p  +  377/ 2) 

are  used  for  •  zero-crossing  comparison ,  the  displacement  resolution  would  be 
A  /8  or  0.08  p.m. 
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6.4.3  Digital  Signal  Processor 


A  digital  electronic  circuit  has  been  designed  to  measure  the  total  phase  difference 
between  the  reference  signal  (4  MHz)  and  the  phase  modulated  heterodyne  signal 
returning  from  targets.  A  zero-crossing  counting  technique  is  used  to  determine 
the  phase  difference  in  a  multiple  of  2  n  radians  which  corresponds  to  a  displace¬ 
ment  resolution  of  about  0.32  p. m.  A  fine  phase  measurement  system  based  on 
the  zero-crossing  timing  has  also  been  built  in  order  to  achieve  a  0.08  fim 
resolution . 

CLOCK  CIRCUIT 


All  the  timing  signals  needed  for  the  digital  circuit  are  derived  from  a  16  MHz 
oscillator.  With  a  4  bit  binary  counter  and  a  4  to  16  line  decoder,  the  clock  cir¬ 
cuit  is  capable  of  generating  16  disjointed  clocks.  Each  of  them  has  a  frequency 
of  1  MHz  and  pulse  width  of  about  60  ms.  The  clocks  used  in  the  signal  processor 
are  illustrated  in  Fig.  2,  where  0  is  used  in  the  coarse  phase  measurement;  0^, 

0  2,0^  together  with  0^,  are  used  in  the  fine  phase  measurement;  0,-  is  used  as 
the  sampling  signal  for  the  output.  In  addition  to  the  clocks  mentioned  above, 
the  clock  circuit  also  provide  a  4  MHz  reference  signal  for  use  as  the  carrier  for 
the  phase  modulated  target  signals. 

COARSE  PHASE  MEASUREMENT 


The  coarse  phase  measurement  circuit  consists  of  a  reference  channel  and  as  many 
target  channels  as  needed.  Figure  3  shows  the  functional  block  diagram  for  a 
two-target  system.  Additional  channels  can  easily  be  added.  As  shown  in  Fig.  3, 
the  reference  channel  is  just  a  14  bit  binary  counter  which  continuously  counts 
the  number  of  zero-crossings  (or  the  number  of  2t radian  phase  change)  of  the 
4  MHz  reference  signal.  A  14  bit  latch  samples  the  content  of  the  counter  at  the 


Channel  1 

J 


Channel  2 


Fig.  3  Vibration  Sensor  Coarse  Phase  Measurement  System 
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at  the  falling  edge  of  clock  0  The  output  of  the  latch  corresponds  to  the  total 
phase  change  of  the  reference  signal  since  the  beginning  of  counting  process. 
Signals  from  each  target  are  first  shaped  by  a  zero-crossing  detector  then  passed 
through  a  circuit  identical  to  the  reference  channel  to  generate  a  binary  number 
corresponding  to  the  total  phase  change  in  the  target  signal.  This  number  together 
with  the  one  provided  by  the  reference  channel  is  fed  into  a  binary  subtractor  to 
calculate  the  total  phase  deviation  of  the  target  signal  due  to  the  displacement 
of  the  target. 


FINE  PHASE  MEASUREMENT 

To  achieve  a  phase  resolution  better  than  2  n  radians,  a  circuit  has  been  designed 
to  determine  the  timing  of  zero-crossing  of  the  target  signal.  As  illustrated  in 
Fig.  4,  the  least  significant  bit  of  the  target  counter  in  each  coarse  system  is 
sampled  at  4  different  instances  which  are  tt/2  radians  apart  in  terms  of  4  MHz 
frequency.  The  output  of  the  latches  are  then  processed  by  a  simple  logical 
circuit  to  determine  when  the  zero-crossing  occurs  and  to  encode  this  information 
into  2  bits.  These  two  bits  together  with  the  14  bits  from  the  coarse  system  are 
sampled  at  the  falling  edge  of  clock  0^  to  provide  a  steady  16  bit  output  for  the 
vibration  sensor  system.  A  graphic  simulation  of  all  the  input  and  output  signals 
of  the  digital  signal  processor  is  summarized  and  shown  in  Fig.  5 
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6.4.4  Experimental  Apparatus 


The  Vibration  Sensor  laboratory  setup  and  the  breadboard  layout  diagram  are 
shown  in  Figs.  6  and  7,  respectively.  The  Channel  1  target  is  a  1-cm  retro- 
reflector  driven  by  a  shaker  located  on  a  separate  Table  3  meters  away  from  the 
sensor  setup.  The  Channel  2  target  is  a  PZT  driven  mirror  set  on  the  same  plat¬ 
form  with  other  optical  components.  A  Hughes  Model  3176H  He-Ne  laser  with 
4  mW  output  at  632.8  nm  is  used  in  the  system.  The  reflectivity  and  transmission 
of  all  the  beamsplitters  used  in  the  system  are  close  to  0.5.  The  efficiency  of  the 
first  Bragg  cell,  e. ,  is  0.8.  The  second  Bragg  cell  efficiency  can  be  expressed 
as  =  0.8/N,  where  N  is  the  number  of  channels.  Because  the  Bragg  cell  can 
only  take  3  watts  RF  power  to  reach  80  percent  efficiency,  the  RF  driving  power 
for  each  channel  is  3/N  watts.  Therefore,  the  Bragg  cell  efficiency  for  each 
channel  is  c2  =  0.8/N.  Two  cylindrical  telescopes  are  used  in  conjunction  with 
the  Bragg  cell  to  generate  a  large  number  of  target  sensing  beams  for  multiple 
targets  vibration  sensing  purposes.  Only  one  beam  expander  is  used  in  Channel  1. 
The  laser  beam  diameter  after  the  beam  expander  is  approximately  1  cm.  The 
return  signal  power  onto  the  photodetector  can  be  calculated  by 

Ps  -  Po  X  T1  x  6i  X  T2  X-  e2  X  €2  x  Tg 

,4  x  10  3  x  0.5  x  0.8  x  .08/N  x  0.8/N  x  0.5  x  0.5 
*  0.128  x  10~3/N2  Watt 

where  T's  and  R's  are  the  transmission  and  reflectance  of  beamsplitters  shown  in 
Fig.  7  and  PQ  =  4  mW  is  the  laser  output  power.  For  50  channels,  N  =  50,  the 
signal  power  is  pg  «  5  x  10®  Watt  for  pQ  =  4  mW  or  pg  ~  5  x  10  ®  Watt  for 
po  =  0.4  mW.  In  our  experiment  N  =  2  and  the  calculated  power  level  is 
32  x  10“ 6  Watt. 


Fig.  7  Vibration  Sensor  (2-Channel)  Breadboard  Layout  Diagram 


The  optical  LO  power  onto  the  detector  is  =  Pq  x  Rj  x  Rg  which  is 

P^i  a  4  x  l(f3  x  0.5  x  0.05  =  10-3  Watt  for  Pq  =  4  raW  or  Pj  =  10~4  Watt  for 

P  =  0.4  mW. 
o 

Two  K&L  Model  5B53-4/2-B/B  bandpass  filters  centered  at  4  MHz,  one  for  each 
channel,  are  used  in  the  system  to  restrict  the  input  signal  to  the  digital  signal 
processor  to  a  bandwidth  of  2  MHz.  A  two-channel  digital  signal  processor  has 
been  assembled  and  tested .  Due  to  the  availability  of  digital  IC  components ,  we 
use  16  bits  instead  of  14  bits  in  the  entire  coarse  measurement  system  thus  pro¬ 
ducing  two  18  bit  outputs.  Two  16  bit  Digital- to- Analog  (D/A)  converters  are 
used  to  convert  the  lower  16  bits  of  each  output  into  analog  form  for  the  display 
purpose. 
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6.4.5  Experimental  Results 


The  oscilloscope  trace  of  the  displacement  of  two  vibratory  targets  is  shown  in 
Figs.  8  and  9.  In  Fig.  8,  the  upper  trace  is  the  sensor  output  for  Channel  1 
target  vibrating  at  30  Hz  and  an  amplitude  of  1.5  mm.  The  middle  trace  is  the 
sensor  output  for  Channel  2  target  which  measures  a  60  Hz  vibration  at  an  ampli¬ 
tude  of  0.9  pm.  It  is  noticeable  that  the  "stair-like"  waveform  is  a  result  of 
digital  signal  processing.  Each  step  of  the  stair  represents  0.08  pm  displacement 
of  target  which  is  the  resolution  of  the  present  system .  The  lower  trace  represents 
the  driving  signal  to  the  PZT  for  Channel  2  target.  Comparing  the  output  of  the 
Vibration  Sensor  with  the  driving  signal  of  the  target  mirror  indicates  a  time 
delay  of  about  1  nsec  between  the  sensor  output  and  the  actual  vibration.  Of 
this,  about  500  nsec  is  contributed  by  the  digital  circuitry  (between  the  falling 
edges  of  input  sampling  clock  <p  ^  and  output  sampling  clock  0,.) ,  the  rest  of  it  is 
due  to  the  settling  time  of  the  D/A  converter.  Power  consumption  for  the  two-channel 
digital  signal  processing  system  is  about  10  watts  with  standard  TTL  components. 

Although  the  output  of  the  sensor  system  closely  follows  the  actual  vibration  in 
general,  we  found  a  large  number  of  glitches  presented  especially  when  the  ampli¬ 
tude  of  the  vibration  is  small.  These  glitches  occurs  at  certain  intervals  of  tar¬ 
get  displacement  corresponding  to  even  2  7 r,  4  7r ,  6  7r ,  and  8tt  radian  phase 
changes.  Most  of  these  glitches  have  been  identified  as  being  caused  by  the  sam¬ 
pling  of  target  signal  counters  while  the  counters  are  changing  states.  We  have 
done  experiments  to  sense  such  ill-conditioned  triggering  cases  in  the  coarse 
system  and  delay  the  latch  action  whenever  it  occurs.  This  scheme  seems  to  be 
able  to  eliminate  glitches  at  2  77,  4n  ,  and  67T  intervals  but  those  with  an  8  it  interval 
still  remain.  At  this  time  no  effort  is  made  to  identify  the  cause  of  this  problem. 

But  we  feel  this  can  be  solved  with  confidence  after  a  careful  examination  of  the 
digital  signal  processing  electronics.  The  DC  response  of  the  sensor  is  shown  in 
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Fig.  8  Vibration  Sensor  Output  for  Two  Vibratory  Targets  - 

Kxhibition  of  Amplitude  Range  and  Sensitivity  of  Sensor 


Fig.  9  Vibration  Sensor  Output  for  Two  Vibratory  Target 
Exhibition  of  AC  and  DC  Response  of  Sensor 


the  upper  trace  of  Fig.  9.  The  only  difference  between  the  Fig.  8  and  Fig.  9 
experiments  is  the  driving  voltage  applied  to  the  shaker  for  channel  1  target. 

In  the  Fig.  9  experiment,  a  0.5  Hz  square  wave  is  applied  to  the  shaker.  The 
sensor  measures  the  steady  state  DC  displacement  (~  1  mm)  as  well  as  the  tran¬ 
sient  behavior  of  the  shaker. 

As  mentioned  earlier,  the  signal  processor  provides  a  18  bit  digital  output  in 
2's  complement  form.  In  other  words,  the  output  digits  represent  a  decimal 
number  between  -131072  to  +131071.  With  the  least  significant  bit  corresponding 
to  0.08  pm  in  displacement,  this  means  a  displacement  range  from  -10.485  mm  to 
+10.485  mm.  It  is  worthwhile  to  mention  that  wider  range  can  be  achieved  by 
simply  using  longer  words  in  the  signal  processing  electronics  and  the  only 
limitation  is  that  the  vibration  amplitude  and  frequency  product  must  be  under 
0.05  m-Hz.  Since  we  are  only  using  the  lower  16  bits  in  the  display,  the  range 
shown  on  the  oscilloscope  is  further  restricted  to  -2.62  mm  to  +2.62  mm.  Larger 
displacement  causes  a  sudden  jump  from  the  largest  representable  value  to  the 
lowest  value  as  can  be  seen  in  Fig.  10. 

A  reset  button  is  built  in  the  electronics  to  signal  the  beginning  of  the  counting 
process.  Hence  the  output  of  the  electronics  represents  the  displacement  of  the 
target  relative  to  the  target  position  when  the  reset  button  is  released.  At  the 
present  time  we  have  no  control  on  the  zero  displacement  position  relative  to  the 
vibration  waveform .  This  causes  a  DC  component  associated  with  the  output  of 
the  sensor. 

The  laser  used  in  the  system  outputs  three  longitudinal  modes  spaced  by  500  MHz 
which  is  shown  in  Fig.  11.  Although  the  unstable  laser  is  used  in  such  a  sensi¬ 
tive  system  we  can  still  obtain  reasonably  good  results.  Because  of  the  500  MHz 
frequency  spacing  between  laser  modes  and  the  destructive  interference  between 
them ,  we  have  observed  that  for  every  30  cm  displacement  of  target  along  the 
optical  path  of  the  vibration  sensing  beam  there  is  a  node  point  for  the  returning 
signal  which  limit  the  full  range  application  of  the  sensor.  The  single  frequency 
laser  will  be  a  solution  to  that  problem. 
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Fig.  10  Output  Discontinuity  of  Vibration  Sensor 
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Fig.  11  Optical  Spectrum  Analyzer  Output  Shows  3  Longitudinal  Modes 
of  the  Laser  Used  in  the  Vibration  Sensor  Experiment 


6  59 


6.4.6  Conclusions  and  Recommendations 


After  this  very  limited  effort  toward  the  development  of  the  Vibration  Sensor, 
we  are  satisfied  with  the  outcome  of  this  research.  The  concept  of  using  zero¬ 
crossing  counting  technique  in  the  sensing  of  low  frequency,  small  amplitude 
vibration  has  been  proved.  The  capability  of  the  sensor  for  simultaneous  sens¬ 
ing  of  multiple  targets  and  wide  frequency  and  amplitude  coverage  of  the  vibra¬ 
tion  has  been  demonstrated.  The  remaining  problems  which  need  to  be  done  in 
the  future  concerning  the  Vibration  Sensor  are  listed  below: 

1.  Design  and  fabricate  or  purchase  a  stable  single  mode  single 
frequency  He-Ne  laser  with  an  output  of  0.4  mW  at  632.8  nm. 

2.  Remove  the  glitches  from  the  sensor  output. 

3.  Develop  a  second  generation  signal  processor  to  handle  larger 
amplitude  and  frequency  range  of  vibration  to  provide  improved 
resolution,  to  automatically  adjust  the  system  zero  to  the  center 
of  vibration,  and  to  derive  velocity  /acceleration  information  of 
targets. 

4.  Incorporate  the  microcomputer  and  the  graphical  display  into 
the  system  for  real  time  display /control  of  vibrating  targets. 

5.  Fabricate  or  purchase  compact  optics  and  electronics  to  minimize 
the  size  of  the  sensor. 

6.  Procure  better  quality  optics  such  as  cylindrical  lenses,  beam¬ 
splitters,  Bragg  cells,  etc.  in  order  to  reduce  the  optical  loss. 

7.  Verify  the  minimum  required  signal-to-noise  ration  experimentally. 

8.  Isolate  the  EMI  problem. 

9.  Develop  and  demonstrate  a  10-channel  Vibration  Sensor  as  a  basis 
for  the  design  of  sensor  with  more  channels. 
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Section  7 

CONCLUSIONS  AND  FUTURE  WORK 

7.1  CONCLUSIONS 

The  integrated  control  design  methodology  presented  in  the  previous  sections 
has  significant  promise  for  large  space  structure  applications.  The  approach  is 
geared  towards  addressing  specific  issues  such  as  low  damping  ratios,  truncated 
modes  and  finite  bandwidth  actuators  and  sensors.  The  control  laws  have  been 
optimized  for  minimal  complexity  to  aid  implementation  on  space  based  systems. 
Controller  robustness  with  respect  to  modeling  errors  and  truncation  effects  are 
discussed. 

Application  of  the  methodology  to  three  simple  structures  illustrates  that  satis¬ 
factory  closed  loop  systems  can  be  obtained.  The  procedure,  however,  is  not 
automatic  but  requires  iteration  based  on  designer's  judgement.  The  techniques 
developed  minimize  the  effects  of  truncation  errors  on  control  and  observation 
spillover ,  but  the  stability  of  the  closed-loop  system  must  be  examined  on  a 
case  by  case  basis.  It  is  also  the  belief  of  the  authors  that  absolute  stability 
will  be  impossible  to  assure  with  realistic  sensors  and  actuators  for  unbounded 
modeling  errors. 

7.2  FUTURE  WORK 

The  work  performed  during  this  phase  has  set  a  sound  foundation  for  large 
space  structures  control  design  methodology.  Further  experience  is  required  with 
some  of  the  techniques  already  developed.  Certain  new  methods  must  also  be 
developed  to  provide  a  more  complete  basis  for  control  design  in  complex  space 
structures  with  densely  packed  modes. 


The  following  specific  tasks  will  aid  the  development  of  an  integrated  control 
methodology : 

(a)  The  question  of  mode  truncation  and  reduced  order  models  is  of 
fundamental  importance  in  the  study  of  dynamic  systems  in  general, 
and  large  space  structures  in  particular.  The  roles  of  system  poles, 
zeros,  and  residues  should  be  further  studied.  In  addition,  projection 
methods  should  be  examined.  The  construction  of  models  based  on 
partial  observed  information  (output)  should  be  studied. 

(b)  The  problem  of  collocated  control  should  be  further  studied  to  yield 
a  better  understanding  of  mode  controllability  through  this  type  of 
control.  Actuators /sensors  dynamics  and  their  effect  on  the  design 
should  be  considered  too. 

(c)  Robustness  properties  of  the  system  must  be  studied  for  various 
system  operating  conditions .  For  example ,  actuators ,  sensors 
dynamics,  structural  damping,  etc. 

(d)  The  question  of  actuators  and  sensors  dynamics  should  also  be  con¬ 
sidered.  Sensitivity  of  actuators /sensors  location  to  their  dynamics 
should  be  examined. 

(e)  The  concept  of  frequency-shaped  cost  functionals  and  its  adaptation 
to  a  design  tool  should  be  further  developed.  In  particular,  a  specific 
guideline  for  the  selection  of  cost  functionals  should  be  developed. 

(f)  Design  of  ground  based  experiments  including  validation  procedures 
for  assuring  adequate  robustness  demonstrations  must  be  investigated 
further. 

(g)  Sensor  and  actuator  designs  to  handle  the  micro- vibration  control 
problems  need  further  development. 
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Appendix  A 

MODAL  STRUCTURE  OF  OSCILLATORY  SYSTEMS 


The  oscillatory  system  is  described  by 


x(t)  =  Ax(t)  +  Bud) 


where 


and  A  =  A^  >  0. 
o  o 

Let  U  be  the  n  x  n  matrix  of  eigenvectors  of  AQ,  i.e. , 

A  U  =  UA. 
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Since  Aq  is  p.d.  the  diagonal  elements  of  Aq  are  all  positive ,  i . e . , 


Aq  =  diag  (a^  &2,  ...  ,  an) 


a.  >0Vl4i4n 
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The  eigenvalues  of  the  matrix  F  are  found  from 
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where 


Let 


i  4 


Uo  =  u2 . V 


Vo  =  Ivl’  v2 . v2n1' 


be  the  2n  x  2p  matrices  of  right  and  left  eigenvectors,  respectively,  of  the 

matrix  F  ;  i.e. , 
o 


F  U  =  U  A 
o  O  o  o 


V  F  =  A  V 
o  o  o  o 


where  A  =  diag  (  A.  ,  ...  ,  \  )  is  the  diagonal  matrix  of  eigenvalues  obtained 

o  a 

via  Eq.  (1).  Then,  it  can  be  verified  that 


the  modal  analysis  performed  above,  in  conjunction  with  the  use  of  Euler's 
equations  for  complex  numbers,  can  be  used  to  obtain  the  transition  matrix  in 
closed  form.  Therefore,  the  2n  x  2n  transition  matrix  is  given  by 


where  the  elements  of  the  four  n  x  n  diagonal  blaocks  are  given  by 
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>  i  =  1,  2 . n 


(5) 


For  evaluating  the  controllability  matrix  defined  in  Eq.  (22)  we  need  0,(0,  t) 

z 

rather  than  0  (t,  0)  ;  however,  since 


V°,  t)  =  0-z1(t,O) 


it  can  be  easily  verified  that 

0Z  (0,  t) 


a  (3 

_  y  a . 


(6) 


where  each  of  the  four  n  x  n  blocks  is  diagonal  with  elements  given  by 


Appendix  B 

INVERSE  OPTIMAL  CONTROL  FOR 
GENERALIZED  COLLOCATED  CONTROL 


Ever  since  Kalman  (Ref  1)  formulated  the  inverse  optimal  control  problem  in  1964, 
the  subject  has  been  of  significant  interest  in  control  theory  (Refs  2  through  5). 
The  inverse  optimal  control  theory  attempts  to  derive  a  class  of  performance  in¬ 
dices  which  are  optimized  by  a  prespecified  control  input.  In  its  simplest  form, 
the  state  and  input  penalty  matrices  of  a  quadratic  performance  index  are  deter¬ 
mined  for  linear  feedback  control  of  linear  dynamic  systems  (Refs  1  and  6  Sec.  4.2) . 

Past  research  has  never  exploited  the  full  potential  of  inverse  optimal  control 
theory .  Once  the  control  law  has  been  specified,  obtaining  the  optimizing  cost 
functional  itself  has  little  direct  value.  However,  when  combined  with  known 
and  previously  proven  properties  of  optimal  control  laws,  the  derived  optimizing 
function  aids  understanding  of  the  feedback  and  the  robustness  properties  of 
the  control  law.  An  optimizing  function  is  derived  for  a  specific  class  of  systems 
and  is  used  to  study  controller  properties. 

This  appendix  addresses  purely  oscillatory  systems  defined  by  all  poles  and  zeros 
on  the  imaginary  axis.  Large  space  structures  (LSS)  can  be  approximately  de¬ 
scribed  by  sufch  systems  because  of  their  typically  low  damping  ratios.  Quad¬ 
ratic  optimizing  functions  are  derived  for  control  laws  which  increase  damping 
ratios  in  such  systems. 


Consider  the  linear  dynamic  system  described  by 

x(t)  =  Ax(t)  +  Bu(t),  xeR2n,  ueRm 
y(t)  =  BTx(t),  yeRm 


(1) 


B-l 


where 


0  'll 

-  * 

i  n 

0 

_ j_ _ 

~  T 

B  = 

— 

-A  i  0 

B 

O  I 

.  °J 

(2) 


and  the  n  x  n  matrix  Aq  is  symmetric  and  positive-definite  (pd)  i.e., 


A 


o 


>  0 


(3) 


Next,  consider  the  set  of  admissible  controls,  °lC ,  defined  by 

A  (ue Rm :  u  =  Ky,  Re  [X(A  +  BKBT)]  <  oj  (4) 

where  X  ( • )  stands  for  the  eigenvalue  of  the  argument ,  and  Re(  • )  is  the  real 
part  of  a  complex  number. 

Let  the  objective  function  to  be  minimized  be  quadratic  in  the  state  and  control 
over  an  infinite  time  period,  i.e., 


J(x,u)  £  f°°  xT(t)Qx(t)  +  uT(t)Ru(t)l  dt 
o  L  -* 


(5) 


with 


Q  a  0  ,  R  >  0 


(6) 
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Problem :  (inverse  optimal  control  problem) 

Given  an  admissible  control ,  u*  ,  under  what  conditions  on  K  ,  Q  ,  and  R  is 
this  control  law  optimal,  ie...  when  is 

J(x,u*)  =  min  J(x,u)  (7) 

If  such  an  optimal  control  u*  exists,  then  (Ref  2) 

u*  =  Ky  =  -R-TBTPx  (8) 

where  P  is  the  2n  x  2n  pd  solution  of  the  algebraic  matrix  Riccati  equation 
(ARE)  given  by 


PA  +  ATP  -  PBR-1BTP  +  Q  =  0 


(9) 


The  inverse  optimal  control  problem  is,  therefore,  the  following.  Given  an  ad¬ 
missible  rate  feedback  control  law  find  for  what,  if  any,  weight  matrices  Q  and 
R  ,  the  ARE  in  Eq  (9)  is  solved 


Let  the  2n  x  2n  symmetric  P  and  Q  matrices  be  partitioned  as 


r  i  i 

1 

p  = 

1 

|M 

H* 

-1 - 

1  ’V 
!to 

Q  = 

T  1 

P  IP 

L 12 ;  *22. 

0  i  ^22 

(10) 


Then  from  the  ARE  in  Eq  (  9) ,  we  find 


P11 

AoP22 

-  P12BoR'lBI 

P22 

=  0 

(ID 

12Ao 

T 

+  AoP12 

+  P12BoR'lBo 

P12 

=  «11 

(12) 
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(13) 


P22BoR'lBoP22  '  P12P12  *  «22 
Since  a  rate  feedback  control  law  is  used ,  we  find  from  Eq  ( 8)  that 

Bo  P12  *  ?  —  P12£-^(Bo)  <14) 

Therefore,  P12  can  be  expressed  as 


P 


12 


Nj 


(15) 


where 


(16) 

$  =  >  0  (17) 

The  n  x  n  matrix  P10  as  given  in  Eq  (15)  is  symmetric  and  from  Eq  (17),  psd. 
The  latter  condition  is  required  to  satisfy  Eq  (12)  since  (-AQ)  is  a  stable 
matrix. 


Lemma  1: 


For  a  dynamic  system  under  a  rate  feedback  control  law,  the  n  x  n  matrix  Pl2 
in  Eqs  (11)  through  (13)  has  to  satisfy 

P12  =  0  (18) 


Proof: 


This  is  verified  by  observing  Eq  (13)  where  Q22  has  to  be  at  least  psd. 


B-4 


Since  =  0  ,  the  resulting  2n  x  2n  matrix  P  that  solves  the  ARE  will  be 
a  block  diagonal  matrix.  This  motivates  the  following  form  for  the  matrix  P  to 
be  checked  as  a  possible  solution  of  the  ARE. 

Let 


T  -1  T  T 

P  =  aBCB^B)  B  +  NEN 


(19) 


where 


Ne^(BT),.  NeR2nx(2n  m) 


T 

N  N  =  I 


2n  -  m 


(20) 

(21) 


cneR, 


E  =  ET  >  0  EeR(2n-m)x(2n-m) 


(22) 


It  is  clear  that  P  in  Eq  (19)  is  symmetric.  To  show  that  it  is  also  pd  we  note 
that 


PB  =  aB 


(23) 


PN  =  NE  (24) 

Therefore,  B  and  N  are  right  eigenvectors  of  the  matrix  P  corresponding 
to  positive  eigenvalues;  therefore,  P  is  pd.  Similar  results  are  obtained  for 
left  eigenvectors. 
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Using  the  form  for  P  as  given  in  Eq  (19)  the  ARE  can  be  written  as 


T  9  -1  T 

PA  +  AaP  -  a;  BR  B  +  Q  =  0 


(25) 


Since 


B 


r  T 


we  find  that  to  satisfy  Eq  (20) 


N 


N  eR 
o 


n  x  (n  -  m) 


where  Nq  is  chosen  such  that 


T 

BN  =  0 
o  o 


(26) 


(27) 


NoNo  =  Jn-m 


(28) 


Let  the  symmetric  matrix  E  in  Eq  ( 19)  be  partitioned  as  follows 


E 


f 

\ 


n  -  m 

n 


(29) 
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Theorem  1: 

A  necessary  and  sufficient  condition  for  P  in  Eq  ( 19)  to  solve  the  ARE  in 
Eq  (9)  is  that 


Q  = 


0  !  A  R_Vf 

o  o 


(30) 


E  = 


flfl 


n  -  m_ 
0 


aA 


(31) 


The  2n  x  2n  pd  solution  matrix  P  is  given  by 


■E 


aA  0 
o  i 

O'Tai. 


n 


(32) 


Proof:  (sufficiency) 


P  = 


T  -11 
aB  B  B  aB 
0  o  .o  o 


aA  |  0 

|__o  | - - 

0  i  ctN  Nr, 
L  i  o  oj 


since,  from  ESqs  (27)  and  (28) 


T 

- 1 

r  , 1  i 

B 

o 

= 

B  btb  )  !n 

Nf 

L 

o 

o' 

, 

o 

o 

1 

L  i  J 

(33) 


•  .% 
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we  get 


P  = 


aA  i  0 
o  I 


0  I  al 


n 


Q  =  a2BR_1BT  -  PA  -  ATP  = 


0  !  0 


-  '  <T  ~  -  _  i  -  t 

0  |  a  \B  R  1B1 


o  o 


(necessity)  from  Eq  (25) 


Since 


Q  =  o^BR'1  -  PA  -  A  P 


P  = 


E,  l 

O  I 
- 1 


eV 

l  O 


|n  E  J  ob(bTb\  +  N  E.N 
L  o  2  I  o\o  o)  o  o  1 


we  find 


Q 


EJN  A  +  A  nJe2 
o  o  o  o  o  2 


A  LbJ  (B  b)_1bJ  +  N  E  N 
_  _o  L  _ o _ o_l_  __ 


1 


Lb  (btb  )  1BT  +  N  E  nT" 
LL  oV  o  ol  o  o  1  oj 


|A0~E* 


a2B  R_1bJ  -  N  E2  -  E  nJ 
o  o  o  2  2  o 


Since  Q  is  block  diagonal,  then  using  Eq  (33)  and  Lemma  1,  we  find 


E 

E 

E 


1 

2 

3 


aIn  -  m 
0 

aA 

o 
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Once  the  solution  to  the  ARE  is  obtained,  one  can  find  the  optimal  control  as 
well.  For  the  case  where  the  sensors  and  actuators  are  co-located  we  have  the 
following. 

Corollary  1; 

When  the  output  feedback  control  law  is  chosen  as 

u  =  Ky  (34) 

where  K  is  a  symmetric,  m  x  m  negative-definite  (nd)  matrix  and  the  output 
is  given  through 


y 


T 

B  x 


(35) 


then 


R  =  -aK  1  (36) 

is  the  m  x  m  control  penalty  matrix  of  the  quadratic  objective  function 

/OO 

^xTQx  +  uTRu|  dt  (37) 

u 

Proof: 

The  optimal  control  that  minimizes  Eq  (37)  is  given  by 

-1  T 

u*  =  -R  BaPx 

using  the  solution  matrix  P  from  Eq  (32)  we  find 

u*  =  -aR_1BTx 
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Appendix  C 

DESIGN  OF  CONTROL  LAWS  WITH  FREQUENCY 
SHAPED  COST  FUNCTIONALS 


The  design  problem  treats  a  linear  problem 


x  =  Fx  +  gu 


y  =  Hx 


and  the  following  performance  index  is  to  minimized: 


OO 

-if 


x*(j w)A(jw)  X(joj)  +  u*(jw)B(ju;)u(jw)  du> 


Note  that  the  weighting  functions  are  functions  of  frequency. 

If  the  weighting  functions  A(jw)  and  B(ju)  are  assumed  to  be  rational 

2 

functions  of  squared  frequency,  •«  ,  a  systematic  control  design  procedure, 
may  be  developed  for  positive  semidefinite  A(jw)  and  positive  definite  B(jw) 
This  is  not  a* serious  limitation  because  a  wide  variety  of  functional  forms  may 
be  approximated  by  ratios  of  polynominals .  To  develop  a  specific  control  de¬ 
sign  procedure,  it  is  further  assumed  that  A(jw)  has  rank  p  ,  and  B(ju>) 
is  positive  definite  with  full  rank,  q 


A(jw)  =  PffjcoJPjfjw) 


B(jw)  =  P|(jw)P2(jw) 


C-l 


(1) 


(2) 


r 


and  P2  are  pxn  and  qxq  rational  matrices.  Define 


P^jotfx  =  xJ 


P2(iw)u  =  u 


If  P  (j oj)  is  a  ratio  of  polynomials  in  ju>  and  the  number  of  zeros  does 
not  exceed  the  number  of  poles,  Eq  (3)  may  be  written  as  a  system  of  differ¬ 
ential  equations  with  output  x1  . 


zi  =  Fizi  +  V 


x  =  HjZj  +  djX 


is  zero  if  the  number  of  poles  is  at  least  one  more  than  the  number  of 
zeros.  Equation  (4)  may  also  be  written  in  terms  of  a  differential  equation, 
again  if  the  number  of  zeros  does  not  exceed  the  number  of  poles. 


Z2  =  F2Z2  +  °2U 


u  =  H2z2  +  °2U 


The  dynamic  Eq  (1)  and  the  cost  functional  (16)  may  now  be  written  in  terms 
of  an  extended  state  vector. 


d  z 

m  zi 


F  0  0 

°i  Fi° 
OOF 


J  =  E  xT  (z 

ss  l 

l  ■?  “T) 

0 

0 

X 

T 

HID1 

0 

0 

Z1 

0 

HX 

T 

H2°2 

Z2 

L° 

D2H2 

T 

D2°2J 

u 

Defining  appropriate  vectors  and  matrices,  Eqs  (7)  and  (8)  become 


(8) 


=  F^1 


G*u 


(9) 


(10) 


The  control  law  is  obtained  by  solving  the  following  modified  algebraic  Riccati 
equation 

-  S  F1  -  FXS  -  A  +  (S  G1  +  N)  B'1  (SG1  +  N)T  =  0  (11) 

and 


u  =  B'1  (S  G  +  N)  V  (12) 

Equation  (12)  is  written  in  an  equivalent  form  as 

u  =  CjX  +  C2z1  +  C3z2  (13) 

The  generalized  controller  structure  is  then  shown  in  Fig.  1.  This  controller 
has  the  form  of  a  general  dynamic  compensator.  The  transfer  function  between 
the  control  input  u  and  the  state  x  is 
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Zj  (]'“)  =  (j«I  - 

z2  (jw)  =  (jwl  - 

C3(jwl  -  F2)_1G2]  u(jw)  =  [Cj 


F^GjX 

F2)_1g2U 

+  C2(jwl  -  F^)_1G j] 


x(jw) 


(14) 

(15) 


(16) 


or 


u(jw)  = 


I  -  c8(j«i  -  f2 


S] 1  [<ci 


+  C2(jwl  -  F 


x(jw) 


The  compensator  may  be  expressed  in  an  equivalent  form  shown  in  Fig.  2. 


Note  that  the  transfer  functions 

between  x1 

and 

u  and 

u1  and  u  are 

x-1(jw)  = 

[Hj  (j«I 

Di] 

x  (jw) 

(17) 

u/dw)  = 

[h2  <1«I 

•vS  + 

D?1 

u  (jw) 

(18) 

Therefore,  the  poles  of 

Pj(ju) 

and  P2(jw) 

show 

up  as  compensator  poles 

and  zeros,  respectively  and  directly  influence  the  closed  loop  transfer  func¬ 
tion.  This  provides  a  direct  relationship  between  the  form  of  the  frequency 
dependence  shaping  used  in  the  cost  function  and  the  structure  of  the  over¬ 
all  compensator. 

If  the  number  of  zeros  in  P^(ju)  exceeds  the  number  of  poles,  Eq  (5)  must 
be  modified.  For  example,  if  Pj(jw)  has  one  more  zero  than  poles,  Eq  (5) 
may  be  written  as 
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Using  Eq  (1),  this  may  be  written  as 


=  PjZj 


Gj  x 


x1  =  +  DjX  +  (Fx  +  Gu) 

=  +  (Dj  +  D^F)  x  +  Dn  Gu  (20) 

Therefore  when  P^joj)  may  have  more  zeros  than  poles,  a  more  general 

form  for  Eq  (5)  is 

z  =  F1z1  +  GjX 

V1 

x1  =  H.z.  +  D.x  +  Di  *  (21) 

i=0 

where  u^  is  the  itk  derivative  of  u  .  Equations  (6),  (7),  and  (8)  are 
also  modified  in  a  similar  way. 

Selection  of  Weighting  Matrices 

It  appears  that  a  solution  can  be  guaranteed  if  B(jw)  is  positive  definite  and 
A(jw)  is  positive  semi-definite  at  all  but  a  finite  number  of  discrete  frequencies 
(though  this  is  not  a  necessary  condition).  It  should  be  pointed  that  that  even 
under  these  constraints,  the  solution  may  not  be  easy  to  find  and  in  fact  may  not 
even  be  causal.  The  total  class  of  weighting  functions  for  which  a  causal  may  be 
found  will  be  subjects  of  future  research. 
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Appendix  D 

MINIMUM-ENERGY  CONTROL 


Consider  the  linear  dynamic  system  described  by 

x(t)  =  A(t)  x(t)  +  B(t)  u(t)  ,  x  Rn  ,  u  Rn  (1) 

Next,  define  the  controllability  matrix  by 

t. 

Wtf^tj)  4  /  1  $(to,T)B(r)B'(r)  $'(to,r)dr  (2) 

where  prime  indicates  transposition  (real  or  complex),  and  $(.)  is  the  transi¬ 
tion  matrix.  Then  we  have  the  following. 

THEOREM 

For  a  controllable  system  described  by  Eq  ( 1) ,  the  controllability  matrix  is  non¬ 
singular  and  the  input  that  transfers  the  initial  state  xQ  to  a  desired  final  state 
x(tj)  in  a  finite  time  tj  >  tQ  is  given  by 

U*(t)  =  B,(t)*(t0,t)W"1(t0.t1)  *(t0,t1)x(t1)-xQ  (3) 

Furthermore,  this  input  requires  the  minimal  control  energy.  This  minimal 
energy  is  given  by 


J*  #  /  1  ||u<r>|J2  dr  =  x'Q  W"1(t0,t1)x0 


for  a  system  whose  final  state  is  the  origin  x(tj)  =  0 


(4) 
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PROOF: 


First  we  prove  that  the  control  in  Eq  (3)  transfers  the  initial  state  xq  to  the 
desired  final  state  x(t^)  .  The  solution  to  Eq  (1)  is  given  by 


x(t)  =  $(t,tQ)xo  +  /  $(t,r)B(T)u(T)d  =  4>(t,tn) 


r 

X  +  J  4>(t  tT)B(T)u(T)dT 


where  we  have  used  the  property 


(5) 


$(t,r)  =  $<t,t0)$(to,T)  . 


Using  Eq  (3)  in  Eq  (5)  we  find 


x( t , )  =  $(t  ,t  ) 
1  1  o 


- 

X  + 

/  1  *(t  ,T)B(T)B’(T)$(t  ,T)dT 

o 

r  o  o 

_o 

Using  the  definition  in  Eq  (2) 


and  using  the  fact  that 


*«rV 


we  obtain  the  result. 

Next  we  show  the  minimal  energy  property.  Let  {u(t)}  be  any  control  sequence 
that  transfers  the  initial  state  xQ  to  x(tj)  ,  then  we  have  to  show  that 


t 


/ 


o 


1 


||u(t)||2  dt 


> 


||u*(t)||2  dt 


(6) 


Define 


x  =  $ 


-1 


(tl>to)x(tl)  " 


-  *(tQ,t1)x(t1) 


then,  by  the  assumption  that  both  control  sequences  transfer  xq  to  x(t^)  we 
have  from  Eq  (5) 


x 


♦(t^TlBOOujOOd 


/  <l>(to,T)B(T)U*(T)d 

to 


Subtracting 


o 


0 


D-3 


which  implies 


/  A 


(  f  *(t0»r)B(r)  Fu1(t)-u*(t)1  dT,  W  1(tQ,t1) 
Vo  / 


x)  =  0 


where  (. ,  ^  is  the  usual  inner  product . 
Since  (x,Aj^  =  (a'x,$  we  have 


(t)— u*(t)  ,  B'(T)0'(to>T)W  1(t0,t1)lN 


dT  =  0 


and  using  Eq  (3)  this  becomes 


f  <^1<T>"U*(T>>  U*(T^ 


dT  =  0 


Consider  now 


tj  tx 

/  Hui'T>n2d  -  /  ||ui(t)  -  U*(t)  +  U*(T)||2dT 


tl  tl 

=  f  ||u1(t)  -  u*(t)||2  d  +  /  ||u*(t)||2  dT 
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**/  4  j(t)  -  u*(T) ,  u*(t^  dr 

t_ 


h  tl 

=  f  llu^r)  -  u*(t)||2  dT  +  J  ||u*(t)||2  dr 


Since 


u*(r)||2  dr  >  0 


we  conclude 


Finally 

J *  ~  J  I|u*(t)||2  dT 

t 

o 


=  J  x'  w'1(to>t1)i'(to,T)B(T)B,(r)i(to,T)W'1(to,t1)x  dr 
4o 

=  x’  W-1  x  =  xQ  W-1  xQ 


n 
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6  MISSION  b 

$  °f  5 

§  Rome  Air  Development  Center  ? 

n  RAPC  plan*  and  execute*  research,  development,  tut  and  2 

*  selected  acquisition  programs  In  support  o{  Command,  Control  2 

k  Communications  and  Intelligence  IC*D  activities.  Technical  f 
%  and  enginee/iing  support  loithln  areas  oi  technical  competence  C 
|  cs  provided  to  ESV  Program  0 Hlc.es  IPOs)  and  other  ESP  $ 

D  elements.  The  principal  technical  mission  areas  axe  V 

k  communications,  electromagnetic  guidance  and  contxol,  Sux-  'h 
*s  velllance  oi  ground  and  aerospace  objects,  intelligence  data  i 

*  collection  and  handling,  iniormation  system  technology,  2 

&  ionospheric  propagation,  solid  state  sciences,  mlcromve  ? 
%  physics  and  electronic  reliability,  maintainability  and  ? 
S  compatibility.  ® 
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